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16362 CENTRIFUGE TESTS AND OVERCONSOLIDATED SLOPES 


KEY WORDS: Centrifugation; Clays; Failure; Overconsolidated 
Overconsolidation; Plasticity; Shear strength; Slope stability; Soil mechanics 


ABSTRACT: method presented for predicting the shear strength mobilized during 
undrained failure model slopes overconsolidated clay. Also long-term failure 
model slopes overconsolidated clay described. drum centrifuge was used 
conduct these model studies both short-term and long-term stability. produce 
overconsolidated conditions, kaolin slurry was consolidated the centrifuge 
acceleration between 150 and 420 After consolidation, the clay deposit was 
allowed rebound and 2-in. (50.8-mm) high, 60° slope was cut around the entire 
circumference the drum. Short-term failure was produced rapidly increasing the 
speed the centrifuge until failure occurred. Long-term stability was investigated 
maintaining acceleration lower than that required cause undrained failure. 
five long-term tests failure occurred after hours continous operation. 


REFERENCE: Fragaszy, Richard (Asst. Prof. Civ. Engrg., San Diego State Univ., 
San Diego, Calif. 92182), and Cheney, James A., “Drum Centrifuge Studies Over 
Consolidated Slopes,” Journal the Geotechnical Engineering Division, ASCE, Vol. 
107, No. GT7, Proc. Paper 16362, July, 1981, pp. 843-858 


16377 MIAMI LIMESTONE FOUNDATION 


KEY WORDS: Construction practices; Dewatering; Florida; Foundation 
construction; Foundation design; Limestone; Porosity; Rock analysis; Rock 
mechanics; Rock tests; Shallow foundations; Strength; Thickness 


ABSTRACT: The techniques and peculiarities are reviewed regarding the design and 
construction shallow foundations bearing the Miami Limestone. The thickness 
and lateral extent the Miami Limestone are illustrated. The geologic characteristics 
this formation are presented, and the types and properties Miami Limestone are 
explored. Methods analysis for shallow foundation design are described, and 
construction problems and difficulties are considered with case histories. Several 
criteria are suggested for evaluating shallow foundations bearing the Miami 
Limestone. Because this limestone variable strength, porosity and thickness, 
adequate investigation these properties should conducted prior foundation 
construction each particular site. 


REFERENCE: Kaderabek, Thomas (Sr. Engr., Law Engrg. Testing Co., 200 San 
Lorenzo Ave., Coral Gables, Fla. 33146), and Reynolds, Richard T., “Miami 
Limestone Foundation Design and Construction,” Journal the Geotechnical 
Division, ASCE, Vol. 107, No. GT7, Proc. Paper 16377, July, 1981, pp. 
859-872 


16374 SHEAR STRENGTH ROCKFILL 


KEY WORDS: Compaction; Estimates; Failure; Friction; Particle shape; 
Porosity; Rockfill dam design; Roughness; Shear strength; Shear tests; 
Strains; Triaxial tests; Uniaxial compressive strength 


ABSTRACT: practical method for estimating the shear strength rockfill 
developed. The peak drained friction angle found closely related 
that rock joints. both cases the values are dependent sample size, 
stress level, surface roughness, and the uniaxial compression strength the rock. 
Friction angles are therefore higher for smaller samples, and very high where stresses 
are low, the toe near the face rockfill dam. Test data reviewed shows that 
the value for rockfill can quantified equivalent roughness (R), and 
equivalent particle strength (S). The value (R) depends the porosity following 
compaction, and the degree particle roundedness and surface smoothness. 
practical method proposed for physically measuring the full-scale shear strength 
in-place compacted 


REFERENCE: Barton, Nick (Sr. Staff Consultant, Terra Tek, 420 Wakara Salt 
Lake City, Utah 84108; formerly Sr. Dam and Rock Group, 
Geotechnical Inst., Oslo, Norway), and Kjaernsli, Bjorn, “Shear Strength Rockfill,” 
Journal the Geotechnical Engineering Division, ASCE, Vol. 107, No. GT7, Proc. 
Paper 16374. July, 1981, pp. 873-891 
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16391 CLAY ANISOTROPY AND BRACED WALL BEHAVIOR 


KEY WORDS: Anisotropic soils; Anisotropy; Braced excavation; Braces; 

Clays; Earth pressure; Excavation; Finite element method; Heaving; 
Movement; Shear strength (soils); Soil stability; Walls 

ABSTRACT: technique presented which modifies the basal heave analysis 
procedures for braced excavations account for the strength anisotropy exhibited 
most natural, soft medium clays. the case braced wall, the clay the 
passive side the wall can develop much lower strength than the clay the active 
side, the clay strongly anisotropic. Since the existence lower strength zone 
not considered conventional basal heave analysis procedures, these techniques 
theoretically lead overestimate basal stability. Finite element analyses are used 
demonstrate the effects anisotropy the behavior braced excavations. 
Analysis results show that the soil presumed isotropic when truly 
anisotropic the following trends develop: (1) The basal heave safety factor may from 
percent larger than the actual one, depending the degree anisotropy; (2) 
lateral wall movements and ground surface settlements can larger than expected; 
and (3) unusually high loads can develop the lower struts cases marginal 
stability. 

REFERENCE: Clough, Wayne (Prof. Civ. Engrg., Stanford Univ., Stanford, 
Calif. 94305), and Hansen, Lawrence A., “Clay Anisotropy and Braced Wall 
Behavior,” Journal the Geotechnical Engineering Division, ASCE, Vol. 107, No. 
GT7, Proc. Paper 16391, July, 1981, pp. 893-913 


16367 FOUNDATIONS AND SAND LAYERS 


KEY WORDS: Bearing capacity; Circular footings; Earth pressure; 
Foundation investigations; Layered soils; Model tests; Sand; Shallow 
foundations; Soil mechanics; Subsoil 


ABSTRACT: The ultimate bearing capacity footings subsoil consisting 
strong sand layer overlying weak sand deposit was investigated. theory was 
developed and then results predicted. These predictions compared well with test data 
model strip and circular footings resting dense sand overlying ioose, well as, 
compact sand. Design charts are presented for use practical application. Based 
these charts, procedure determine the ultimate bearing capacity strip and 
circular footings strong sand layer overlying weak sand deposit described. 


REFERENCE: Hanna, Adel (Assoc. Prof., Concordia Univ., Dept. Civ. Engrg., 
1455 Maisonneuve Boulevard West, Montreal, Quebec, Canada H3G 
“Foundations Strong Sand Overlying Weak Sand,” Journal the Geotechnical 
Division, ASCE, Vol. 107, No. GT7, Proc. Paper 16367, July, 1981, pp. 
915-92 


16369 SLIP CIRCLES AND CRITICAL SHEAR PLANES 


KEY WORDS: Mohr circle; Safety factor; Shearing; Shear planes; Slope 
stability; Slope stability analysis; Soil mechanics; Soil stability; Stability 


ABSTRACT: Assuming linear stress distribution, the normal effective stress 
determined the bottom vertical interslice boundaries along slip circles through 
earth slope. circle stress constructed for each point, and the slope 
the critical shear planes and the factor safety these planes are determined. These 
values are compared, respectively, with the slope the slip circle and the factor 
safety against failure the slip surface. The correlation between the values found 
poor for arbitrarily chosen slip circles but good for critical circles, except the 
upper end where the soil tension. 


REFERENCE: Spencer, Eric (Principal Lect., Dept. Civ. Engrg., Sheffield City 
Polytechnic, Pond Street, Sheffield, 1WB, England), “Slip Circles and Critical Shear 
Planes,” Journal the Geotechnical Engineering Division, ASCE, Vol. 107, No. GT7, 
Proc. Paper 16369, July, 1981, pp. 929-942 


16366 DYNAMIC PROPERTIES SOILS 


KEY WORDS: Cohesionless soils; Cyclic loads; Dynamic properties; Field 
data; Field tests; Rayleigh waves; Soil dynamics; Soil-structure interaction; 
Soil tests; Stiffness tests; Test procedures; Test results; Vibration tests 


ABSTRACT: Various field test (namely vibration tests blocks plates, steady-state 
vibration Rayleigh wave tests, wave propagation tests, and cyclic load tests) were 
conducted number sites India determine the dynamic shear modulus, 
Data obtained different sites are described. The values obtained from the 
different tests given site vary widely. The rational approach for selecting the value 
from field tests for use the analysis and design soil-structure interaction 
problems under dynamic loads must account for the factors affecting The suggested 
approach, which provides possible answer, suitable cohesionless soils below the 
water table where rather difficult, not impossible, obtain undisturbed samples. 


REFERENCE: Prakash, Shamsher (Prof., Civ. Engrg. Dept., Univ. Missouri-Rolla, 
Rolla, Mo. 65401; formerly Prof. Civ. Engrg. and Head Geotechnical Engrg. 
Div., Univ. Roorkee, Roorkee, UP, India), and Puri, Vijay Kumar, “Dynamic 
Properties Soils from In-Situ Tests,” Journal the Geotechnical Engineering 
Division, ASCE, Vol. 107, No. GT7, Proc. Paper 16366, July, 1981, pp. 943-963 


16379 NEW DESIGN CORRELATIONS FOR PILES SAND 


KEY WORDS: Bearing capacity; Cohesionless soils; Design criteria; 
Friction; Pile driving; Pile-driving formulas; Pile foundation design; Pile 
foundations; Residual stress; Residual stress measurement; Sand; Soil density 


ABSTRACT: critical review the literature clearly shows that existing theories 
which are used determine the bearing capacity piles driven sand are not 
satisfactory. Careful examination the results experimental studies indicates that 
existing theories fail consider all significant parameters. Field load test data are used 
determine which the pile geometry and soil property parameters are significant. 


Special attention given the problem residual stresses. Bearing capacity factors 
are calculated from measured unit point and side resistances. The bearing capacity 
factors are correlated with relative depth (depth penetration diameter ratio) and 
friction angle relative density sand. New design correlations are developed for 
use predicting the bearing capacity axially loaded piles sand. 


REFERENCE: Coyle, Harry (Prof. Civ. Engrg., Texas Univ., College 
Station, Tex. 77843), and Castello, Reno R., “New Design Correlations for Piles 
Sand,” Journal the Geotechnical Engineering Division, ASCE, Vol. 107, No. GT7, 
Proc. Paper 16379, July, 1981, pp. 965-986 


U.S. Conversion Factors 


accordance with the October, 1970 action the ASCE Board Direction, which stated 
that all publications the Society should list all measurements both U.S. Customary and 
(International System) units, the following list contains conversion factors enable readers 
compute the unit values measurements. complete guide the system and its 
use has been published the American Society for Testing and Materials. Copies this 
publication (ASTM E-380) can purchased from ASCE price $3.00 each; orders must 
prepaid. 

All authors Journal papers are being asked prepare their papers this dual-unit format. 
provide preliminary assistance authors, the following list factors and guides 
are recommended the ASCE Committee Metrication. 


Multiply 

convert 
inches (in.) millimeters (mm) 25.4 
feet (ft) meters (m) 0.305 
yards (yd) meters (m) 0.914 
miles (miles) kilometers (km) 1.61 
square inches (sq in.) square millimeters 645 
square feet (sq ft) square meters 0.093 
square yards (sq yd) meters (m’) 0.836 
square miles (sq miles) square kilometers 2.59 
acres (acre) hectares (ha) 0.405 
cubic inches (cu in.) cubic millimeters 16,400 
cubic feet (cu ft) cubic meters 0.028 
cubic yards (cu yd) cubic meters 0.765 
pounds (Ib) mass kilograms (kg) 0.453 
tons (ton) mass kilograms (kg) 907 
pound force newtons (N) 4.45 
kilogram force (kgf) newtons (N) 9.81 
pounds per square foot (psf) pascals (Pa) 47.9 
pounds per square inch (psi) kilopascals (kPa) 6.89 
U.S. gallons (gal) liters (L) 3.79 
acre-feet (acre-ft) cubic meters 1,233 
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DruM CENTRIFUGE STUDIES OVERCONSOLIDATED 


Richard Fragaszy,' ASCE and James Cheney,’ ASCE 


INTRODUCTION 


Civil engineers have been concerned with the stability overconsolidated 
clay slopes for over century. Because overconsolidated clay deposits are 
common geological features much the industrialized world, slope failures 
these clays can cause considerable damage both economic and human 
terms. Despite long history research the behavior overconsolidated 
clay, there still uncertainty the mechanisms producing slope failures. 
Neither the maximum strength which can mobilized failure nor the delay 
between slope formation and failure are well-defined. 

This paper describes series centrifuge tests designed model both 
short-term and long-term stability overconsolidated clay slopes. drum 
centrifuge, type machine used only one previous study slope stability, 
was used test 2-in. high, 60° model slopes. The objectives this 
study were develop method predict the strength mobilized during undrained 
failure overconsolidated slopes and determine long-term failure model 
slopes would occur centrifuge. 


Review Previous Work 


The literature long-term failures overconsolidated clay slopes covers 
period over 100 yr. The oldest papers, such Refs. and were mainly 
descriptive; however since the publication Ref. there have been many 
theories explaining this phenomenon. These theories have been based several 
types analysis: (1) Total stress with undrained shear strength and 
(2,20); (2) effective stress analysis with (8,9,11,21,22); (3) residual strength 
(3,4,17,23); (4) creep (14,15); and (5) critical state (18). extensive literature 
review has been presented elsewhere (6); however the theories based creep 

the October 22-26, 1979, ASCE Annual Convention and Exposition, held 
Atlanta, Ga. 

Prof. Civ. Engrg., San Diego State Univ., San Diego, Calif. 92182. 

Civ. Engrg., Univ. California, Davis, Calif. 95616. 

Note.—Discussion open until December 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication February 12, 
1980. This paper part the Journal the Geotechnical Engineering Division, Proceedings 
the American Society Civil Engineers, ©ASCE, Vol. 107, No. GT7, July, 1981. 
ISSN 0093-6405 /0007-0843 /$01.00. 
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and the critical state concept are briefly examined the following. 

1968, Schofield and Wroth (18) proposed thoery concerning the failure 
mechanism overconsolidated clay slopes based the Cam-Clay theory. They 
hypothesized that after peak strength reached the failure zone, dilation 
occurs very thin zone (perhaps millimeter thick) the expense soil 
either side the failure zone. The soil this zone becomes wetter and 
loses strength until the critical state reached. this point, the strength 
the soil reduced the which are numerically 
equal the peak strength parameters normally-consolidated sample 
the soil. Two years after this theory was proposed, Skempton (25) stated that 
the basis field observations the limiting strength first-time slide 
overconsolidated clay the fully-softened strength, not the residual strength 
had previously stated (23). 

Acreep theory failure has been recently proposed Nelson and Thompson 
(14,15). They have postulated that fully-drained conditions exist the slope 
and, under sustained shear loading, creep occurs which gradually reduces the 
strength the clay its residual values. Taking data from Skempton (23) 
they were able find creep parameters for their model which resulted calculated 
times failure within the actual times failure. 

The short-term stability overconsolidated clay slopes has not been subjected 
the same scrutiny long-term stability. and LaRochelle (24) describe 
short-term failure deep excavation London Clay. Since the slide occurred 
only days after the excavation was completed, undrained analysis was 
performed. The results the analysis indicated that only 55% the calculated 
strength was mobilized. The authors attribute this discrepancy the migration 
pore water and the presence fissures. 

Lyndon and Schofield (12) reported the centrifuge modeling prototype 
trench (33 ft) deep with side slopes and 1/2:1. failure occurred 
after min the centrifuge, corresponding weeks the prototype 
scale. clearly-defined slip surface was found and stability analysis was 
carried out. The shear strength mobilized was approx 60% the average measured 
shear strength prior the excavation. closely matched the shear strength 
the soil (determined immediately after failure) behind the failure wedge. 
increase water content near the failure surface was also observed. 

series centrifuge tests overconsolidated boulder clay was carried 
out English and Schofield (5) and reported Schofield (19). Use 
basket centrifuge allowed the modeling very long slopes. Although their machine 
was capable maximum acceleration 1,000 the slope height was restricted 
ft) maximum acceleration. Their results showed that the use undrained 


cohesion short-term stability analyses unconservative for overconsolidated 
deposits. 


AND TESTING 


Previously, virtually all centrifuge testing has been conducted with bucket 
centrifuges which the model placed bucket the end one 
the centrifuge arms. The other arm contains either counterweight identical 
bucket and model. major difficulty with bucket centrifuge testing slopes 
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the limited amount soil which may placed the bucket. Although 
the depth the soil may sufficient model very high slope, the length 
the slope severely limited, resulting typical length height ratios 
six less. When the length height ratio this low, friction against the 
sides the bucket will cause change the stress distribution within the 
slope and may actually the slope. Reducing the scale the model 
will lessen this problem, but the. maximum speed the centrifuge will place 
upper limit the scaling factor used. addition, construction the model 
becomes more difficult the model size reduced and the deformation 
measurements become less accurate (13). 

method centrifuge modeling which particularly advantageous for slope 
stability studies has been proposed Schofield (19). 
centrifuge employed which soil can placed around the entire circumference 
form deposit. slope cut such centrifuge would have 
rigid boundary the circumferential direction. The potential this type 


FIG. Centrifuge 


machine was demonstrated English and Schofield (5) the study 
the aforementioned short-term stability overconsolidated clay slopes. 

The centrifuge used for the research reported here, shown Fig. was 
constructed modifying existing 4-ft (1.2-m) diam bucket centrifuge. 
steel drum in. (230 mm) wide, in. (130 mm) deep and (1.2 diameter 
was attached the centrifuge arms. its maximum speed 1,000 rpm (650 
g), the centrifuge can carry approx 450 (200 kg) soil, capacity excess 
145 g-tons. 

form overconsolidated deposit, clay slurry poured inio the drum 
while spinning approx 120 rpm. After the drum has been filled, the 
speed increased produce radial acceleration between 150 and 450 
porous blanket the bottom the drum allows the clay drain both 
the top and the bottom the deposit. After primary consolidation, the water 
the surface the clay flows down through the clay causing further 
consolidation. 

When there standing water the surface, the centrifuge turned 
off and allowed stop. Two samples are then removed determine the variation 
water content with depth. Water content data are subsequently used with 
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laboratory consolidation results determine the variation maximum effective 
vertical consolidation stress, with depth. 

After water content samples are removed, 60° slope cut around the 
circumference the drum. minimize soil disturbance, stiff wire loop 
first used roughly cut out the slope. form the final slope geometry, 
cutting blade raised into the slope while the drum slowly spinning. 

For short-term testing, the centrifuge rapidly accelerated until failure 
observed (less than 2.5 min all cases). stroboscope connected 
photoelectric pick-off and flash delay used observe any location around 
the circumference the drum. complete sweep around the circumference 
takes about sec. 

For long-term stability tests, the centrifuge rapidly accelerated level 
between 40% and 60% the level required cause quick failure. The speed 
then maintained and the slope observed until failure occurs. prevent 
the clay from drying out, fine mist water sprayed into the drum 
frequent intervals. 

Two series seven short-term tests were conducted (Q2 and series). 
For each series the consolidation acceleration was varied between 150 and 
450 50-g intervals. The only difference between the two series was that 
the series the top 1/2 in. (13 mm) kaolin was removed prior cutting 
the slope; this was not done the series. 

Only five out ten long-term tests were successful. the remaining tests, 
fluctuations the centrifuge speed invalidated the results. 


Properties 


The clay used this study mixture two commercial grades kaolin 
manufactured California, i.e., Snow-Cal and Mono-90. Mono-90 the finer 
the two with all particles finer than No. 325 sieve. All Snow-Cal particles 
are finer than No. 250 sieve. For all testing, three parts Snow-Cal were 
mixed with one part Mono-90. The powder was then mixed with distilled 
water form slurry. 

Atterberg limits tests were performed the kaolin. The results these 
with calcium the exchangable ion. very low plasticity clay which 
barely falls into the classification the Casagrande Chart. 

The average two oedometer tests with maximum stress kg/cm? 
(230 psi) was used determine the Cam-Clay parameters (18) 0.086 and 
0.01. Although Cam-Clay theory described terms isotropic 
compression, possible show that oedometer test produces virgin 
compression line with the same slope, isotropic test; only the intercept 
given pressure different. The compression index, uniquely related 
and equal 0.198. Also determined were the coefficient 
consolidation, c,, and the average permeability, These values are 1.0 
respectively. 

Drained triaxial tests were performed determine both the Mohr-Coulomb 
strength parameter and the parameter, The values these 
parameters are 31° and 1.24, respectively. 


DRUM CENTRIFUGE STUDIES 
Description 


Under stroboscopic illumination, the failure each slope could easily 
observed. The failure mode observed was quite similar for all tests. The first 
sign failure was the appearance crack approx cm-2 (0.4 in.-0.8 
in.) behind the top the slope. Almost immediately after the appearance 
the crack, failure, defined observable movement the slope, began. The 
failures were neither sudden nor catastrophic. Movement was slow and vertical 
displacements were usually 1/4 1/3 the slope height after minute 
movement. Failure always occurred through the toe; movement occurred 
front the toe only after substantial slope displacement. The failure surface 
extended vertically downward from the top the slope and intersected the 
toe nearly horizontal line. The failure surface was not fully visible from 
top toe, but likely that the surface formed curved arc between them. 


FIG. 2.—Short-Term Failure 


Fig. shows failed slope, which several important features common 
all test failures can seen. The long length embankment failure clearly 
shown, the scarp formation the top the slope. The dark line 
front and parallel the slope was formed during slope excavation and was 
not caused the failure. The failure appears have taken place 
aforementioned block movement. 

Failures frequently occurred almost simultaneously several locations around 
the circumference the drum. The circumferential length these failure zones 
varied between approx and 200 (12 in. and in.) compared 
total circumference 380 (150 in.). When single failure was observed, 
acceleration the drum was stopped and constant speed was maintained 
until additional failures were observed. These additional failures occurred within 
minute two; occasionally secondary failures were also observed. Cusping 
failures, reported English and Schofield (5), were not observed any 
the tests. Because the failure zones were always much longer than the 2-in. 
(50-mm) slope height, plane strain analyses are appropriate. 

The slopes consolidated low accelerations had two distinguishing features. 


— 
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First, they appeared flow somewhat failure rather than fail solid 
block. The flow was not very pronounced, but was clearly observable. Second, 
these slopes failed the entire circumference almost simultaneously. 


INTERPRETATION SHORT-TERM RESULTS 


Before explanation for short-term behavior can suggested, some 
edge the stress history the slope necessary. Fig. shows plot 
the stress history typical element soil the slope. The stresses are 
given terms shear stress, and mean normal effective stress, 
and minor principal effective stresses acting the element, respectively. The 
initial consolidation the kaolin represented the line AB. During this 
portion the test, the vertical (radial) stress the major principal stress. 


/ 
State Line 


Shear Stress q 


Mean Normal Effective Stress, p 


Zero Vertical Stress Lines 


FIG. 3.—Typical Yield Surface Showing Effective Stress History and Loading Path 
for Undrained Tests 


The intermediate and minor principal stresses act the circumferential direction, 
o,, and the direction normal the plane the drum (horizontal direction), 
o,. For the condition zero lateral strain, these two stresses will both 
equal the vertical stress multiplied the coefficient earth pressure 
rest, the centrifuge, this condition only approximately met. the 
circumferential direction some strain inevitable because, the clay consoli- 
dates, moves radially outward new position where the circumference 
greater. Although the side walls the drum are quite rigid, they are inclined 
slight angle inward from the vertical. Therefore, the clay consolidates 
and moves radially outward, small positive horizontal strain induced. The 
actual circumferential and horizontal strains, however, are much smaller than 
the strains which would result the clay were allowed expand freely, i.e., 
were not constrained the drum. The partial constraint against lateral 
movement causes compressive stresses the circumferential and horizontal 
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directions. Although the actual coefficient lateral stress could not deter- 
mined, must bounded the coefficient earth pressure rest, K,, 
and the coefficient active earth pressure, K,. Assuming conditions, the 
therefore the stress path will inclined angle 39° from the p’-axis, 
following examination, conditions will assumed because felt that 
the actual stress state closer conditions than conditions. 

When consolidation complete, the centrifuge stopped, causing one-dimen- 
sional unloading the clay. The stress path induced this unloading can 
approximated using the method proposed Wroth (28). Wroth has shown 
that the unloading path becomes straight line the stress ratio 
plotted against the logarithm p’. The equation this line 


which the gradient the rebound line; and the values 
mean normal effective stress and stress ratio the end consolidation; and 
and the values mean normal effective stress and stress ratio, 
during rebound. Using Wroth’s graph function the plasticity index, 
the value determined 1.45. Assuming conditions during 


This equation can used find two points particular interest the 
rebound line. The first point Fig. where the shear stress, 

The second point interest the limiting value p’, shown Fig. 
point assumed that the limiting value occurs when either 
zero vertical stress reached when the unloading path intersects the yield 
surface, causing reverse yielding, whichever comes first. When the circumferential 
and horizontal stresses are equal, the condition zero vertical stress described 
line g-p’ space with slope —1.5. was found, however, that 
reverse yielding does not occur, because the rebound line does not intersect 
the yield surface. The limiting value can, therefore, found setting 
condition. 

Point represents the state stress the soil after the centrifuge has 
come complete stop and pore pressure equilibration has occurred. The 
final change stress which occurs before the centrifuge restarted caused 
cutting the 60° slope. The cutting allows the soil expand horizontally 
resulting the loss almost all horizontal stress. This change stress 
shown Fig. the line DE. assumed that the circumferential stress 
not changed and the movement space horizontal the line representing 
zero vertical stress when one horizontal stress also zero. The slope this 
the end consolidation and the beginning the stability test was always 
longer than much longer than the time required for pore pressure equilibration. 
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estimate the shear strength the clay undrained loading, the Cam-Clay 
theory used. Undrained loading requires that the stress path space 
between point and the yield surface vertical. determine the shape 
the yield surface, the general equation for yield surface is, i.e. 


Mp’ 

Although the yield surface for each element soil will different, the ratio 
1.0, will 0.8. Substituting these values and setting 1.24 
yields constant equal 0.65. then set equal zero, the value 
1.92. 


Using this information, the yield surface for any element soil the slope 
can obtained from 


which found from the water content data; and 0.48. setting 
0.0 and the constant Eq. can determined and the yield 
surface can drawn, shown Fig. 

major assumption, that the vertical effective load path will extend the 
yield surface (point F), made this point the analysis. Point lies 
region space where one the major principal stresses must tensile, 
3p’; therefore order reach point the yield surface, some 
tensile strength required. this analysis, the clay assumed possess 
the required tensile strength. 

The value the intersection the load path with the yield surface 
(point can found from 


constant and equal 108 Because point the left 
side the critical state line, any small amount drainage this point would 
sufficient cause unstable yielding. Although the tests were run quickly, 
some drainage inevitable, thereby causing dilation and loss strength when 
the yield surface reached. The stress path from point with local seepage 
would progress toward the critical state line the yield surface decreasing 
values shear stress; thus because the applied shear stress cannot reduced, 
the slope becomes catastrophically unstable. 

The maximum shear strength which can mobilized any point the 
slope, 1/2 can now calculated from the variation with depth. 
The equation for shear strength, is: 
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TABLE 1.—Predicted Safety Factors using Cam-Clay Theory 


Consolidation, Failure Predicted 
acceleration, acceleration, safety 


Data Used Regression Analysis 
© Data Not Used in Regression Analysis 


Equation Least-Squares Fit: 
N,= 8.8 18.5 


Average (kg/cm?) 


FIG. 4.—Graph Average Consolidation Stress Versus Level Failure for Short- 
Term Tests 
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Q34 300 0.8 
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confirm the aforementioned hypothesis, slope stability analyses using the 
Modified Bishop method analysis were performed for each the short-term 
tests. The results these analyses, presented Table strongly support 
the proposed failure mechanism. The range predicted safety factors 
and the average the safety factors approx 1.0. plot acceleration 
failure versus average effective consolidation stress shown Fig. The 
consolidation stress was determined from the water content data. least-square 
regression analysis performed these data indicate that the points fit the linear 
equation: 


which the level failure. The three tests with the largest deviation 
from predicted safety factor 1.0 were excluded from the regression analysis. 
These points are shown the open circles Fig. The regression equation 
fits the data very well, indicated the highly significant Pearson product- 
moment correlation coefficient obtained (r,, +0.96, and 0.01). 

Because the entire slope could not viewed simultaneously, the reported 
level failure must considered upper bound. Based the time 
takes make complete sweep the circumference the drum, and 
the angular acceleration the centrifuge, estimated that the actual levels 
failure were less than 90% the reported values. This error leads 
underestimation the predicted safety factor perhaps 5%, and can 
least partially account for the nonzero intercept the regression equation with 
zero effective consolidation stress. 


TABLE 2.—Summary Long-Term Test Results 


Average con- 


solidation 
stress, 
kilograms 
per square 
Con- centimeter 
solidation Failure (pounds Time Predicted 
acceleration, acceleration,| per square failure, safety 


(2) (3) (5) (6) 


“Shear strength calculated from Eq. 


Test 
(1) 
450 5.3 4.5 1.5 
(75) 
QL2 450 5.5 1.0 2.2 
(78) 
QL3 450 6.1 8.3 2.4 
(86) 
450 6.0 7.2 1.7 
(85) 
450 3.2 2.1 
(69) 
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The major finding from the long-term tests that long-term failures occur 
centrifuged model slopes. five tests, failure occurred after h-8 
operation. these tests, the centrifuge speed was closely monitored and 
major fluctuation speed was noted. The surface the slopes after failure 


FIG. 6.—Long-Term Failure with Cracks Slope Face Extending Behind the Slope 
(Test QL2) 


was moist and felt similar the slope surfaces after short-term failures. 
summary the data from long-term tests presented Table The failures 
which occurred three the five successful long-term tests were quite similar 
appearance failures short-term tests. each these tests, failures 
occurred one more locations the drum, and each was about 


FIG. 5.—Long-Term Failure (Test QL5) 
| 
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The failures produced the other two long-term tests were quite different 
from short-term failures. these tests, the first sign failure was horizontal 
cracks which appeared the face the slope rather than behind the slope. 
Failures eventually occurred, beginning these cracks. Some the failures 
did not extend behind the slope but were contained the slope face itself. 
Other failures extended the top the slope shown Fig. 

The average vertical consolidation stress for each the long-term failures, 
plotted versus the level failure, shown Fig. along with corresponding 
data for short-term failures. anticipated, the long-term data not match 
the trend evident the short-term data. The level failure for each long-term 
test lower than would expected based the analysis short-term data. 
determine the difference was statistically significant, dependent student’s 
statistic was calculated the actual level failure versus the level 


Fit Short-Term Data 


FIG. Average Consolidation Stress Versus Level Failure for Long- 
Term Tests 


predicted the linear regression equation fitted data from short-term tests. 
The difference between predicted and actual values was significant (dependent 


INTERPRETATION 


The stress history (from point point Fig. the model slopes 
used long-term stability experiments identical the preceding for short-term 
tests. The difference the two types tests the effective stress load path 
(point point Fig. 3). the long-term tests, the load path vertical 
only until the desired acceleration reached. Pore pressure increases during 
the initial minute two the test while the centrifuge accelerating because, 
although the total stresses increase the soil, the mean normal effective stress 
remains constant, and therefore, the pore pressure rises. The excess pore pressure 
will dissipate with time, causing increase. The increase important 
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difference between model slopes and the prototype slopes. the model, the 
soil equilibrium after the slope cut, just prior the acceleration 
failure. The acceleration causes increase effective stress the slope. 
the prototype, cutting the slope unloading process (although shearing 
increase). the field, negative pore pressures are induced when 
cutting made. These negative pore pressures dissipate over time causing 
reduction p’. The reduction has been cited the cause delayed 
failure (27). Although this phenomenon may contribute the failure prototype 
slopes, can rejected the failure mechanism the centrifuge tests. 

Creep second possible failure mechanism, aforementioned. Laboratory 
creep tests were not performed the soil used this study; however Paduana 
(16) has reported relationship between creep rate and plasticity index for 
wide variety clays, including kaolinite. His findings clearly indicate that 
kaolinite exhibits very low creep rate. The low creep rate kaolinite and 
the equivalent creep rates the centrifuge model and the prototype slope 
both make unlikely that creep-induced failures occurred these centrifuge 
model tests. 

The third and most promising explanation failure the critical state theory. 
zones local overstress, such the toe the slope, possible for 
the state stress reach the yield surface. When this occurs water sucked 
into the soil and the soil loses strength. The rate which this dilation occurs 
function the permeability the soil. The permeability the soil used 
these experiments quite high. Consequently, the process dilation 
occurred simultaneously over the entire failure arc, failure would occur quickly, 
perhaps min less. progressive failure, however, would take much 
longer. first, only small fraction the potential failure arc subjected 
shearing stresses high enough cause plastic yielding. these regions 
yield and lose strength, the shearing stress the remaining portion the 
failure arc increases, causing additional soil yield. This process will continue 
until either the slope fails the shear strength has been reduced the 
fully-softened strength along the entire potential failure arc. 

reasonable expect that the preceding gradual reduction strength 
would take several times longer than the time required for 90% consolidation 
the soil, due the progressive nature the failure mechanism. For 
the model slopes, this would mean time failure several hours 
min). Also, because the process governed least partially the flow 
water through the soil, the centrifuge tests should model prototype events 
which take the order times longer occur when the level 
during the test. This because both the distances and the pore pressure gradients 
involved are scaled the factor the centrifuge. 

The long-term stability experiments did not provide sufficient data prove 
disprove any the aforementioned hypotheses. However, they show 
that long-term stability can modeled the centrifuge. Future centrifuge 
experiments which pore pressures the slope are measured should provide 


the quantitative data necessary for more rigorous analysis long-term failure 
mechanisms. 


Summary 


The objectives set forth the beginning this research have been met. 
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method predict the shear strength mobilized during undrained failure 
overconsolidated model slopes has been postulated. Second, long-term failures 
were observed the model slopes, indicating the potential drum centrifuge 
studies long-term stability problems. 

The prediction shear strength mobilized failure accomplished three 
steps. First, the variation water content with depth used determine 
the variation vertical effective consolidation stress clay deposit consolidated 
the drum centrifuge. Second, the effective stress unloading path caused 
stopping the centrifuge and cutting the slope postulated. Third, the intersection 
vertical effective stress load path (caused accelerating the drum until 
slope failure observed) with the Cam-Clay yield surface determined. The 
shear strength the intersection the predicted strength failure. Using 
this method strength prediction, slope stability analyses were performed. 
The predicted safety factors the time failure for tests ranged between 
0.7 and 1.3 with average approx 1.0. 

Long-term failures were observed the centrifuge maintaining accelera- 
tion significantly lower than that required cause short-term failure. five 
tests, failures occurred after continuous operation. Successful modeling 
long-term stability centrifuge implies that the delay failure caused 
process such hydraulic flow which occurs more rapidly the centrifuge 
model than the prototype. The critical state theory long-term stability 
thereby supported these test results. However, more work required 
accurately model long-term stability problems the drum centrifuge. 
future work will important reproduce the prototype stress path more 


accurately. This may require the ability cut slope while the centrifuge 
spinning prevent complete rebound the condition zero vertical 
stress. 
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LIMESTONE FOUNDATION 
DESIGN AND 


INTRODUCTION 


The purpose this paper consider the techniques and peculiarities 
involved designing and constructing shallow foundations bearing the Miami 
Limestone. The thickness and lateral extent the Miami Limestone are 
illustrated. The geologic characteristics this formation are described, are 
the types and properties Miami Limestone. Methods analysis for shallow 


foundation design are presented. Construction problems and difficulties are 
considered with case histories. 


The Miami Limestone soft, light tan, porous, sometimes sandy, fossiliferous, 

oolitic grainstone. The term Miami Limestone rather than Miami oolite will 
used because the formation encompasses sparsely oolitic fossiliferous 
limestone (1). The Miami Limestone exposed along the southeastern portion 
Florida, the Miami area, and also appears the lower Florida Keys. 
The lateral extent the Miami Limestone shown Fig. 

The Miami Limestone obtains maximum thickness approx (14 
along the Atlantic Coastal Ridge and thins abruptly near the Atlantic coastline. 
West the ridge, thins gradually and has average thickness about 
m-4.6 m). generalized geologic cross section constructed east-west 
through Dade County shown Fig. 

The Miami Limestone was formed approx 120,000 ago. that 
time, the sea level stood (7.6 higher than does today, and the Miami 
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area was tidal shoal bounded the east barrier sand bar. This environment 
facilitated the formation concentrically layered sand size carbonate grains 
called oolites. These grains formed repeated precipitation calcium carbonate 
about nucleus sand shell. Two geologic facies are present the Miami 
Limestone formation: the barrier bar oolitic facies and the tidal shoal limestone 
facies. 

The barrier bar oolitic facies consists both cross-bedded and chalky limestone. 
Photographs these two Miami Limestone types are shown Fig. The 
cross-bedded grainstone consists sets crossbeds separated intermittent, 
thin (25-mm) less] cemented layers mollusk fragment grainstone. 
The beds generally dip the southeast the seaward side and the northwest 
the inland side the ridge. The dip the crossbeds can steep 


~ 
- 


CHALKY 


CROSS-BEDDED 


~ 


FIG. Showing Cross-Bedded, and Porous Miami Limestone 
(Barrier Bar Oolitic Facies) 


Layers uncemented oolite have high intergranular porosity and are planes 
20°-30°. Zones high porosity are parallel bedding plains this facies. 
selective dissolution. This gives rise systems small tubular interconnected 
voids the rock. 

The second significant facies the Miami Limestone formation the tidal 
shoal limestone facies. The tidal shoal deposits generally underlie the oolitic 
barrier bar and outcrop some areas the formation. example tidal 
shoal limestone presented Fig. The tidal shoal limestone facies separated 
from the barrier bar distinct contact sometimes marked thin 


— 
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(25-mm) less], cemented layer shell fragment grainstone. The tidal shoal 
facies composed essentially the same material the barrier bar facies 
with several important differences. Oolites are smaller, less abundant, and 
relatively less developed. Fossils are more abundant and diverse. Bedding 
indistinct nonexistent. The absence primary sedimentary structures can 
attributed the action burrowing organisms which have churned the 
sediments. addition disturbing bedding, the burrowing has affected the 
porosity the facies. Burrowing organisms rework the unconsolidated settle- 
ments the environment deposition. The abandoned burrows were later 
filled with coarse calcareous material (10). Relict burrows act permeable 
channels the rocks the tidal shoal facies. Above the water table, these 
channels remain with uncemented shell fragments and oolites. Below the water 
table, ground water has often removed the filling material and enlarged the 


FIG. Diameter Core Very Porous Miami Limestone (Tidal Shoal 
Facies) with Length in. (610-mm) 


burrow channels, giving the rock the appearance Fig. 
example very porous Miami Limestone. 


Properties 


Rock Coring.—The writers have found that through the use large 4-in. 
(100-mm) diam rock coring equipment, relatively intact samples Miami 
Limestone can recovered, even where the porosity fairly high the 
cementation moderately low. The percent core recovery generally ranges 
between 75% and 100% upper two-thirds the Miami Limestone formation, 
and between 20% and 50% the lower portion. The core recoveries are generally 
function the rock porosity. the solution channels increase with depth, 
the core recoveries decrease. The Rock Quality Designation (RQD) values are 
generally above 50%-80% the upper part the formation and usually near 
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the lowest portion the formation. The RQD value dependent the 
solution channel size, porosity, the extent cross-bedding the Miami limestone, 
and the method coring used. 

Unit Weight.—Dry unit weight values for the Miami Limestone can vary 

Unconfined Compression Tests.—A histogram Miami limestone unconfined 
compressive strength presented Fig. 

Rock Modulus.—A plot stress versus strain was used determine the 
modulus elasticity. Rock modulus values are presented the form 
histogram Fig. The stress-strain curves are concave upward with the modulus 
values calculated from the straight line portion the plot. 

Splitting Tension Tests.—The difficulties associated with performing direct 
uniaxial tension tests rock have led number indirect test methods. 


NUMBER OF SAMPLES 663 
MEDIAN STRENGTH 31 KSF 


PERCENTAGE OF TEST RESULTS 


UNCONF INEO COMPRESSIVE STRENGTH. KSF 


FIG. 5.—Histogram Ultimate Unconfined Compressive Strength for Miami Lime- 
stone from 4-in. (100-mm) Diam Rock Core Test 


The method used the writers collect tensile strength data the Miami 
Limestone was diametral compression test. The rock core placed between 
the platens compression test machine and failed splitting tension along 
the vertical diameter. Tensile strengths measured this way are reproducible 
and are reasonable agreement with values obtained from uniaxial tension 
tests. Methods proposed Jaeger and Cook (2) are used derive the rock 
splitting tensile strength from the measured compression load. Splitting tension 
values for the Miami Limestone are presented the form histogram 
Fig. 

seen Fig. the Miami Oolitic Limestone often has considerable tensile 
strength. Consequently, when reasonable thickness this formation overlies 
less stiff material, such loose cohesionless sand, the limestone can act 
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analyzing the Miami Limestone for its ability support shallow 
foundations, the writers suggest that least four. criteria considered: (1) 
Punching diagonal tension failure the limestone mat; (2) local crushing 
the contact the foundation and the limestone; (3) beam tension 


NUMBER OF SAMPLES 688 
MEDIAN MOOULUS 12,000 KSF 


PERCENTAGE OF TEST RESULTS 


ROCK MODULUS VALUES 1000 


FIG. 6.—Histogram Rock Modulus Test Values for Miami Limestone from 4-in. 
(100-mm) Diam Rock Core Test 


NUMBER OF SAMPLES 1088 
MEDIAN STRENGTH 


PERCENTAGE OF TEST RESULTS 


SPLITTING TENSION VALUES, 


FIG. 7.—Histogram Splitting Tension Test Values for Miami Limestone from 4-in. 
(100-mm) Diam Rock Core Test 
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failure the limestone mat; and (4) settlement the limestone mat and the 
underlying soils. 

Punching Diagonal Tension Failure.—Fig. shows punching diagonal 
tension failure the Miami Limestone formation due excessive foundation 
load. This type foundation failure more likely occur the ratio 
H/B small. This can occur when the foundations are placed deep into the 
Miami Limestone, even though the total thickness the limestone large 
comparison the footing’s size. 

For computing factor safety against punching type failure, suggested 
that cylindrical failure surface assumed. The surface area the cylinder 
should equal the perimeter the footing times the limestone thickness below 
the foundation. The ultimate resistance punching shearing this cylinder 
through the Miami Limestone formation may computed multiplying the 
shear strength the rock the surface area the sheared cylinder. 

example punching type failure was given Sowers (11,12). this 
example, footings for building the South Campus the Miami-Dade 
Community College were placed deep into the Miami Limestone formation. 


Insufficient punching shear resistance the limestone produced diagonal tension 
cracks adjacent excavations. 


Oolitic 
Limestone 


Loose 
Sand 


FIG. 8.—Punching Shear Failure Miami Limestone 


Local Crushing.—This type failure occurs when the confined compressive 
strength the Miami Oolitic Limestone exceeded the contact pressure 
from the footing. The writers have observed this type failure small test 
footings when the contact pressure the order twice the unconfined 
compressive strength the rock. 

For footings with both axial and moment loading, Teng (13) provides method 
for theoretically calculating the peak contact pressure under the footing. The 
actual contact pressure, however, reduced significantly from the value computed 
theoretically for footings the Miami Limestone formation because: (1) The 
moduli elasticity the concrete and the rock are the same order 
magnitude and therefore there considerable bending the footings, which 
reduces the peak contact pressure; and (2) footings bearing the Miami Limestone 
are usually cast integrally with the rock around the perimeter the foundation, 
and this allows load transfer shear between the concrete footings and the 
limestone. 

possible that the actual peak contact pressure can the order 
only 1/4 1/3 that computed theoretical basis. Nevertheless, the 
writers suggest that the theoretical contact pressure limited the unconfined 
compressive strength the rock. Thus, minimum factor safety against 
local crushing type failure least would available. 
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example local crushing failure was documented Proctor (6) wherein 
footings the newly constructed Dade County Courthouse were believed 
have crushed the Miami Limestone formation. This crushing, however, was 
also associated with site dewatering and the removal portions weakly 
cemented material from the formation the footing bearing level. 

Beam Tension Failure.—Fig. shows failure the beam tension mode. This 
type failure occurs when the H/B ratio large and the flexural strength 
the limestone small. 

There are known documented failures this type the South Florida 
area. attempt was made induce this type failure during full-scale 
foundation load test the writers; however such failure could not 
accomplished. was found that such testing, the theoretically computed 
maximum tensile strength the bottom the limestone mat exceeded the 
laboratory splitting tensile strength test results 4-in. (100-mm) diam rock 
cores factor about This finding not completely unexpected, because 
common other cemented materials for the modulus rupture exceed 
the tensile strength. The theoretically computed maximum tensile stress the 
bottom the Miami Limestone can computed using slab stress analysis 
methods suggested Westergaard (15). These methods are summarized Roark 


(8). 


9.—Beam Tension Failure Miami Limestone 


Settlement Due Compression Rock and Underlying Soil.—For estimating 
settlement footings the settlement total structure bearing the 
Miami Limestone Formation, the writers suggest using elastic solutions for 
calculating deformations. Where there reasonable certainty that there will 
not beam tension punching shear type failures, mat (beam) elastic 
foundation analysis considered appropriate. Poulos and Davis (5) provide 
charts, tables, and formuli readily calculate these deformations. The writers 
suggest that the elastic moduli from the Miami Limestone obtained from 
laboratory compression tests the rock. Local experience indicates wide 
range variability for South Florida rock conditions. 

For the underlying cohesionless soils, the writers have found reasonable 
correlation between actual and calculated settlements using the correlations 
developed Webb (14). Webb’s correlations compare elastic moduli compres- 
sion for cohesionless soils with static cone resistance. 

computing settlements foundations bearing the Miami Limestone, 
the effects loading adjacent foundations must always considered. Because 
the beam action the limestone, considerable stress overlap occurs even 
for widely-spaced foundations. For heavy structures, such multistory buildings 
storage tanks, the writers have found that the beam action the Miami 
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Limestone has little effect the total amount settlement that occurs. Such 
analysis can done using Westergaard stress distribution; 
moduli elasticity developed from laboratory test rock cores, and static 
cone penetration tests cohesionless soils. 

example settlements under heavily loaded structure was described 
Kaderabek and Reynolds (3). this case, settlements earth structure 
with average loading approx ksf (190 were monitored and 
compared with settlements calculated using the conventional stress-strain method 
noted previously. 


noted previously, reduction the contact pressure between the bottom 
foundation and the Miami Limestone can achieved casting the 
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FIG. 10.—Mohr Diagram for Miami Limestone Laboratory Test Values 


foundations against the side rock excavations. this case, the contact pressure 
reduction due load transfer through shear between the sides the 
footing and the rock. The Miami Limestone typically has considerable shear 
strength. The writers have found that even where rock core recovery poor, 
considerable shear strength often still present. Where rock cores can 
recovered and tested the laboratory, estimate the shear strength 
the rock can obtained performing splitting tensile and unconfined compres- 
sion tests the rock cores, and constructing Mohr envelope representing 
the results series these tests. Fig. shows such test results. 

Where impossible recover testable rock core samples using 4-in. 
(100-mm) diam larger core barrel, the writers suggest that in-place shear 
strength tests are appropriate. Schmertmann (9) used test called 
Tension-Shear obtain the shear strength rock the Florida Keys. 
The writers have had considerable experience using similar test the Miami 
Limestone. Fig. shows schematically the test arrangement. Such tests have 
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been instrumented determine the stress distribution along the anchor and 
approximate stress-strain relationship during the test. Typical test results 


are shown Fig. 12. 
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FIG. 11.—Anchor Pullout Test Arrangement 
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Because the Miami Limestone has considerable shear strength this property 
may incorporated into the design foundations bearing the limestone. 
Examples foundations would shallow foundations with rock anchors located 
along the perimeter increase the overturning capacity the foundation, and 
anchors under shear walls building where wind pressures can result 
uplift forces. The writers also believe that most the stresses transferred 
the subsurface piles, such augercast piles and pressure injected footings 
(when installed through the Miami Limestone), are carried shear the 
limestone. Their opinion about load transfer these pile types substantiated 
very low settlements during pile load tests. 


Construction 


Construction the Miami Limestone quite unique. Advantage can taken 
the limestone’s properties and this can lead significant construction cost 
savings. the other hand, lack knowledge the limestone’s properties 
can lead severe cost overruns. 

the positive side, the material’s tensile and shear strengths often permit 
deep excavation walls cut into the Miami Limestone stand vertically. Fig. 
photograph vertical cut the Miami Limestone. 

less fortunate situation occurred the new Dade County Zoo, however, 
where vertical excavations were cut into and through the entire Miami Limestone 
formation, loose sand was encountered and excavated from below the limestone, 
and subsequently block failures the limestone walls occurred. Fig. 
photograph such block failure. 

Although the Miami Limestone rock, the material can excavated for 
footings with conventional hydraulic backhoes and can ripped and bulldozed 
larger excavations. Drilling and blasting the Miami Limestone are rarely 
required. 

some area, the upper surface the Miami Limestone has been severely 
eroded. This erosion has caused surface undulation and the rock 
surface. the eroded areas the top the rock surface can vary more than 
(1.5 over distances less than m). The amount surface 
variation quite dependent location and appears occur sporadically even 
within general area. The writers have found that backhoe test pit excavations 
are about the best method for determining the existence the eroded limestone 
prior construction. 

Even where such preconstruction studies have been made, greater than expected 
surface variations have been encountered. For 4-story structure located just 
north downtown Miami, footing excavations uncovered 15-ft (4.6-m) diam 
irregularly shaped hole which extended deep into the Miami Limestone. The 
solution this site was dewater the hole, thoroughly clean the sides and 
bottom the hole, and backfill with lean concrete. 

Deep and isolated rock slots have been observed depths (4.6 
m). Some these rock slots with diameter about in. (100 mm-200 
mm) have been observed contain peat and other organic materials. Usually 
these conditions occur closer Biscayne Bay. Engineers not familiar with 
subsurface conditions South Florida have mistakenly interpreted this condition 
site underlain with vast amounts organic material. When isolated deep 
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rock slots are encountered, they should cleaned depth 
diameters and backfilled with lean concrete. 

The ability dewater excavations into the Miami Limestone depends 
the type limestone present. Figs. and illustrated high and low porosity 
Miami Limestone. the high porosity Miami Limestone present site 
above the desired level dewatering, such dewatering will almost impossible 
conventional pumping methods. This was the case project western 
Dade County where undercutting was necessary. The contractor elected utilize 
drag line remove the unsuitable soils, which extended about (1.2 


FIG. 13.—Vertical Cut Miami Limestone 


FIG. 14.—Block Failure Miami Limestone 


beneath the ground-water table. order construct footings and place fill 
soils the dry, was decided dewater the 100-ft (30-m) square excavation. 
20-in. (508-mm) diam pump was used remove water from the excavation. 
With this pump operating maximum capacity, drawdown the water 
surface was possible. This fact necessitated extra expense assoicated with 
setting forms and placing concrete (1.2 water. 

Very porous Miami Limestone oiten encountered near the ground-water 
table, but this not always the case. project located downtown Miami 
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contractor needed establish temporary concrete batch plant site. The 
contractor felt sure adequate source water could obtained excavating 
shallow sump just below the water table. Unfortunately, the excavation 
was pumped dry with 4-in. (100-mm) pump matter minutes. 

The ability dewater the Miami Limestone can determined accurate 
rock identification and performing pump tests backhoe excavated test 
pits drilled holes. These tests yield estimations inflow rates. Caution should 
used dewatering cross-bedded Miami Limestone since the sand and silt-sized 
particles can removed, affecting the performance nearby foundations. 
The pump discharge water should checked frequently determine fines 
are being removed from the formation. 


The lateral extend and thickness the Miami Limestone along with its 
engineering properties have been presented. The writers suggest several criteria 
for evaluating shallow foundations bearing the Miami Limestone. Utilizing 
these criteria, evaluation shallow foundations possible. 

Since the Miami Limestone variable strength, porosity, and thickness, 
adequate investigation these properties should carried out prior 
foundation construction each particular site. 
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SHEAR STRENGTH ROCKFILL 


INTRODUCTION 


important question which arises during stability analysis rockfill dams 
the relative frictional resistance the rockfill and the underlying rock 
foundation interface. Attempts estimate the shear strength smooth, ice- 
polished interfaces has led examination the related shear behavior 
rock joints. found that rockfill, interfaces, and rock joints have several 
features common, including dilatant behavior under low effective normal 
stress, and significant crushing contact points with reduced dilation high 
stress. each case failure resisted strongly stress dependent friction 
angles. example, the peak drained friction angle rockfill near the 
base high dam may low 35°, while close the toe the same 
rockfill might exhibit value high 60°. Similar stress dependency 
observed with rock joints and interfaces. 

this article shown how the value for rockfill can estimated 
from knowledge the following parameters: (1) The uniaxial compressive strength 
the rock; (2) the particle size; (3) the degree particle roundedness; 
and (4) the porosity following compaction. The degree particle roundedness 
and the porosity determine the magnitude the structural component strength. 
This component increases shear resistance much the same way interlocking 
asperities rough joint surface. The structural component strength 
strongly stress dependent. added the basic angle friction flat 
nondilatant, i.e., sawn, surfaces the rock obtain 

The methods developed here are simple practical nature, and allow 
dam designer obtain preliminary estimate the peak drained friction angle 
rockfill, whether consists angular quarried rock, moraine, well-rounded 
fluvial gravels. Simple large scale tilt tests in-place rockfill are suggested 
for checking the shear strength different lifts dam, during construction. 
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SHear Rock Joints 


practical method for describing the peak shear strength three typical 
classes rock joints shown Fig. (2). The general equation follows: 


Cc, 

which peak shear strength; effective normal stress; JRC joint 
roughness coefficient; JCS joint wall compression strength; and residual 
angle friction. 

JRC.—The joint roughness coefficient dimensionless number, approxi- 
mately equivalent the ratio roughness amplitude and sample length. 
Completely planar smooth joints, sawn surfaces have JRC equal zero. 
The most reliable method measuring JRC conduct tilt tests, preferably 
with rock blocks natural size containing through-going joints. Sliding occurs 
due the effect gravity the upper block. The JRC back-calculated 


which tilt angle when sliding occurs (typically 50°-75°); 
effective normal stress acting across the joint when sliding occurs. 


joints, rough sheeting, rough | sheeting, planar shear joints, planar 
bedding. | foliation and bedding. | foliation, bedding. 
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FIG. 1.—Method Estimating Peak Shear Strength Rock Joints, Based Joint 
Roughness Coefficient JRC (20, 10, 5), and Joint Wall Compression Strength 
JCS (100 MPa, MPa, MPa, MPa) (2) MPa 145 psi) 
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Later this paper shown how equivalent roughness (R) can 
estimated for rockfill, based the porosity, and the degree particle 
roundedness and surface smoothness. Alternatively, can back calculated 
from large-scale tilt tests in-place 

JCS.—Each envelope Fig. numbered with the appropriate joint wall 
compression strength (JCS) value units megapascals MPa 145 psi). 
The joint wall compression strength can measured directly with Schmidt 
hammer, estimated from the uniaxial compression strength the rock. 
Weathered joints have JCS values less than due the weakening effect 
the thin skin weathered rock along the two walls the joint. the 
joints are saturated, both (unweathered) and JCS (weathered) are somewhat 
reduced. This reduction compression strength due water saturation 
important factor explaining the long-term settlement rockfill dams, 
explained Terzaghi (25) and the second writer (12). 

Later this paper shown how equivalent compression strength (S) 
can obtained for rockfill, based the uniaxial compression strength 
the rock and the particle size. Large particles prove weaker than 
smaller ones. similar scale effect for joint asperities has been demonstrated 
recent work reported Bandis, al. (1). 

and residual angle friction rock joint may reached 
after large shear displacements when joint asperities are worn down. shown 
30° Fig. for simplicity. practice, may range from about 15°-30°, 
depending the degree weathering and the mineralogy the rock. 
When joints are unweathered, JCS the minimum shear strength 
angle some degrees higher and called the basic angle friction 
usually ranges from about 25°-35°. can estimated indirectly conducting 
simple tilting tests, using dry, flat, nondilatant, i.e., sawn, surfaces the rock 
question. The value fundamental parameter intact rock, and 
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FIG. 2.—Triaxial Test Results for Rockfill (15) MPa 145 psi) 
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seems key parameter describing the shear strength 

The residual friction angle weathered rock joints very difficult 
determine experimentally, due the large displacements required. can 
nevertheless estimated with acceptable accuracy using Schmidt hammer 
and tilt test measure the parameters given the following empirical relationship: 


which Schmidt (L-type) hammer rebound saturated weathered joints 
interfaces; and Schmidt (L-type) hammer rebound dry, unweathered 
rock surfaces (for example fresh fractures sawn surfaces). 


25°. Experimental details are described greater detail the first writer 
and Choubey (4). 

Comparison with Rockfill Behavior.—Direct shear tests performed 130 
rock-joint specimens selected from seven Norwegian rock types gave the following 
mean values these key parameters: (1) JRC 8.9; (2) JCS 13,350 psi 
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(90 MPa), and (3) 27.5°. Eq. shows that the resulting value 
arctan under effective normal stress 145 psi MPa) would 
45°. 

review numerous triaxial tests rockfill materials Leps (15) shown 
Fig. indicates that ‘‘average displays value nearer 40° 
under the same stress level. The aforementioned rock joints seem display 
average peak shear strength comparable Leps’ packed, well-graded, 
strong rockfill. Note that the log-linear trend versus 
assumed Leps also predicted and measured for the case rock joints 
(Eq. 1). 

Mohr envelopes obtained from triaxial tests rockfill are nonlinear, and 
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FIG. Approach Shear Strength Estimation for Rock Joints, Rockfill, 
and Rock Interfaces 


very similar shape the peak strength envelopes obtained from shear box 
tests rock joints. This clear when comparing Fig. with Fig. (see Class 
joints). 

fundamental reason for this similarity behavior can found the 
observations contact area during shear tests rock joints. Shear tests that 
are interrupted the instant peak strength display apparent visible area 
(A,) damaged asperities that small fraction the total joint sample 
area usually assumed when converting force stress. The ratio 
failure occurring the points contact the instant peak strength, 
whether the stress very low, 0.001, very high, i.e., 


rock 
ROCKJOINT 
ROCKFILL 
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0.1. This similar principle the yielding contact points observed 
metal friction. 

extremely high stress JCS increases due the confining effect caused 
approaching A,. this stage the value JCS most accurately 
represented the triaxial strength the rock itself, and substitution 
this term Eq. gives very accurate fit high stress data. 

The maximum possible shear resistance rocks occurs when 
according detailed review data the first writer (3). This horizontal 
point the Mohr envelope reached when and 30°. 
implication, dilation totally suppressed, and behavior essentially 
ductile. 
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FIG. 5.—Measured Effect Particle Size (8) and Comparison with Estimated 
Values in. 25.4 mm) 


The apparent link between rock joint and rockfill behavior seen when 
triaxial test data are analyzed using the breakage factor’’ (B). This 
term was defined Marsal (18) the percentage weight particles that 
have undergone crushing following triaxial shear test. Data collected Marachi 
(16) indicate that approaches limiting value 3.0, approaches 
minimum 30° the particle breakage factor (B) extrapolated 100%. 
This limiting state cannot reached unless exceptionally high confining stresses 
are applied. practice the particle breakage factor seldom exceeds 40%. 


The empirical approach shear strength estimation followed the 
remainder this paper shown diagramatically Fig. will seeking 
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equivalent strength particles (S) and equivalent roughness (R) represent 
fully possible the complicated behavior rockfill. The following form 
Eq. therefore envisaged for rockfill: 


(=) 


Some important aspects the behavior rockfill will summarized 
provide frame reference for assessing the performance Eq. 

Crushing Particles: Scale al. (16,17) and Becker (8) 
observed that the following factors increased the amount particle crushing: 
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FIG. Effect Particle Size (17) and Comparison with Estimated 
Values in. 25.4 mm) 


(1) Increased water content; (2) increased uniformity; (3) increased angularity; 
(4) reduced particle strength; (5) increased effective confining pressure 
(6) increased shear stress under given confining pressure; (7) testing 
triaxial cell compared plane strain testing; and (8) increased particle sizes. 

Reduced shear strength with increased particle size particularly well 
documented result for rockfill, following the detailed studies the University 
California, Berkeley. Some the results these scale effect investigations 
are reproduced the upper half Figs. and The experimental data 
incorporates particle sizes ranging from 0.08 in. mm) using 
triaxial ceils with diameters ranging from 2.8 in.-36 in. (70 mm-910 mm). 

Scale effects appear marked this range particle sizes but apparently 
die out either end. For example, triaxial tests reported al. 
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(26) indicated noticeable scale effects for samples with less than 0.2 
in. mm), and crushing tests single sand grains with diameters from 0.016 
in. (0.4 mm-2 mm) reported Billam (9) suggest limited scale effect 
for sand below 0.028 in. (0.7 mm) size. 

upper bound suggested comprehensive review for uniaxial compres- 
sion test data for rock reported Lama and Gonano (14). For the most part, 
rock samples with volumes larger than (610 in.) displayed relatively 
insignificant scale effect, while smaller samples displayed considerable scale 
effect. The statistical volume theory Weibull (27) appears the explanation 
for these effects. The probability sample containing microfractures increases 
the volume increases. The fact that the larger rock samples may contain 
numerous microfractures and sand grains probably none, perhaps the reason 
for the S-shaped relation between crushing strength and particle size shown 
Fig. 

the light the preceding appears logical introduce size-dependent 
equivalent (S) for rockfill. This would consistent with scale effects 


TRIAXIAL TEST PLANE TEST 


size (mm) 


FIG. 7.—Method Estimating Equivalent Strength (S) Rockfill, Based Uniaxial 
Compression Strength and Particle Size 


observed from shear tests rock joints. Bandis al. (1) demonstrate how 
the joint wall compression strength (JCS) reduces with increasing size joint 
sample, due the mobilization larger asperities peak strength. Larger 
displacements are required reach peak strength, the length sample 
increased. 

Trial and error fitting numerous triaxial data with our prediction model 
(Eq. has led the empirical solution particle size effects shown Fig. 
should noted that the size will not necessarily the most relevant 
particle size for estimating the equivalent strength (S) given grading. The 
size however convenient description given rockfill. Shear strength 
differences between widely and narrowly graded rockfill with equal sizes 
will presumably arise chiefly from the different porosities attainable, and this 
will affect the structural component strength independent any possible 
changes equivalent particle strength (S). 

The example points and the triaxial test curve Fig. represent 
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sizes in.) and 240 (9.5 in.), respectively. These can 
considered representative laboratory samples and in-place rockfill, respec- 
tively. The following equivalent strengths (S) would appropriate the uniaxial 
compressive strength the parent rock was found 150 MPa (21,800 
psi): (1) Grading mm, 150 0.3 MPa (6,500 psi); 
and grading 240 mm, 150 0.2 MPa (4,350 psi). 

The great majority shear-test data for rockfill have been obtained from 
triaxial testing. However, plane-strain tests are more representative the strain 
conditions for many practical situations. common experience with such 
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FIG. 8.—Method Estimating Equivalent Roughness Based Porosity Rockfill, 
Origin Materials, and Degree Roundedness and Smoothness Particles 


tests that failure occurs relatively small strains, that less cruching particles 
occurs, and that values are least 2°-4° higher than the case for triaxial 
tests the same materials. The different scale effects for triaxial and plane-strain 
tests envisaged Fig. empirical but effective method accounting 
for these differences The proposed method will generally result about 
differences over the range stresses interest. 

Degree Particle Roundedness.—The cohesionless materials used rockfill 
dam will generally come from one more the following sources: (1) Quarried 
rock; (2) talus; (3) moraine; (4) glacifluvial deposits; and (5) fluvial deposits. 
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Each these potential sources produced characteristic degree particle 
roundedness and surface smoothness, which general dependent the 
length transport from the source origin. However, the mineralogy the 
rocks question sometimes complicates attempts correlate length transport 
with the degree roundedness. Softer rocks obviously attain rounded edges 
more rapidly than hard rocks both are transported equal distances. 

Investigations Pellegrino (21) nevertheless indicate that particles originating 
from talus are generally angular, those from glacial deposits are partly angular 
partly rounded, and those from fluvial deposits are rounded well rounded. 
Morris (20) was able induce degree roundedness artificially tumbling 
crushed particles basalt water over long period. The value measured 
special tilt apparatus thereby reduced from 53.8°-46.2°. 
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FIG. 9.—Effect Variable and Shear Strength Estimates Fig. 


Several classification systems for particle shape and roundedness have been 
suggested the past. Those Krumbein (13), Powers (23), and Pittijohn 
al. (22) were conveniently reviewed Chilingarian al. (10). The classification 
system presented Fig. combines several aspects these writers’ work, 
but specifically tailored rockfill rather than sand particles. The addition 
porosity and equivalent roughness (R) make possible estimate the 
peak drained friction angle given rockfill, using Eq. The method 
was developed through trial and error, backanalysis numerous triaxial 
(and some plane strain) tests rockfill materials reported the literature. 
Unpublished data from tests the Norwegian Geotechnical Institute were also 
included. 

Example Applications Equivalent Roughness Method.—The aforementioned 
method for estimating the peak strength rockfill requires the following input: 
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particle size 


equivalent strength (S) (Fig. 
uniaxial compressive strength 
rock 


description: origin 
roundedness 
smoothness 


equivalent roughness (R) (Fig. 
porosity after compaction 


basic friction angle (usual range 25°-35°) 


Figs. and show the potential the method. Estimated values, and 
their variation with confining pressure and particle size, are compared with 
measured data for Pyramid dam argillite (quarried rock) and Oroville dam 
amphibolite (fluvial gravel). Uniaxial compression strengths were 15,805 psi (110 
MPa) and 28,565 psi (200 MPa), respectively. 


TABLE Varies Degrees for Each Ten-Fold Change Pressure 


Equivalent 
megapascals overburden triaxial 


(4) 


Note: 25.4 in.; MPa 145 psi. 


shown Fig. can now study the effects variations porosity and 
rock strength since these variables have direct effect and 
respectively. Fig. some typical values and are plotted evaluating 
Eq. with 27.5°. The packed, well graded, strong 
Fig. are readily simulated high values and while low values 
both parameters provide the best fit the packed, poorly-graded, 
weak appears that Leps was correct drawing straight line envelopes 
inversely proportional log but may have been incorrect drawing 
parallel upper and lower envelopes. 

Shear Strength Rockfill under Extremely Low Stress.—The triaxial data 
reviewed Leps (15) included tests with less than about 0.04 MPa 
psi). Yet even lower values operate the toe and close the face 
rockfill dam. The associated high values are probably responsible for 
the extremely favorable resistance ravelling. Accurate determination 
under extremely low values clearly practical interest. 

Special plane-strain shear box tests angular particles crushed granite 


0.0001 67° 65° 
0.001 61° 59° 
0.01 500 53° 
0.1 49° 47° 
1.0 43° 41° 
10.0 500 37° 
ver 
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reported Seed and Goodman (24) indicated values high 67°, 
under effective normal stress approximately MPa (0.0145 psi). 

Extremely high values are also predicted the equivalent roughness 
method (Eq. these low stress levels. assume the following input 
data representative Seed and Goodman’s crushed granite: 7.5 
(0.3 in.); 150 MPa (21,750 psi); 35%; and 30°, the results 
presented Table are obtained. The equivalent roughness (R) six this 
example obtained from Fig. based assumed angular very rough 
crushed granite. The predicted values appear realistic 
over five orders magnitude stress, judging the high stress data shown 


25%) 


ASSUMPTIONS: 


*/.(R= 8. 
1SOMPa 1. ne (R= 8.0) 


2. n= 35% (R=5.5) 
200mm 


FIG. 10.—Estimated Variation under Toe and beneath Downstream Slope 
Dam 


The triangular shape rockfill dams results approximately triangular self- 
weight vertical stress distributions, and correspondingly high values close 
the toe, shown Fig. 10. Values excess 50° and 60° are 
estimated for the assumed porosities 35% and 25% these examples. 
also interesting see the high values close the downstream face 
the dam. this slope was analyzed infinitely long slope, with safety 
factor tan (in which slope angle), then clear that 
would quite high close the face. 

Variation Mobilized with Deformation.— While the main purpose this 
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paper the estimation peak for rockfill, clear that the strain required 
mobilize peak strength also considerable importance. 
analysis some forty triaxial tests has indicated relationship between 


(the axial strain peak strength), and ,). The latter can termed 
the structural component strength (i), which 


The test data analyzed indicated the following approximate ranges 

Thus general relationship for axial strain peak the following: 


The following example illustrates typical case: 


7.5 


quarried granite 


28% 


21,750 psi (150 MPa) 
7,250 psi (50 MPa) 
0.8 in. (20 mm) 


30° 
72.5 psi (0.5 MPa) 


peak 


TABLE 2.—Mobilization Equivalent Roughness (R) Rockfill with Increasing Axial 
Strain Axial Strain Any Instant) 


The mobilization shear strength rockfill strained towards peak shear 
strength can characterized the term Thus, the friction angle mobilized 
any given strain may expressed as: 


mov R nov’ 10g 


This concept was used Barton (6) generate shear stress displacement 
curves for rock joints for any desired loading unloading. the case 


(1) (2) 
0 —o,/i 
0.25 0.25 
0.50 0.75 
1.00 1.00 
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rock joints the term was used. Table shows how even extremely 
simple table values can used generate realistic stress-strain 
curves for 

example application this method shown Fig. 11. Measured data 
derived from vacuum triaxial tests the Norwegian Geotechnical Institute (Janbu, 
11) are compared with estimated data. Two different gradings are investigated, 
one representing the laboratory tests with d,, (0.24 in.), the other 
(stippled) representing the grading the dam where 250 (9.8 in.) 
Realistic values (saturated) 150 MPa (21,750 psi) and 25° were 
assumed for the schistose gneiss. Note the reduced strength and increased 
deformability predicted for the full-scale material. The important influence 


NGI vacuum triax. 

Svartevanndam | 

(schistose gneiss) 250mm 


MOBILIZED FRICTION (tan. 


0 5 10 Ss 5 
AXIAL STRAIN E (%) 


FIG. 11.—Comparison Measured and Estimated Mobilization with Increasing 
Axial Strain MPa 145 psi, in. 25.4 mm) 


compaction, porosity, strength and deformability clearly reflected 
both the measured and estimated results. 


Practicat 


The proposed method for estimating the peak strength rockfill requires 
the following input: 


particle size 


equivalent strength (Fig. 
uniaxial compressive strength 


rock 


description rockfill: origin 
roundedness 
smoothness 
equivalent roughness (Fig. 
porosity after compaction 


basic friction angle 


12 50° 
| Omod. 
/ %,=25 
/ Ww b 
of 
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The value can estimated from simple index tests such the Schmidt 
hammer point load test. The value can estimated from tilt tests 
sawn surfaces the rock question. usually ranges from 25°-35°, depending 
rock type. 

Preliminary Parametric equivalent roughness method recom- 
mended for obtaining preliminary estimates prior detailed strength 
investigations. particular advantage the ease with which parametric studies 
can performed. Shear strength envelopes can generated from Eq. 
investigate the effects varying the particle size, varying the compactive effort, 
varying the type rockfill. Eq. and Table can used generate 
approximate mobilized strain curves. 

Extrapolation Laboratory Shear Strength Data.—During the more detailed 
strength investigation that usually precedes construction rockfill dam, the 


PLACE THT SOX OW LEVEL 


ba 
BACK- CALCULATE (®) 


ISOLATE FROM SURROUNDING Fri MEASURE TILT ANGLE (¢) AT WHICH FAURE OCCURS 


FIG. 12.—Method for Checking Full-Scale Strength (and Porosity) Different Lifts 
Dam Trial Fill 


method can used extrapolate triaxial (or plane strain) strength data 
stress levels above and below those normally employed. This will particularly 
advantageous for the case large vacuum triaxial tests. general, 
recommended that the equivalent roughness (R) back calculated from the 
available test data, and the values and estimated from simple index 
tests. However, there may occasions when Fig. considered adequate 
for estimating (R), and therefore the value equivalent particle strength (S) 
would back calculated basis for extrapolating available test data. The 
assumed dependency the particle size, which obvious 
simplification, would then unnecessary. Equivalent particle strength (S) would 
back calculated rearrangement Eq. Thus 


i/R 


REMOVE SHEAR PLANE STIFFENERS AND BEGIN 
(8) 
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Extrapolation Large-Scale Tilt Tests In-Place Rockfill.—A serious limitation 
laboratory shear strength investigations the inability test full-scale rockfill 
samples with as-built gradings and porosities. Development this equivalent 
roughness method, with its apparent abilities for accurate extrapolation strength 
over many orders magnitude stress, has provided means interpreting 
large-scale tests performed under extremely low stress. 

practice, robust rectangular open box split horizontally could placed 
trial fill the partly completed dam, ready for sampling the next 
lift rockfill. would embedded, filled, and compacted the customary 
number passes vibrating roller. Subsequent excavation fill from around 
the outsides the box would release and allow the whole box, preferably 
several meters length, lifted hydraulically one end. Fig. shows 
the principle. 

The maximum angle tilt (a°) tolerated before failure, would probably 
the order 55°-65° under the extremely low effective normal stress operating 
across the failure surface. The angle would extrapolated design 
stresses estimating values and from index tests, and back calculating 
the value from Eq. Note that errors estimating and are automatically 
compensated the values back calculated. Tilt tests have the added 
advantage that the nonuniform strain and progressive failure common conven- 
tional shear box tests much reduced, due the inherently more uniform 
nature gravitationally induced shear and normal stress. 

important by-product the test would the relative ease with which 
in-place porosity could measured. This would calculated from the known 
volume the tilt sample, and from the force required lift one end, i.e., 
one half, the tilt box, before and after filling. 


Summary 


Shearing between the two interlocked walls rock joint, between 
the interlocked particles rockfill, involves significant degree overriding 
(dilation) under low effective stress and significant degree crushing 
contact points under high effective stress. The peak shear strength the two 
materials surprisingly similar, although rock joints generally reach their peak 
strength much smaller strains than rockfill. 

The peak drained friction angles rock joints and rockfill can quantified 
equivalent roughness (JRC R), equivalent asperity particle 
compression strength (JCS and the residual basic friction angle 
respectively. The latter are generated nondilatent surfaces 
the given rock types. Each these six parameters can estimated simple 
index tests, classification. the case rockfill, the porosity, the degree 
particle roundedness, the size, and the uniaxial compression strength 
the relevant rock are all essential the estimation the peak drained 
friction angle. 

The value for both rockfill and rock joints found inversely 
proportional the log effective normal stress. Thus for rockfill increases 
some 4°-8° for each ten-fold reduction effective normal stress. Extremely 
high values are observed for very low effective stress levels. This stress 
dependency very positive factor the critical toe region rockfill 
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dam. High values this region help explain the high resistance ravelling 
during extreme leakage. 


The satisfactory fit between estimated and measured values over several 


orders magnitude stress, suggest that large-scale tilt testing in-place 
compacted rockfill may promising economic method checking full-scale 
strength and porosity different lifts dam, trial fill prior dam 


construction. 
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The following symbols are used this paper: 


apparent (visible) area damaged asperities when rock joint sheared 
peak strength; 
total area rock joint sample assumed when converting force 
stress; 
particle breakage factor after Marsal (18); 
maximum particle size; 
particle size with 50% passing; 
safety factor; 
structural component shear strength for rockfill, where 
joint wall compression strength; 
roughness coefficient rock joints; 
porosity, percentage (after compaction); 
equivalent roughness rockfill; 
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Schmidt (L) hammer rebound number, measured saturated, weath- 
ered joints; 

Schmidt (L) hammer rebound number, measured dry, unweathered 
rock surfaces; 

equivalent strength rockfill particles; 

maximum angle inclination recorded when upper half jointed 
block rock (or rockfill sample) slides under influence gravity 
loading, during tilt test; 

axial strain, percentage; 

uniaxial compression strength measured cylindrical rock samples; 
effective normal stress; 

effective normal stress acting when sliding occurs during tilt test; 
effective major principle stress; 

effective confining pressure (applied triaxial test); 

peak shear strength; 

basic friction angle smooth, planar unweathered rock surfaces; 
residual friction angle weathered rock joints; and 

peak drained friction angle (for rock joints rockfill). 


Subscripts 
mob mobilized value (of JRC any given strain shear 
displacement; 
inst instantaneous value (of axial strain and 
peak value peak strength (of R). 
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CLAY ANISOTROPY AND BRACED WALL BEHAVIOR 


INTRODUCTION 


Over the past yr, the nature and significance anisotropy the undrained 
behavior soft medium clays has been studied detail. now recognized 
that strength and modulus values can vary much factor four 
for loading under vertical and horizontal directions. The typical undrained test 
results Fig. show the variation undrained shear strength for four soft 
Norwegian clays with the angle the major principal stress reorientation during 
loading The condition for represents the normal testing condition, 
with the major principal stress failure, acting vertically. 90°, the 
major principal stress, which begins vertical plane, rotated that 
failure horizontal. The results suggest that the undrained strength 
maximum 0°, and decreases steadily increases. The ratio the 
plastic clay and 0.2 for Drammen lean clay, with the Vaterland and Studenter- 
lunden clays yielding values between these extremes. Test results reported for 
clays located the continental United States show similar trends, but not 
dramatic strength changes (10,16). 

The behavior indicated Fig. and related test results has potentially 
significant implications for braced tied-back wall supporting excavation 
clay. terms classical concepts frictionless walls, the soil the 
passive side the wall undergoes 90° reorientation principal stress directions 
during excavation, while that the active side remains unchanged (2,15). 
this true general, the data Fig. suggest that the soil strength the 
passive side during excavation may low one-fourth that the 
active side. Bjerrum, Clausen, and Duncan (2) showed that for the few cases 
where earth pressures braced walls clay have been measured, the anisotropic 
behavior has considered order predict theoretical values consistent 
with those observed. 

view this information would appear important able evaluate 
the effects anisotropy basal heave potential, movements, and loads 

Civ. Engrg., Stanford Univ., Stanford, Calif. 94305. 
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VATERLAND CLAY 

CLAY 
PLASTIC CLAY 
ORAMMEN LEAN CLAY 


UNDRAINED STRENGTH RATIO, 


ANGLE MAJOR PRINCIPAL STRESS CHANGE 
THE VERTICAL, 


FIG. Undrained Shear Strength with Angle Stress Reorientation 
for Four Norwegian Clays (after Ref. 19) 


the wall-soil system. addition, several important questions arise: 


How have conventional designs based upon assumed isotropic behavior 
been successful the past? 
Under what conditions anisotropy important consider design? 


This paper attempts answer these questions, well provide techniques 
which can used for designer consider anisotropy predicting the potential 
for basal heave and movements braced wall systems soft medium clays. 
Limit equilibrium and finite element analysis techniques are utilized the primary 
investigative tools. 


excavation deepened clay, the shear stresses below the bottom 
the cut may increase the point that they exceed the shear strength 
the soil. Terzaghi (24) likened this shallow bearing capacity failure, and 
proposed the mechanisms and analysis procedures shown Fig. describe 
the basal heave condition. 1955, Bjerrum and Eide (3) showed that for the 
special case deep and narrow excavation underlain clay, different 
failure mechanism was applicable. They assumed the basal heave occur 
deep bearing failure which existed locally, near the bottom the cut. Neither 
the Terzaghi nor Bjerrum and Eide approaches considered anisotropy the 
clay behavior. 

For this study, the Terzaghi approach basal heave analysis extended 
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include anisotropy. particularly useful since accounts for the effects 
underlying rigid base consistent manner (see Fig. and found 
correlate strongly with excavation system movements. Later herein will 


(b) 


FIG. Heave Failure Mechanisms Proposed Terzaghi (24) 


used means predicting the movements supported excavations 
anisotropic clay. 
order consider anisotropy effects base stability analysis, necessary 
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allow for strength variation along the potential failure planes demonstrated 
Fig. Along the vertical plane extending through the overburden, the undrained 
strength for 45° stress reorientation Below the bottom the 
cut, shallow bearing failure occurs and the stress reorientations vary from 

For these conditions, and assuming undrained strength variation with 
which follows the patterns indicated the data Fig. Davis and Christian 
(9) developed bearing capacity factors account for anisotropy. Denoting the 
anisotropic bearing capacity factor distinguish from the normally 
used for isotropic soils, bearing capacity obtained multiplying s,, times 
Using all the foregoing information, the Terzaghi expression for basal 
heave factor safety, FS, may modified for anisotropic strength effects 


FIG. 3.—Rotation Principal Stresses Following the Terzaghi (24) Failure Mechanism 


the total unit weight the soil; and s,, and s,,, are the undrained shear 
strengths the clay and 45°, respectively. 

inconvenient for practical application because the strength 
45° not easily determined. more useful express s,,, terms 
and and this can done first defining ratio Next, 
(1/2) K,) (13). Using these definitions, Eq. may written 

FS A N? Suo 


(1) 
which the depth excavation; the width the excavation; 
Suo 
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which expresses the factor safety against basal heave for homogeneous soil 
conditions terms the undrained anisotropic strength variation clay. 
Parameter may found performing triaxial compression and extension 
tests and taking the ratio the strengths, may selected from data 
the literature (14,16). Term function K,, and values for are 
given Davis and Christian (9). Eq. may also modified for the effects 
underlying rigid layer which limits the depth failure approximate 
fashion simply substituting (see Fig. for 

order evaluate the impact anisotropy base stability, useful 
consider the ratio the factor safety for anisotropic soils and that for 
isotropic soils. Assuming underlying rigid base influence the failure 
mechanism, the ratio, given 


Su0 
Syo(1 K,) 


defining the factor safety for isotropic conditions given terms 
the undrained strength from triaxial compression test, the 
conventional practice. Eq. 4(a) shows that function two terms; the 
first which simple ratio the anisotropic isotropic bearing capacity 
factors. The second term more complex, but approaches one for soft 


oO 
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FIG. 4.—Ratio Factor Safety Against Basal Heave for Anisotropic Soil That 
for Soil for Wide Excavation 


inwhich and (4b) 
K,) 
100 
075 
050 
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medium clays the excavation width greater than (50 ft) and 
0.25 greater. Under such conditions, Eq. becomes Since 

Importantly, the data Fig. show that long less than one 
the factor safety against basal heave failure for isotropic behavior 
less than that for anisotropic behavior. The effect becomes larger 
decreases, with the factor safety for anisotropic conditions about 30% less 
than that for isotropic conditions for 0.25. For cases where the and 
parameters are outside the ranges where the term approximately one, 
the factor safety ratio below that indicated Fig. further accentuating 
the influence anisotropy. These results suggest that the basal heave factor 
safety defined conventional practice overestimates the factor safety 
for anisotropic soil, effect that becomes more important the degree 
anisotropy increases. The implications this information for general practice 
are considered subsequently this paper. 


Finrre 


addition being able evaluate basal heave potential, necessary 
for design predict excavation system movements and loadings. order 
quantify these parameters, the finite element method was used. beyond 
the scope this paper describe the details the procedures; the interested 
reader referred Hansen (13), and Hansen and Clough (14). Basically, the 
analyses attempt model the braced wall problem with minimum idealiza- 
tions. This involves simulating the soil-structure system and the actual step 
sequence excavation followed strut installation. Initial stresses are assumed 
characterized at-rest conditions with for all elements. each 
stage construction, the new stress level and degree stress reorientation 
for each element the finite element mesh are calculated and used determine 
new modulus and shear strength. 

The nonlinear elastic soil model used the program allows for changes 
response the shear stresses increase well anisotropy effects shear 
strength and modulus where stress reorientation occurs other than 0). Undrained 
shear strength assumed vary with the angle principal stress reorientation 
following the equation proposed Bishop (1): 


which and parameters used fit the equation experimental 
data. For soft medium clay, Eq. fits experimental data well using the 
any value are fitted hyperbola, which allows the tangent modulus 
defined mathematically. For the extreme cases and 90°, 
the tangent modulus equations are 

-2 

max 


2(S.0 qo) 
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—2 
max 


which and tangent moduli for and 90°, respectively; 
the maximum shear strain any stage the analysis; one-half 
the initial principal stress difference; and and the initial tangent 
modulus values for and 90°, respectively. Parameter used fit 
the hyperbola the actual stress-strain curve; typically varies from 


R, 


° 


ROTATION ANGLE, 


FIG. Variations Undrained Modulus Values with Angle Stress 
Reorientation 


for undrained response (11). Eqs. and were developed extending previous 
work Duncan and Chang (11), and Simon (23). 

values other than and 90°, the calculated using expression 
proposed Dunlop, Duncan, and Seed (12), i.e. 


which undrained modulus for stress reorientation undrained 
fitting parameter. The variation the modulus with shown Fig. 
Depending upon the ratio initial tangent modulus for that for 


GT? 899 
= ¢ 
=-0.25 
‘ 


900 JULY 1981 GT7 


90°, parameter termed the modulus may increase decrease with 
increasing Most test data suggest that the modulus should increase with 
B(K, 1.0), trend just opposite that found for undrained strength (see 
Fig. 1). 

Using 5-8, incremental finite element analysis, allows the strength 
and modulus values for each element the finite element mesh adjusted 
each stage construction account for nonlinearity and stress 
reorientation. The results the program were checked against theoretical and 
field evidence insure that reasonable behavior was obtained (13,14). 

should noted that previously Palmer and Kenney (20) performed finite 


SOFT 
MEDIUM, 
ANISOTROPIC 
CLAY 


BOTTOM OF EXCAVATION 


FULLY 
PENETRATING 
WALL 


RELATIVELY STIFF SOIL 
GREAT DEPTH 


FIG. 6.—Parametric Study Problem 


element analyses braced walls anisotropic clay. However, their work assumed 
only two types response, one for the soil the active side the wall 
and one for the soil the passive side the wall. the present study the 
soil free choose wide range responses depending upon the degree 
stress reorientation dictated the analysis. 


Parametric 


The basic excavation problem considered shown Fig. its width 
(40 ft) and the maximum depth (50 ft). The soil assumed 
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layer encountered. braced wall with four levels struts supports the 
excavation. The wall assumed penetrate the full depth the clay and 
pinned the rigid base level. 

Values utilized for the key soil and support system parameters are given 
Table The basic value listed Table that used for most the 


TABLE 1.—Basic Conditions and Variations for Parametric Study Braced Excavation 
Anisotropic Cley 


Parameter Basic values Variations 
(1) (2) (3) (4) 


(a) Soil Conditions 


Total unit weight, kilonewtons 

per cubic meter 20.0 None 
At-rest earth pressure coefficient 0.6 None 
Poisson’s ratio 0.49 None 
Undrained shear strength 0°), 

Anisotropic strength ratio 0.5, 0.75 
Initial tangent modulus 0°), 

kilonewtons per square meter None 
Anisotropic modulus ratio 0.5, 2.0, 4.0 
Anisotropic modulus variation 

parameter 0.0 
Hyperbolic correlation parameter None 


Soil-to-wall adhesion, kilonew- 
tons per square meter 

Shear stiffness, kilonewtons 
per square meter per meter 


Wall stiffness, kilonewtons 

per square meter per meter 7.93 
Strut stiffness, kilonewtons 

per square meter per meter 1.96 
Number strut levels 
Vertical strut spacing, meters 3.5 


Excavation width, meters 
Depth firm layer, meters 
Final excavation depth, meters 


analyses; variations considered parameter are given the column next 
the basic value. Initial conditions for all the analyses consist level ground 
surface, with the wall place, and at-rest stresses the soil. 

The undrained shear strength increases with effective vertical stress, 
manner consistent with the known behavior soft medium clay. Its 


(b) Interface Conditions 
0.0 2.0 
(c) Support Conditions 
None 
None 
None 
None 
(d) Excavation Geometry 
12.0 None 
30.0 None 
15.0 None 
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initial tangent modulus 1,200 times the undrained shear strength, value 
shown number studies appropriate for soft medium clay 
The basic wall stiffness that PZ-32 sheetpile, and the basic strut stiffness 
falls between values backcalculated for field cases with similar support systems 
(13). This latter value allows for the effects connectors and shims which 
reduce the effective strut stiffness from its theoretical value. 

the finite element mesh, interface element used between the wall 
and the soil allow for relative movements. The interface assumed 
frictionless have strength following results field tests reported 
the Norwegian Geotechnical Institute between sheetpiles and soft clays (17). 
The interface mobilized its full resistance relative displacement 0.25 
(0.01 in.). 


System 


The effects anisotropy braced wall behavior can illustrated 
comparing the results two the finite element analyses, where the only 
difference that one case uses isotropic soil while the other uses soil 
which has anisotropic strength variation. The problem considered conforms 
the basic values listed Table the anisotropic soil utilizes 0.5 
0.5). The analyses described are typical excavations soil having 
shear strength representative medium clay. The value employed 
for the anisotropic soil intermediate the limiting values reported (13). 

Movements.—Lateral wall movements and ground surface settlements 
excavation depths and (30 and ft) for the two analyses 
are shown Fig. General forms the curves are similar for both cases, 
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with maximum lateral wall movements occurring slightly below excavation 
bottom. There settlement trough that forms the ground surface short 
distance behind the wall and settlements thereafter decrease with distance from 
the wall. The settlement the ground surface the wall relative movement 


40-80 25-40 


75-85 STAGE FS=1.3 
=0.5 


FIG. 8.—Angles Stress Reorientation, for Anisotropic Analysis 


along the frictionless interface because the wall undergoes settlement 
all since its tip founded top the underlying rigid base. 
Despite the similarity form the isotropic and anisotropic movements, 
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the absolute movements associated with the anisotropic soil (K, 0.5) are 
always larger than those the isotropic soil (K, 1.0). Importantly, the 
differences become larger the excavation deepens. The maximum lateral 
wall movement for the anisotropic soil reaches 200 (7.9 in.) full depth, 
while that the isotropic soil only (2.8 in.). 

The reason for the different levels movement related the effect 
stress reorientation the strength the anisotropic soil. the next section 
this paper will shown that 90° reorientation stresses occurs 
the passive zone front the wall; the anisotropic soil this leads the 
use one half the strength the isotropic soil that area. This reflected 
the relative factors safety against basal heave; depths and 
(30 and ft), 1.7 and 1.3 for the anisotropic case, but 
2.3 and 1.7 for the isotropic case. Since movements wall systems accelerate 
rapidly after the factor safety drops below 1.5 (8), the anisotropic soil shows 
much larger movements than the isotropic soil for the depth (50 ft). 

Stress Reorientation Patterns.—One the results the analyses primary 
interest relative the anisotropic model the degree principal stress 
reorientation, since this has strong influence the soil behavior. Fig. 
shows calculated degrees stress reorientation, for the basic excavation 
problem using the anisotropic soil with value 0.5 for the clay. Two 
excavation stages are depicted, one depth (30 ft) and the other depth 
(50 ft). Corresponding basal heave factors safety are 1.7 and 1.3, 
respectively. The results show rather complex pattern which only moderately 
affected going from depth (30 ft) and (50 ft). the passive 
side, large, typically between 85° and 90°. This consistent with the 
response expected based simple frictionless retaining wall theory. the 
active side however, calculated values range from 90°, contrast 
classical retaining wall theory which suggests that should 0°. close 
look Fig. reveals that the values which are significantly different than 
are largely confined area near the underlying rigid base, and area 
just behind the wall from m-9 (30 ft) depth. The stress reorientation 
both these cases caused boundary shear effects which are not considered 
classical retaining wall theory. The lower zone produced the movement 
the clay laterally with respect the unmoving underlying rigid base. This 
produces shear and stress reorientation. The upper zone develops due the 
sharp downward movements the clay the area the settlement trough 
(see Fig. 7). Here, the clay moves downward relative the soil either 
side it. Below this zone, depths and greater, the soil moving 
more laterally and less vertically, thus, stress reorientation concentrated 
the upper areas the soil mass behind the wall where the relative movements 
create boundary shear effect. 

Ojher than the two boundary shear regions, vaiues are small (0°-25°), 
suggested classical theory. The patterns stress reorientation observed 
Fig. are found essentially independent changes degree anisotropy 
the soil. There some accentuation the effect the interface between 
the wall and the soil assumed have friction, since the shear transfer 
increased. all cases, the average degrees stress reorientation are reasonably 


consistent with those assumed for the limit equilibrium basal heave analysis 
(see Fig. 3). 
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Yielding.—The zones soil undergoing yielding each excavation 
depth for analyses and are shown Fig. Yielding defined 
condition where the shear stresses element are equal within 
the undrained shear strength. expected, the yielding zones for the anisotropic 
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FIG. 10.—Comparison Earth Pressures for Anisotropic and Analyses 


soil (analysis are larger every stage than for the isotropic soil (analysis 
18). Yielding the passive side the wall begins for the anisotropic soil 
when the excavation only deep, and the entire passive side soil mass 
has failed when the excavation reaches depth 
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Earth Pressures.—Lateral earth pressures acting the wall for the isotropic 
and anisotropic analysis the final depth excavation are shown Fig. 
10. The anisotropic soil yields considerably lower passive pressure front 
the wall than the isotropic soil because its strength this area half that 
the isotropic soil. the active side the wall the anisotropic soil shows 
lower pressure than the isotropic soil below the bottom strut and higher 
pressure the level the upper struts. This reflects the larger lower wall 
movements for the anisotropic soil which lead greater redistribution 
arching effect the manner described Bjerrum, Clausen, and Duncan (2). 


Movement 


order examine the effect strength and stress-strain anisotropy 
excavation behavior the results series the parametric studies will 
examined. all cases, the wall taken PZ-32 sheetpile section which 
penetrates the full depth the clay. 

Fig. depicts the effect the anisotropic strength parameter, K,, 
movements. The maximum lateral wall movement the four excavation stages 
(depth 5.5 12.5 and predicted for the anisotropic soil 
normalized dividing the corresponding value for the isotropic case, 
and this ratio plotted against The results show that the maximum movement 
for the anisotropic soil always larger than the isotropic soil, with the difference 
increasing with lower values (greater anisotropy) and with increasing 
depth excavation (lower basal heave factors safety). the full depth 
the cut, the anisotropic soil yields maximum lateral wall movements that 
are 2.8 and 1.4 times those the isotropic soil for values equal 0.5 and 
0.75, respectively. Similar trends are also observed for surface settlements. 
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FIG. 11.—Effect Anisotropic Strength Ratio, Maximum Lateral Wall Move- 
ment 
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The effect the anisotropic stress-strain parameter, much more subtle 
than that the strength parameter, Fig. shows the ratio the maximum 
lateral wall movement for seven analyses with 1.0 that for 
1.0 plotted against the value used the analysis for excavation depths 
5.5 and 12.5 (18 and ft). The parameter varied from 0.5-4.0 
while varied from 0.5-1.0. The results show that for any value 
the value has only relatively small effect the movements. Slightly 
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FIG. 12.—Effect Anisotropic Modulus Ratio, Maximum Lateral Wall Move- 
ment 


larger movements result less than one, and slightly smaller movements 
occur greater than one. 

The interrelationship between lateral wall movement, excavation depth, aniso- 
tropy and factor safety against basal heave defined Fig. 13. Maximum 
lateral wall displacement normalized dividing excavation depth, and 
plotted against factor safety against basal heave for anisotropic and isotropic 
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analyses. The factor safety defined terms the Terzaghi (24) approach 
modified earlier this paper for effects anisotropy. Significantly, the 
results collapse single curve. This demonstrates that the factor safety 
against basal heave, properly modified for anisotropy, controlling variable 
determining the amount movement supported excavation clay will 
undergo. Below factor safety 1.5, the nondimensional movements increase 
rapidly, reflecting the strong effect the subsoil yielding and nonlinear soil 
behavior. Clough, Hansen, and Mana (8) have demonstrated similar trends 
field data. The data Fig. also suggest that there inherent difference 
between wall behavior isotropic and anisotropic clay once the difference 
factor safety considered given excavation depth. 

should noted that maximum ground surface settlements follow the same 
patterns shown Fig. for lateral wall movements. The settlements are 
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FIG. 13.—Relationship Between Maximum Lateral Movement and Basal Heave Factor 
Safety for and Anisotropic Analyses 


about half the lateral wall movements when 2.0 and above and about 
equal them below 1.5. 


common practice design excavation support walls using apparent pressure 
diagrams. These are diagrams pressure from which design strut loads are 
derived; commonly used examples are given Peck (22), and Peck, Hanson, 
and Thornburn (22). Fig. 14, the apparent pressures calculated six 
the fimie element analyses are compared values recommended Peck, 
Hanson, and Thornburn (22). The finite element results were chosen demon- 
strate the effect shear strength with values 57.9 
28.4 and s,, 8.8 For each shear strength, isotropic 
and anisotropic case considered, with the parameters, 0.5 and 
2.0, for the anisotropic soil. The lowest shear strength case was analyzed 
specifically demonstrate extreme condition, since the factor safety 
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against basal heave reaches 1.0 for the anisotropic soil before full depth 
achieved and about 1.0 for the isotropic soil full depth. 

predicted the Peck, al. (22) theory, the apparent earth loads from 
the finite element results increase shear strength decreases. And, almost 
every instance the apparent pressures for the isotropic analyses lie within the 
bounds defined the Peck, al. (22) pressure envelopes. The same also 
holds true for the anisotropic analyses long the factor safety against 
basal heave above about 1.4. However, the factor safety drops below 
this, the cases the lower excavation stages with s,, 28.4 
and 8.8 the apparent pressures clearly exceed the Peck, 
al. (22) limits. these instances, the soil the passive zone failing completely 
(see Fig. causing the struts have carry more load. For the lowest 
shear strength case with anisotropic soil the large inward soil movements below 
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FIG. 14.—Apparent Earth Pressure Diagrams for Selected Anisotropic and 
Analyses 


the excavation causes enough rotation the upper part the wall away from 
the upper strut generate negative loading it. 

Field evidence for strut load response similar that observed Fig. 
rare, since most cases appear bounded the design apparent pressure 
diagrams. However, unusually large loads lower struts were reported for 
several excavations Norway where low basal heave factor safety was 
reached (17,18,21). The reasons for this behavior have never been fully explained, 
but given the apparent strong anisotropy the Norwegian clays (see Fig. 
and the results predicted herein, would not unexpected. 


From the information presented effects anisotropy braced wall 
behavior, consistent trend emerges. Namely, that anisotropy clay becomes 
increasingly important the apparent factor safety against basal heave, 
based s,, only, falls below about 1.5, and the degree strength anisotropy 
increases. The results show that for strongly anisotropic clay (K, 0.5), 
the basal heave factor safety may reduced 20%-50% from its apparent 
value, depending upon the size the excavation and the shear strength 
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the soil. Also, movements can doubled, earth pressure distributions altered, 
and strut loads significantly increased relative those for presumed isotropic 
condition. Given these circumstances the question arises, have conven- 
designs managed successful the past not considering 
There are two responses the question. First, conventional designs have, 
fact, not been free from problems situations where stability conditions are 
marginal. There are more than few published and unpublished case histories 
failures such instances (4,6,20,21). Second, unusual for condition 
arise where: (1) The soil highly anisotropic; and (2) the basal heave factor 
safety low. The evidence available today suggests that few clays are 
anisotropic the data from the Norwegian clays Fig. imply (14,16). 
Furthermore, because inherent limitations our sampling and testing, the 
strength soil often underestimated, leading predictions low factor 
safety which often does not exist. Thus, the conditions where anisotropy 
can play significant role are not common. However, given the potential for 
anisotropy lead behavior that less desirable than expected, the engineer 
should especially cautious marginal stability situations. Existing data 


anisotropy soils, and the information presented herein can used guide 
design. 


Summary 


recognized that most natural clays are some degree anisotropic, and 


there growing trend incorporate this aspect behavior into geotechnical 
analysis techniques. the case braced tied-back wall the behavior 
the clay the passive side the wall can decidedly different than 
that the active side the wall the clay strongly anisotropic. the 
passive zone, failure the soil will occur with the major principal stress directed 
horizontally, 90° reorientation from the condition normally applied conven- 
tional testing. The strength soft medium clay loaded this manner 
known less, some cases four times less, than that for clay 
loaded simply increasing the vertical stress failure. This important 
since the existence lower strength zone the soil not considered 
conventional basal heave analysis procedures which assume isotropic soil 
behavior. Thus, the conventional techniques theoretically lead overestimate 
basal stability. 

this paper, technique presented which modifies basal heave analysis 
procedures account for anisotropy. Also, finite element analyses are used 
demonstrate effects anisotropy general system behavior. The results 
show that the soil presumed isotropic when truly anisotropic 
the following trends develop: 


The basal heave factor safety may much 50% above the actual 
one special cases where the clays are strongly anisotropic, and 10%-30% 
above the actual one moderately anisotropic clays. 

Lateral wall movements and ground surface settlements can larger than 
expected. 

The distribution strut loads can changed with unusually high loads 
developing the lower struts cases marginal stability. 
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The potentially detrimental effects which can occur result neglecting 


anisotropy appear most strongly when the factor safety against basal heave 
below about 1.4. 


also shown that the primary influence the anisotropy increase 
the potential for basal heave relative isotropic conditions. When accounting 
made for the actual factor safety against basal heave, isotropic and anisotropic 
movement patterns can reconciled. There apparently inherent difference 
wall behavior result soil anisotropy. 

Past records performance for supported excavations soft soil conditions 
indicate many cases where unexpectedly large movements and strut loadings 
have occurred. These problems can be, part, explained terms the neglect 
the effects anisotropy. would appear important that consideration 


design given anisotropy clay the future, particularly cases 
marginal stability conditions. 
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The following symbols are used this paper: 


™ 


curve fitting parameter for modulus anisotropy; 

curve fitting parameter for undrained shear; 

excavation width; 

curve fitting parameter for undrained shear strength anisotropy; 
ratio bearing capacity factors; 

ratio terms anisotropic basal heave equation; 

depth from bottom excavation underlying rigid base; 
initial tangent modulus; 

tangent modulus; 

modulus stress reorientation; 

modulus 90° stress reorientation; 


(3, 
max 
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modulus stress reorientation; 

factor safety; 

depth excavation; 

ratio initial tangent modulus value for that for 
90°; 

ratio undrained shear strength for that for 90°; 
bearing capacity factor for isotropic soil; 

bearing capacity factor for anisotropic soil; 

ratio factor safety for anisotropic soil that for isotropic 
soil; 

curve fitting parameter; 

undrained shear strength stress reorientation; 

undrained shear strength 45° stress reorientation; 

undrained shear strength 90° stress reorientation; 

undrained shear strength stress reorientation; 

one half the principal stress difference for at-rest soil stresses; 
angle stress reorientation; 

maximum lateral deflection; 

maximum shear strain; 

total unit weight; 

effective vertical stress; 

initial major principal stress; 

major principal stress failure; 

normal stress; and 

shear stress. 
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FOUNDATIONS STRONG SAND 
OVERLYING WEAK SAND 


INTRODUCTION 


often the case that foundation rests soil consisting thin, 
strong sand layer overlying weak sand deposit. The strong sand layer may 
naturally existing artificially provided increase the bearing capacity 
and decrease the settlement the foundation. 

available solution such problem the literature (13,14) considers 
the bearing capacity the lower weaker layer, this layer surface 
stratum; assumed that the upper layer serves principally spread the 
load, thereby reducing its intensity the lower layer. The predicted ultimate 
bearing capacity obtained this method depends mainly the assumed load 
distribution through the upper layer; however, the results seem the 
conservative side since the shearing resistance the upper layer ignored. 
fact, some rupture surface must develop the upper layer deform 
with the lower layer failure. 

The writer recently presented rational approach solve cases footings 
layered soil (4,5,6,10). This paper extends the writer’s previous theory 
cover the case footings resting subsoil consisting strong sand 
layer overlying weak sand deposit. Design charts are presented for this purpose. 
The charts are simple use; knowing the shear strength the sand layers 
terms the plane-strain angle shearing resistance, the coefficient punching 
shear can determined. This coefficient allows the ultimate bearing capacity 
strip circular footings predicted will described later. 


The assumption involved predicting the theoretical ultimate bearing capacity 
that, the ultimate load, soil mass the upper sand layer, roughly 
truncated pyramid shape, pushed into the lower sand layer (9). The forces 
the assumed vertical punching failure surface the upper layer (of thickness 
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below the footing base) can taken the total passive earth pressure 
P,, inclined average angle acting upwards [Fig. Thus, for 
strip footing width and depth the ultimate bearing capacity approximately 
given 


which and the ultimate bearing capacities the strip footing 
very thick bed the lower and upper sand layers, respectively (10). Further 


which and the bearing capacity factors corresponding 
the angles shearing resistance, and the upper and 


Actual plone of failure 
H Assumed plone of failure 


strong sand 
weak sand 


(a) Theory (b) Test (scale 


FIG. Failure: Strip Footing Strong Sand Layer Overlying Weak Sand 
Deposit 


lower sand layers, respectively (8); and the unit weights the upper 
and lower sand layer, respectively. Further 


which the coefficient passive earth pressure. 
determine the theoretical value the angle shearing resistance, 
mobilized the assumed failure planes and the corresponding coefficient 


passive pressure, Kp, for any two layers sand with the upper layer being 
the stronger, the following arguments are considered (4,6): 


the analysis made the actual curved planes failure [Fig. 1(5)], 


B 
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the angle will equal If, however, the analysis made the assumed 
vertical planes, the angle, mobilized must less than failure has 
not taken place the assumed planes (11). 

Based the fact that the failure strain the upper strong sand layer 
less than that the lower weak sand layer, simultaneous occurrence 
the shearing failure both layers cannot take place and more strain after failure 
required the upper layer reach the lower layer failure strain value. 
Thus, the mobilized angle shearing resistance the upper sand layer may 
less than the peak value and may approach the residual value. This will 
reflected the value and thus Kp. 

The mobilized passive earth pressure the assumed vertical failure planes 
the upper sand layer decreases with decrease the lower weak sand 


assumed 
plane failure 


upper sand layer 


lower sand layer 


FIG. 2.—Local Angle Shearing Resistance Assumed Plane Failure 


layer strength. This can explained the fact that with decreasing lower 
sand layer strength, the vertical displacement the upper sand punching column 
increases and the lateral movements decrease, resulting decrease the 
passive pressure. This lateral movement may not sufficient for the maximum 
mobilization the passive pressure that would generated the full value 
the angle shearing resistance 


mathematical verification these arguments rather difficult. Also, 
equally difficult separate the effect explained herein evaluating the 
average mobilized angle shearing resistance, and consequently the mobilized 
passive pressure, Kp, the assumed failure planes. However, these difficulties 
may circumvented expressing the angle the dimensionless form 


X=H 
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the value the coefficient passive pressure varies with the angle 
(2), can stated that the total passive earth pressure, P,, (Eq. depends 
large extent the relative strength the upper and lower sand layers. 


Based the fact that the locally mobilized angle shearing resistance, 
the assumed failure plane point [Fig. must equal the 
angle shearing resistance the upper sand layer, since located 
both the assumed and actual failure planes. Moving down from the point 
along the assumed failure plane, the angle decreases with increasing 
depth the distance the actual failure plane increases (11) minimum 
value Ad, the point the soil interface. The value depends 
not only the lower layer strength, but also the relative strength the 
two layers. The variation the locally mobilized angle shearing resistance, 
with the depth, can linear curved (Fig. 2). 

order predict the distribution the locally mobilized angle shearing 
resistance the assumed failure planes, together with the value the 
coefficient, computer program was developed perform such analysis. 
this study different distributions were assumed (linear, circular, parabolic, 

ete.) and with each, the value the coefficient was varied from zero 
total passive earth pressure, was calculated from the following equation 


The value that was determined from Eq. was later equated the 
value given Eq. order determine the average mobilized angle 
shearing resistance, The solution was conducted trial and error, assuming 
value and then determining the corresponding value from the 
tables Caquot and Kerisel (2) Eq. until both sides agreed. 

comparison between the deduced angle shearing resistance, from 
experimental data with the theoretical values will enable the determination 


the existing distribution together with the corresponding value the coefficient 


order verify the theory presented, tests were conducted model strip 
and circular footings resting dense sand layer overlying loose well 
compact sand. 

The material used this investigation was well-graded medium course 
angular silica sand. The grain size distribution had uniformity coefficient 
2.76. 

The angles shearing resistance the dense, compact, and loose sand 
were determined from plane-strain tests, and triaxial tests. The results showed 
that the angle shearing resistance from the plane-strain test was 10.4% higher 
than the one from the triaxial test for the dense sand while were almost equal 
the case loose sand. These values are given Table 

order assure reproducibility the sand density throughout the testing 
program, procedure was developed raining the sand from certain height 
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give uniform and desired density. Dense packing was achieved raining 
the sand from height in. (910 mm) for each 3-in. (76-mm) layer 
means metallic sieve, in. (460 mm) diameter; compact packing was 
achieved raining the sand from height in. (150 mm) for each in. 
(25 mm) layer means funnel with rubber tube and end sieve. Loose 
packing was obtained pouring the sand slowly from very low height, using 
the same funnel used for obtaining compact sand. The average dry density, 
porosity, and relative density the sand used this investigation were also 
given Table 


TABLE 1.—Properties Sands 


Angle Shearing 
Resistance, 
degrees 


Dry 
density, 
pounds per Porosity, 
cubic feet asa Relative Plane-strain Triaxial 
Packing (kilonewtons per per- density, test test 
state cubic meter) centage results results 


(1) (3) (4) (5) (6) 


Dense 
Compact 


Loose 


203 mm 
width of the testing box 


FIG. Strip Footing in. 25.4 mm) 


Two main boxes were used this investigation: (1) Rectangular for strip 
footing tests; and (2) circular for circular footing tests. The rectangular box 
was designed simulate the plane-strain conditions strip footing tests. The 
box had glass sides having length, width, and depth in. (610 mm), 
in. (200 and in. (510 mm), respectively. The circular box was 
steel drum in. (480 mm) height and (510mm) diameter. compression 
machine was used apply the load the footings. order maintain the 
load the vertical direction during the tests, the loading ram was passed through 
lubricated ball bearing guide. 

The strip and circular model footings used this investigation in. (51 
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104.0 0.369 0.691 47.7 43.2 
(16.3) 
95.5 0.420 0.465 42.8 42.1 
(15.0) 
87.8 0.467 0.218 34.0 33.6 
(13.8) 
Aluminum 
“32mm 
Sand 
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Theory 


Test results 


bearing capacity , q, (kPa) 


Ultimate 


Depth of dense sand below footing base/footing width, H/B 


FIG. Results: Strip Footing Dense Sand Overlying Loose Sand 
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Test results 


Ultimate bearing capacity (kPa) 


Depth dense sand below footing base/footing 
width, H/B 


FIG. Results: Strip Footing Dense Sand Overlying Compact Sand 
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mm) were machined from aluminum sections. threaded hole 
the center the footing allowed rigid connection loading ram through 
which the loads were applied. The bases the footings were roughened 
cementing fine grain sandpaper onto them using epoxy resin glue. simulate 
plane-strain condition, the strip footing had in. (200 mm) long (width 
the testing box). minimize friction the interface the strip footing 
and the glass sides the testing box, cut 1/5 in. mm) was made 
the footing material both ends and replaced with flexible foam. One side 
the foam was glued the footing, while the other side was covered with 
polyethylene (Fig. 3). 


Test 


The test results are given Figs. and for strip footing and Fig. 
for circular footing. These figures show the variation the ultimate bearing 


a 
~ 
c 
= 
a 


Depth dense sand below footing base/Footing 
H/B 


FIG. Test Results: Circular Footing Dense Sand Overlying Loose 
Sand 


capacity with the ratio the upper layer thickness below the footing base 
the footing width diameter would expected, the observed 
ultimate bearing capacity increased rapidly with increasing thickness the dense 
sand below the footing base from minimum value (at H/B equal 
the ultimate bearing capacity the footing homogenous lower weak sand 
layer maximum value (at equal the ultimate bearing capacity 
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the footing homogeneous upper strong sand layer. Practically the maximum 
value the ultimate bearing capacity achieved value ratio beyond 
which further change the ultimate bearing capacity was observed, 
H/B ratio which the existence weak sand layer does not influence 
the measured ultimate bearing capacity, where the upper layer can accommodate 
classical failure. This ratio depends the relative strength the layered 
system. The lower and upper limits for this ratio are zero for footing 
two sand layers having the same shear strength (case homogenous sand) 
and infinity for footing dense sand layer overlying soil layer hypothetically 
having shear strength. 

The test results the present investigation were analyzed employing 
and The solution obtained follows. 


The left hand side Eq. can evaluated from the test data, again 
assuming value the angle the value can evaluated from Eq. 


TABLE 2.—Comparison between Theoretical and Experimental Average Ratios 


Angle Shearing 
Resistance Average 

(Plane-Strain) Ratio 

case degrees degrees 
(1) (2) (3) 

Loose 47.7 34.0 

Dense /Compact 47.7 42.8 


and thus the right hand side Eq. The solution continued trial and 
error until both sides Eq. were equal. 

The deduced average angle shearing resistance from the present test 
results agreed well with the theoretical values (Table determined from Eq. 
if: 


The coefficient replaced with the ratio which and 
are the ultimate bearing capacities the same footing the surface 
homogeneous lower and upper sand layers, respectively (8). Further 


The distribution the locally mobilized angle shearing resistance 
taken parabola with the following form: 


with the origin the soil interface. The plane parallel the assumed failure 
plane (at distance and the interfacial plane between the two layers 


ore- 
ical 
6) e 
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66,99 


were assigned the ‘‘x’’-axis and respectively, (Fig. 2). Further, 


which depth the upper layer thickness below the footing 
base, and 


should noted that the two aforementioned assumptions reasonably 
represent the actual situation, since the ratio expresses the relative strengths 
the two layers, and the proposed parabolic distribution reflects the difference 
between the assumed failure plane and the actual one, which has direct effect 
the values. 


The computer program referred earlier was used again utilizing Eq. which 
expresses the distribution the locally mobilized angle shear resistance 


resistance the Dense/compact sands 
upper sand layer sands 


shearing resistance the lower 
layer 


FIG. Ratios for Strip Footing Strong Sand Layer Overlying Weak Sand 
Deposit 


924 JULY 1981 


and coefficient equal the ratio g,/q,. The purpose this procedure 
was determine the average mobilized angle shearing resistance for given 
angles shearing resistance the upper and lower and respectively. 
The results this analysis are presented Fig. 

interest note Fig. that the values tend unity when 
which the case homogeneous soil, where the assumed vertical 
planes are real planes failure the radial shear zones (7). Thus, the average 
mobilized angle shearing resistance these vertical planes must equal 


Angle shearing resistance 


~ 
~ 
© 
= 
2 
= 
3 
a 
° 
° 
oO 


shearing resistance the 
lower sand layer 


FIG. 8.—Design Chart for Determining Coefficient Punching Shear Resistance, 


the angle shearing resistance the soil, Setting 


which the coefficient punching shear resistance (10), and substituting 


q w 
u 1 H s B 


Knowing the ratio for given and values (Fig. 7), the coefficient 
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passive pressure corresponding the angle can determined (2). Finally 
the coefficient punching shear resistance can calculated using Eq. 12. 
For convenience, the coefficient punching shear resistance, was determined 
for the ratios given Fig. The results are presented Fig. 
the form design charts. Eq. together with the design charts will allow 
the prediction the ultimate bearing capacity strip footing strong 
sand layer overlying weak sand deposit. Comparison between experimental 
and theoretical values the ultimate bearing capacity are given Figs. 
and 

The analysis for strip footings can extended circular footings 
determining the passive resistance inclined angle vertical cylindrical 
surface through the footing edge (Fig. 1). Thus, Eq. becomes 


tan 


which the shape factor for punching shear resistance cylindrical 
surface. 

The present test results circular footing (Fig. dense sand layer 
overlying loose sand layer were analyzed determine the experimental shape 
factor S,, using the design charts given Fig. The deduced values 
were slightly higher than unity the plane-strain angles shearing resistance 
and were used the analysis. If, however, the triaxial values were 
used, the deduced values were between According the failure 
mechanism observed this investigation for circular footing tests [Fig. 
the shear box test may appropriate for determining the shear strengths 
the sand layer. trial calculation using the angles shearing resistance deduced 


from shear box tests (6, 45.6° and 34°) showed shape factor 
about unity. 


theoretical and experimental investigation the ultimate bearing capacity 
footings strong sand layer overlying weak sand deposit has been 
presented. The theory developed was verified with results tests conducted 
strip and circular footings dense sand overlying loose well compact 
sand. 

From the design chart introduced (Fig. 8), the coefficient punching shear 
resistance can determined from known values the angles shearing 
resistance and the upper and lower layers These angles 
may determined from plane-strain shear test results; however, for conservative 
design, shear box test results may used. The coefficient punching shear 
resistance can used conjunction with the proposed bearing capacity equations 
(Eqs. and 14) predict the ultimate bearing capacity strip circular 
footings. For circular footings, shape factor unity recommended the 
shear strength the sands were determined from either plane-strain shear 
box test results. If, however, triaxial shear tests results were used, shape 
factor two suggested. 

The theory given this paper may also utilized determine the thickness 
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upper sand layer provided artificially increase the ultimate bearing 
capacity desired design value. 


The financial support from the Natural Science and Engineering Research 
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S, 
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relative density sand; 

thickness upper layer below footing base, 

coefficient passive earth pressure; 

coefficient passive earth pressure depth 

coefficient punching shear resistance; 

bearing capacity factors; 

bearing capacity factors for homogeneous upper sand layer; 
bearing capacity factors for homogeneous lower sand layer; 
total passive earth pressure; 

ultimate bearing capacity; 

ultimate bearing capacity footing interface two layer 
soil; 

ultimate bearing capacity footing thick upper sand layer 

ultimate bearing capacity surface strip footing homogeneous 
upper sand layer; 

ultimate bearing capacity surface strip footing homogeneous 
lower sand layer; 

shape factor; 

distance; 

depth point; 

bulk density soil; 

bulk density upper sand layer; 

bulk density lower sand layer; 

parameter (coefficient); 

angle shearing resistance soil; 

locally mobilized angle shearing resistance assumed failure 
planes; and 

average mobilized angle shearing resistance assumed failure 
planes. 
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CIRCLES AND CRITICAL SHEAR PLANES 
Eric 


INTRODUCTION 


Mohr’s circle stress drawn for any point mass soil subject 
stress, the geometry the figure will show, firstly, that there are two planes 
which the shear stress maximum relation the normal effective 
stress and, secondly, that the factor safety, against shearing these 
two planes the same. 

has been suggested that when the stability earth slope determined 
using the Limit Equilibrium Method, potential slip surface might intersected 
series critical shear planes, aforementioned, which the values 
the factor safety would less than that for the assumed slip surface. 
this were true, could inferred that the Limit Equilibrium Method 
overestimates the stability earth slope. 

The risk overestimation would appear particularly high when, 
cross section, the slip surface assumed some simple geometric shape, 


such circular arc. For this reason, was decided that slip circle analyses 
should used check the validity the suggestion. 


The procedure used this investigation can summarized follows: 


Earth slope stability problems were analyzed order obtain: (a) The 
value, the factor safety against failure slip circle; (b) the magnitude 
and position the force due effective stress normal each (vertical) interslice 
boundary; and (c) the normal effective stress and mobilized shear stress 
the base each slice and the slope, the base. 

estimate was made the (horizontal) effective stress the bottom 
each interslice boundary. 

Mohr’s circle was constructed for element adjacent the slip circle 
each interslice boundary order determine the slope, the critical 

Lect., Dept. Civ. Engrg., Sheffield City Polytechnic, Pond Street, Sheffield, 
1WB, England. 

Note.—Discussion open until December 1981. extend the closing date one month, 
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shear planes through the element and the (common) value for these planes. 
Comparisons were then made between the values and and between 
corresponding values and 


Further details will not given the method analysis, the method used 
estimate the distribution stress the interslice boundaries, the determination 
the slopes the critical shear planes, and their factor safety. 


The method used analyze the problems was developed the writer (3). 
resembles the method Morgenstern and Price (1) the following respects: 


The mass soil enclosed the slip surface divided into narrow, vertical 
slices. 

All the relevant conditions equilibrium are satisfied. 

The method applies failure surfaces any shape. 

The method permits variation the slope, the total interslice forces 


using 


Like variable whose value adjusted until the conditions equilibrium 
are satisfied. range values, one for each interslice boundary, chosen 
arbitrarily for the coefficient k,. 

The method uses the same definition for the factor safety, for the 


slip surface, i.e., the factor which the shear strength parameters, 
and tan must divided bring the earth mass enclosed the slip surface 
into state limiting equilibrium. 


The method has been programmed Fortran for solution computer. The 
program includes instructions for large amount data printed out 


The original program gave the position and magnitude the force, Z’, due 
the effective stress normal each interslice boundary. The position was 
expressed terms the ratio L’/H, which the distance the point 
action above the bottom boundary height Knowing these 
values, possible estimate the value the horizontal effective stress 
the bottom boundary assuming the form the distribution stress 
the boundary. first approximation, was assumed that the distribution 
would linear. can then shown that the values the (horizontal) effective 
stresses the top, and bottom, interslice boundary are 
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The derivation these expressions given Appendix 

There however problem L’/H less than one-third—which often 
the case. then found that the stress, the top the boundary 
negative—or tensile. fact, the stress distribution then resembles that for 
structural member subject both compression and bending. 

Few geotechnical engineers would satisfied rely the tensile strength 
soil. One way out the difficulty assume that tension crack will 


A 


s 
= 
= 
2 
2 
Ss 
= 
a 


B 


FIG. 1.—Stress Distribution Interslice Boundary which 1/3 


form the top any boundary where less than one-third and that 
below this crack the pressure distribution will triangular, shown Fig. 
clear from the figure that 


=z 
ray 
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This expression has been used determine the value the (horizontal) effective 


stress the bottom those interslice boundaries which L’/H less than 
one-third. 


sketch element one side the bottom interslice boundary 
shown Fig. 2(a); variables o’, and the slope the base 


the slice, and the normal effective and mobilized shear stresses respectively; 
and the normal effective and shear stresses the bottom the 


Interslice 
boundary 


FIG. Circle for Stresses Element Bottom Interslice Boundary 


boundary, respectively. The Mohr’s circle for the element shown Fig. 
description the construction this figure given Appendix II. 

the factor safety against shearing the most critical planes, F,, were 
equal the value for the slip circle, then the stress circle would tangential 
line PA, which has slope and position determined the values c’, 


SLIP CIRCLES 933 


were extented would intersect the stress circle. The critical shear planes 
will parallel XO, and YO, and the angle between these planes, will 


tan 
F,= 
tan 


which the slope the tangent the circle, PD. 

Instructions determine and print out the values F,, the slope, 
one the critical shear planes, and the angle, between the planes were 
incorporated into the computer program. 


TABLE 1.—Details Problems 


Depth 

tension 

Problem per square crack, 
number meter tan meters distributions 


(1) (5) 


throughout 

throughout 

throughout 

throughout 

throughout 

Varies linearly from 0.94 
toe 1.06 tension crack 

Varies linearly from 1.05 
toe 0.95 tension crack 


Note: each problem: Angle shearing resistance 30°; pore pressure coefficient, 


FIG. 3.—First and Second Slip Circles Analyzed 


11.55 
11.55 

11.55 7.0 

28.87 4.0 

5.77 12.5 

28.87 4.0 

12.5 

35 m 
20.7 n>, 
| 
3 
& 1 
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TABLE 2.— 


Solution 
number 
(1) (3) 
70.6 


& 


a 


ANALYZED 


Details the problems solved are summarized Table 

The first slip circle analyzed was chosen arbitrarily, i.e., was not 
critical slip circle (on which the factor safety, would have had its 
lowest value). The curve marked Fig. The value the pore pressure 
coefficient, r,, was taken 0.5 and the soil properties this uniform earth 
slope were: (1) Effective angle shearing resistance, 30°; (2) effective 


934 
(7) (9) 
69.3 69.8 69.8 70.0 70.3 
1.53 1.57 1.57 1.59 1.61 
0.46 0.36 0.32 0.30 0.28 
68.3 68.3 68.3 68.3 68.3 
1.45 1.45 1.45 1.45 1.45 
0.46 0.39 0.37 0.35 0.34 
68.2 68.2 68.1 68.1 68.1 
1.44 1.44 1.44 1.44 1.44 
0.47 0.41 0.38 0.36 0.35 
62.8 62.8 62.8 62.7 62.8 
1.13 1.13 1.12 1.12 1.12 
0.49 0.41 0.38 0.36 0.36 
73.1 73.1 73.1 73.1 73.1 
1.90 1.90 1.90 1.90 1.90 
0.46 0.40 0.38 0.36 0.35 
63.0 63.0 62.9 62.9 62.9 
1.13 1.13 1.13 1.13 1.13 
0.47 0.39 0.36 0.34 0.33 
73.1 73.1 73.0 73.0 73.0 
1.90 1.90 1.90 1.89 1.89 
0.47 0.41 0.39 0.37 0.36 
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Solutions 


slice Boundary 


Details the solution obtained for this problem are given the first five 
lines the Table The values given the first, second and fourth lines 
indicate that the slip circle intersected series planes which the 
factor safety against shearing, F,, substantially less than the value 
(1.64) for the slip circle. The Mohr’s circle many the boundaries resembled 
that shown Fig. the 12th, 13th, and 14th boundaries, the circle 
actually enclosed point this indicated that the values the stresses 
were inadmissible. 

The problem was then solved again identifying the position the most critical 


935 
(10) (11) (17) (19) 
3.4 22.9 28.1 33.6 39.4 45.8 
1.6 4.0 6.2 7.8 8.2 5.5 73.7 
1.63 1.64 1.63 1.60 1.52 1.16 1.05 
8.9 12.5 23.8 27.7 31.8 36.1 40.6 
5.8 8.6 11.3 13.8 16.4 18.7 20.2 14.3 
1.45 1.45 1.45 1.45 1.44 1.43 1.39 0.96 
0.32 0.31 0.30 0.29 0.24 0.19 0.06 0.30 
3.6 6.9 20.3 23.8 27.4 31.1 35.0 
5.5 8.2 10.9 16.3 19.3 22.7 27.0 33.6 47.1 
68.1 68.1 68.2 68.1 68.2 68.2 68.1 66.3 
1.44 1.44 1.44 1.44 1.44 1.44 1.44 1.44 
0.34 0.33 0.33 0.31 0.31 0.30 0.30 0.37 
6.7 9.7 12.8 22.2 25.4 28.8 32.3 35.9 
9.4 11.7 16.1 18.5 21.0 23.9 27.6 33.0 45.3 
62.8 62.9 63.0 63.0 62.9 62.9 63.0 63.0 61.6 
1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.13 1.07 
0.35 0.34 0.33 0.33 0.33 0.34 0.42 
3.1 6.6 13.6 17.1 20.8 24.5 28.3 32.3 
1.7 5.0 8.1 11.3 14.6 18.3 22.5 28.1 36.7 
73.1 73.1 73.1 73.1 73.1 73.0 72.4 69.5 
1.90 1.90 1.90 1.90 1.90 1.90 1.90 1.89 1.82 1.54 
0.34 0.32 0.31 0.30 0.29 0.29 0.30 0.41 
6.7 9.7 12.8 19.0 22.2 25.4 28.8 32.3 35.9 
7.7 10.0 14.5 16.9 19.4 22.4 26.0 31.6 44.2 
63.0 63.0 63.0 62.9 62.8 62.8 62.8 63.0 61.9 
1.13 1.13 1.13 1.13 1.12 1.13 1.13 1.08 
0.33 0.32 0.31 0.31 0.31 0.31 0.32 0.40 
3.1 6.6 10.1 13.6 17.1 20.8 24.5 28.3 32.3 
73.0 73.0 73.0 73.1 73.1 73.0 72.9 72.3 69.5 
1.89 1.89 1.89 1.90 1.90 1.89 1.88 1.81 1.54 
0.35 0.34 0.33 0.31 0.31 0.30 0.31 0.41 
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slip circle, which shown curve Fig. The solution for this circle 
given Table (solution 2). 

There now much better agreement between the values and indicating 
that one the critical shear planes the bottom the interslice boundaries 
almost tangential the slip circle. Furthermore, except the upper end, 
the agreement between the values and quite remarkable. However, 
even this solution, the stress circles the upper end the slip circle were 
inadmissible. clear from the low values L’/H these points that there 
tension the soil. Further solutions were therefore tried which was 
assumed that the slip circle would terminate, its upper end, the bottom 
vertical tension crack full water. 


Tension Crack 


The depth the tension crack sometimes obtained formula. This involves 
difficulty where the formula includes the shear strength parameters; 
consistent with the other parts the analysis essential that the 
shear strength parameters used. determine the mobilized shear strength 
parameters necessary know the value the factor safety, but 
the determination one the purposes the analysis. One therefore 
the position needing know the answer the problem before attempting 
solve it. 

way out this difficulty has been suggested the writer (4) who showed 
that the depth the tension crack depends upon: the slope, the pore water 
pressure coefficient, and the parameter, tan none these depends 
upon and usually all will known begin with. chart was given 
facilitate the determination the depth the tension crack from the values 
these three parameters. 

alternative method determining suitable depth for the tension crack 
try successive solutions which the depth gradually increased until 
satisfactory solution obtained. This method has been adopted here. passing, 
can said that the inclusion tension crack the analysis (even with 
water pressure it) has remarkable little effect either the value 
the position the critical slip surface. 

the next problem solved, the height, slope, values the soil properties, 
and the position the circle were all the same problem the only 
difference was the inclusion tension crack which was assumed full 
water. The depth the crack was taken, first, This reduced the 
value from 1.45 1.44, but the low values L’/H obtained solution 
were replaced higher values. 

further improvement the values L’/H was obtained increasing 
the depth the tension crack This solution number Tables 
and The values this solution lie between 0.47 and 0.3 implying 
that any tensile stress the soil would low intensity. Furthermore, the 
agreement between the values and and particularly between and 
are now even better than before and none the stress circles was inadmissible. 
the whole, this can said satisfactory solution. should noted 
that the inclusion the tension crack has reduced the factor safety for 
the slip circle less than (1.45 
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point some importance about these solutions that, the absence 
details about the stress circles, would not have been known that the values 
for the factor safety for the slip circle would have been based inadmissible 
stresses some points the soil. 

the aforementioned cases, the (total) interslice forces were assumed 
parallel i.e., the value the coefficient, (in expression 1), was taken 
one throughout. Further problems were then solved for other types 
distribution, though without finding any real improvement solution 


the cases examined now, was taken that the 
value the important parameter tan was 11.55. Other problems 


FIG. Circle for Case which 


were then tackled using higher and lower values for and various types 
distribution and depth tension crack. Details the problems and the most 
satisfactory solutions are shown the tables (see problems 

Reasonably satisfactory solutions, showing good agreement between the values 
and were obtained these problems assuming the interslice forces 
parallel (cases and 5). Some slight improvement the agreement between 
the values and was obtained using the linear distributions described 
Table every case studied, however, the main divergencies between 
and and between and occurred the upper end the slip circle 
where values for less than one-third were obtained indicating the 
presence tension the soil. 

For most the boundaries the solutions shown Table the Mohr’s 
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circles resembled that shown Fig. The principal features this diagram 
are, firstly, that the stress circle tangential lines and PA’—whose 
position and slope are determined the values c’, and Secondly, 
and consequently the angle between the critical shear planes 90° 


The following conclusions can made: 


has been shown that slip circle through earth slope may intersected 
family shear planes which the factor safety against shearing, 
less than that for the slip circle, Where the slip circle chosen 
arbitrarily, the discrepancy may substantial some points. the other 
hand, where the slip circle critical one the difference between the values 
and negligible except, perhaps, the upper end the slip circle 
where tension begins play important part. appears therefore that, from 
this point view, the factor safety for critical slip circle does not seriously 
underestimate the stability earth slope. 

Except the upper end critical slip circle, where there some 
divergence, there close relationship between the slope, any point 
Critical slip circle and the slope, one the critical shear planes. 
The other critical shear planes intersect the slip circle angle which 
almost constant and agrees closely with the value 90° 


remains discovered whether there critical slip surface any 
shape through earth slope uniform soil which the factor safety 
constant and exactly equal the value the critical shear planes which 
intersect it. appears that, such surface does exist, its shape cannot differ 
very greatly from circular arc. 

cannot assumed without proof that the critical shear surfaces which 
intersect the principal slip surface extend planes ground level, seems 
likely that they would curved. Neither should assumed that the factor 
safety against shearing would constant along these surfaces. This raises 
interesting point where the slip surface under consideration not constant 
radius. Any movement along such surface would involve distortion the 
soil mass enclosed it. This distortion could occur only the principal slip 
surface were intersected family shear surfaces which the factor 
safety was constant and equal that for the principal slip surface. 

Cases this type have been analyzed Sarma (2) assuming interslice 
boundaries which, general, were not vertical. trial-and-error identifies 
boundaries which the average value the factor safety equal that 
the principal slip surface. does not, however, appear have proved 
that the factor safety constant along these surfaces. not, then 
reserve strength would remain some points the soil when failure 
was imminent the slip surface. 

linear stress distribution was assumed the interslice boundaries obtaining 
the solutions given Tables and Other types distribution are, course, 
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possible, but were not tried since the more simple assumption has led reasonable 
results. 


The writer wishes express his appreciation Josselin Jong the 
Technical University Delft, The Netherlands, for having drawn his attention 


the presence critical shear planes intersecting slip surface through 
earth slope. 


normal effective stress the boundary and acts height above the 


bottom it; and the values the normal effective stress the 
top and bottom the boundary, respectively. 


B 


FIG. 5.—Stress Distribution Interslice Boundary, AB, when 2/3 


Then equal the area the stress distribution diagram 


Taking moments about the bottom the boundary, 
2 2 


Combining and 6a, obtain 


2 
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, 
, 


Similarly 


Eq. identical with Eq. and Eq. with Eq. 


The procedure followed the construction the Mohr’s circle shown 


The shear strength axis drawn vertically and the normal effective stress 
axis horizontally. Points and are located the horizontal axis from the 
values and which are, respectively, the normal effective stresses 
the base and vertical faces the element. Vertical lines are drawn through 
and 

Point located the vertical line through from the value 
the mobilized shear stress the base (and top) the element. 

The variable the slope the base the element. line slope 


940 
Zz’ , 
Substituting Eq. from Eq. 
2— 
2— 
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drawn through intersect the vertical line through O,. The 
perpendicular biscector AO, constructed intersecting the horizontal axis 

circle drawn with the center passing through and O,. The 
vertical line through intersects the opposite side the circle Variable 
often called the line drawn through this point 
parallel any plane through the element, the coordinates the point where 
this line intersects the opposite side the circle give, respectively, the normal 
effective stress and shear stress that plane. Thus, the coordinates 
correspond the stresses the base (and top) the element and the coordinates 
the stresses its vertical faces. 

The horizontal axis extended the left and the position point 
will found that this line has slope equal 

Line drawn tangentially the circle Point can located 
either drawing second tangent the circle or, Fig. 2(b), from the 
intersection vertical line through with the lower side the circle. Both 
and are joined O,. Lines XO, and YO, give the slopes the critical 


shear planes; the slope one these planes; the angle between them 
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The following symbols are used this paper: 


effective cohesion; 

mobilized cohesion; 

factor safety slip circle; 

factor safety against shearing critical shear planes; 

height interslice boundary; 

value coefficient which variation slope (total) interslice 
forces obtained; 

height force, Z’, due normal effective stress interslice 
boundary; 

force due effective stress normal interslice boundary; 

slope base slice; 

unit weight soil; 


« 
a = 
= 
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siope (total) interslice force; 

angle between critical shear planes; 

variable relating slope interslice forces; 

normal effective stress bottom interslice boundary; 
normal effective stress top interslice boundary; 
shear stress; 

mobilized shear stress; 

shear stress vertical plane; 

effective angle shearing resistance; 

mobilized angle shearing resistance (tan tan 
slope tangent from point Mohr’s circle; and 
slope one critical shear planes. 
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PROPERTIES SOILS FROM TESTS 
Shamsher ASCE and Vijay Kumar ASCE 


INTRODUCTION 


The dynamic stiffness properties soils must determined for the purpose 
analyzing machine foundations and substructures subjected earthquakes. 
The significant advances that have been made recent vears the mathematical 
techniques for analyzing dynamic soil-structure interaction problems necessitate 
realistic determination pertinent soil stiffness properties. This determination 
can quite involved because the dependence these properties large 
number parameters. The stiffness characteristics soil under dynamic loads 
can represented best terms the dynamic shear modulus, which 
may determined laboratory field tests. laboratory test commonly 
used for the purpose the column test.’’ The field tests may include 
vibration tests blocks plates, steady-state vibration Rayleigh wave 
tests, cyclic load tests, and wave propagation tests. 

general, the values dynamic soil moduli determined from laboratory 
and in-situ tests may differ significantly (4,11,19). This may due problems 
associated with sample disturbance, the nonhomogeneity the soil deposits, 
and the effects time. The latter important the case clays, and the 
correlation between the shear moduli determined from laboratory and in-situ 
tests clayey deposits better after accounting for the effects time (1,9,24,25). 
There need account for other factors affecting dynamic shear modulus 
interpreting field test data. Although difficult establish the superiority 
laboratory tests over field tests vice versa, there may situations 
which undisturbed samples may rather difficult not impossible obtain, 
such saturated cohesionless soils below the water table. Thus, greater 
reliance may have placed in-situ testing. The writers have attempted 
determine the realistic values the dynamic shear modulus, from in-situ 
tests accounting for the effect the significant parameters affecting it. 
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The approach suggested this paper makes possible study the variation 
with strains over range interest based upon field test data only. 
The utility the suggested approach has been effectively demonstrated 
analyzing data from different sites India. 

India, one more the following field tests are commonly used 
determine the dynamic stiffness properties soil: (1) Vibration tests blocks 
plates; (2) steady-state vibration Rayleigh wave tests (known shear 
modulus tests India); (3) wave propagation tests; and (4) cyclic plate load 
tests. 

The procedures for conducting these tests are well-known (3,12,27). The Indian 
standard dealing with in-situ dynamic properties soils (7) gives detailed 
procedure for conducting these tests and interpreting the test data. Because 
this standard may not readily available most readers, brief description 
the relevant test procedures and interpretation the test data each case 
included here. 

Block vibration tests are conducted cast in-situ concrete blocks. The block 
size commonly adopted India 1.5 0.75 0.7 high. Variations 
the size the test block may made depending upon the soil conditions 
and type testing equipment available. The block size may thus increased 
the soil hard consists boulder deposits. The object bring the 
observed natural frequency the block-soil system within the range the 
vibration exciting equipment. The sensitivity the recording equipment also 
taken into consideration deciding the size the test block. Because 
the stiffness parameters soil-foundation system may affected the 
shape the foundation contact area, may preferable conduct block 
vibration test circular, square, rectangular block depending upon 
the shape the prototype foundation. The test block essentially surface 
footing cast the proposed depth the foundation. The test conducted 
exciting the block either vertical horizontal mode, and monitoring 
the frequencies and corresponding amplitudes vibration. The natural frequency, 
i.e., the frequency which maximum amplitude occurs the pertinent mode, 
can thus obtained. The test may repeated varying the level dynamic 
force. This done changing the angle setting the eccentric masses, 
which denotes the change distance the center gravity the rotating 
masses from the axis rotation. Dynamic force given rotating mass oscillator 
circular frequency. change represents change effective eccentricity. 

From the observed natural frequency (damped natural frequency), the case 
the vertical vibration test, the value the coefficient elastic uniform 
compression, C,, which defined the ratio stress intensity elastic 
settlement for base given area, can determined treating the 
foundation-soil system single degree-of-freedom system and using 


which mass the block and mounted equipment; observed natural 
frequency vertical vibrations, hertz; and area contact the block 
with the soil. From the value C,, the value Young’s modulus, may 
computed using (3) 


7 
2 


DYNAMIC PROPERTIES 


which Poisson’s ratio; and constant depending upon the ratio 
length width the foundation. Typical values are given Table 
The value may then computed from 


When horizontal exciting force used the block vibration test, becomes 
case coupled rocking and sliding, and the natural frequency vibrations, 


expressed terms the limiting natural frequencies rocking 
and sliding (3) 


2 2 
4 (@ + 


which natural frequency (frequencies) coupled rocking and sliding; 
limiting natural frequency sliding; limiting natural frequency 


TABLE 1.—Values for Rigid and Flexible Footing (3) 


Rigid footing Flexible footing 
(2) 


rocking; and ratio where M,, represents the mass moment 
inertia the foundation about axis through its center gravity and 
perpendicular the plane vibrations; and the mass moment inertia 
the foundation about axis through the projection the center gravity 
the base contact area, and perpendicular the plane vibrations. The 
foundation. 
The values w,, are given (3) 


which coefficient elastic uniform shear; coefficient elastic 
nonuniform compression; and moment inertia the contact area 
foundation with the soil about the axis vibration. 

The terms C,, C,, and are related follows (3) 


945 
2(1 
(1) 
1.08 1.06 
1.10 1.09 
1.15 1.13 
1.24 1.22 
1.41 1.41 
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and 


From the field test data, the values observed natural frequency combined 
rocking and sliding are known, and the values w,, can expressed 
terms making the necessary substitutions Eq. one can obtain 


— 4+ = 0 @ Wie (9) 
m M ano ¥ ym M,,, 


rearranging Eq. and solving for one obtains 


The plus sign used the observed natural frequency, corresponds 
first mode vibration, and the minus sign used corresponds 
the second mode vibration. From the value thus obtained, and 
thus, and can computed. 

The vibration tests (vertical vibrations only) may conducted plates 
manner similar that blocks. mechanical oscillator attached 
the plate, and the mass the system can varied adding weights. The 
method conducting the test and interpreting the test data the same 
described earlier for the block vibration tests for vertical vibrations. 

free vibration test block may conducted exciting the block 
free vibrations either hitting it, pulling and suddenly releasing it. The 
vibrations the block are then recorded, and the natural frequency obtained 
directly from the record (12,18,21,25). The values the dynamic soil properties 
can then evaluated using Eqs. 1-8, depending upon the particular case 
vibration being considered. 

The steady-state vibration test (shear modulus test) conducted exciting 
base into vertical vibrations known frequency and finding the phase 
difference between the arrival the Rayleigh waves two points known 
distance apart (7,22,27). knowing the phase difference and the distance between 
the points observation, the wave length the Rayleigh waves can 
determined. degrees the phase difference between the arrival the 
Rayleigh waves the two observation points distance, apart, then 
the wave length, given 


The velocity the Rayleigh waves, may obtained from 


The Rayleigh wave velocity, and shear wave velocity, are nearly 
the same for the Poisson’s ratio between 0.25 and 0.5 (for 0.5, 
0.955 and for 0.25, 0.919 (22). 
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TABLE 2.—Soil Description Different Sites Investigated 


Project 
description 
(2) 
Forging hammer 
foundation 
HAL, Koraput 


Aero-engine test 
beds AEF, 
Chandigargh 


Diesel power- 
house, Sirhind 


Diesel power- 
house, Nakodar 


Cement factory, 
Rajban 


University 
Roorkee, Roor- 
kee 


Note: void ratio; and angle internal friction. 


Soil 
characteristics 
(3) 
Silty sand (lateritic soil) 
2.67 
1.93 g/cm’ 
0.59 
15.2% 
Silty sand 
2.61 
0.72 
30.5 
18.6% 
Fine medium sand with 
some silt 
1.75 
0.72 
30° 
22.2% 
Medium sand 
2.58 
1.79 
0.718 
19.2% 


Boulder deposits with ma- 
trix medium coarse 
silty sand; properties 
the matrix material 

2.70 

2.15 

0.605 

32.5° 

28.0% 

Poorly-graded fine silty 
sand 5.0m 

2.63 g/cm’ 

1.83 

0.69 

30° 

17.6% 

Depth water table 4.2 
layer silty clay 
17.75 


Types tests 
conducted 
(4) 

Forced vertical vibration 
test; shear modulus (Ray- 
leigh wave) test; and cy- 
clic plate load test. 


Forced vertical vibration 
test; forced horizontal 
vibration test; and cyclic 
plate load test. 


Forced vertical vibration 
test; forced horizontal 
vibration test; and cyclic 
plate load test. 


Forced vertical vibration 
test; forced horizontal 
vibration test; shear 
modulus test (Rayleigh 
wave test); and cyclic 
plate load test. 

Forced vertical vibration 
test; forced horizontal 
vibration test; free verti- 
cal vibration test; and cy- 
clic plate load test. 


Forced vertical vibration 
test; forced horizontal 
vibration test; free verti- 
cal vibration test; and 
free horizontal vibration 
test. 
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The value may therefore obtained from: 


which mass density the medium. 

The wave length the Rayleigh waves depends upon the frequency 
excitation. The depth penetration waves, i.e., the depth the sampled 
soil zone, which depends upon the wave length Rayleigh waves, may thus 
altered changing the frequency the vibrator (22). Further, the average 
soil properties determined may considered representative 

wave propagation tests, the seismic waves are generated either impact 
blasting with small charges, and the time taken the waves travel 
known distance measured (2,10,26). The velocity either the compression 
waves, V,, the shear waves, may thus obtained. case has 
been observed, may computed follows: 


From may determined using Eq. has been measured 
the field, may determined using Eq. 13. 

Cyclic plate load tests are conducted with plate load test device. static 
load first applied, and the total settlement noted. The load then released, 
and the elastic recovery noted. The load increased increments, and the 
process repeated. The value determined the ratio load intensity 
corresponding recoverable settlement (assumed elastic recovery) per 
load increment. The test conducted the proposed depth the foundation, 
Circular plates can also used. 

The writers have had experience determining dynamic soil moduli under 
different soil conditions number sites India (13-17,21). most 
the sites, several tests different types were conducted. The analysis this 
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FIG. 1.—Typical Grain Size Distribution Plots 
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paper confined only those sites where the soil was sand silty sand. 
The soil conditions occurring each the sites are described Col. 
Table Fig. shows typical grain size plots for some sites. Information pertinent 
the tests conducted listed Cols. and Table The instruments 
used the tests are briefly described Appendix 


INTERPRETATION Test 


The test data obtained from the different field tests were analyzed determine 
the value the dynamic shear modulus using the pertinent procedure for 
each case described earlier. The values the dynamic shear modulus obtained 
from the test data are listed Col. Table from which may observed 
that the values any one site generally show wide variation. the 
block vibration tests, the value decreases with increase the angle 
the setting the eccentric masses, i.e., with increase dynamic 
force levels. The dynamic force given rotating mass oscillator 
frequency. change indicates change Free vibration tests and 
steady-state vibration (Rayleigh wave) tests give higher values compared 
those obtained from block vibration tests. The values obtained from cyclic 
load tests also show similar variation relation the results from other tests. 

designer thus faced with the problem choosing appropriate parameters 
for analysis and design. The wide scatter the data poses large problem. 
The question is, should choose the lowest value the average value, 
should guided some other considerations? Stiffness properties the 
soil are the basic input parameters, and the accuracy achieved using any 
advanced technique analysis can better than the precision with which 
these parameters are chosen. must mentioned here that, contrast 
static problem, selecting low values dynamic soil constants may not always 
mean error the safe conservative side. illustration, machine 
foundation has been designed that the natural frequency the soil-foundation 
system higher than the operating speed the machine, the selection 
value for that lower than the actually occurring value, the result 
actual response will better than the predicted response. However, the 
machine foundation designed based erroneous value corresponding 
natural frequency lower than the operating frequency the machine, 
the predicted response the foundation will smaller magnitude than 
the actual performance. The results will still more complicated analysis 
deeply embedded structure, such well foundation reactor building 
under seismic loading conditions, investigated. There is, therefore, necessity 
choose rational value that takes into account the various factors 
affecting this parameter, namely the mean effective confining pressure, void 
ratio, shear strain amplitudes, degree saturation, time effects, and number 
cycles loading per second (5,22). Some these parameters may significant 
certain cases only, e.g., time effects are important case clays undergoing 
secondary consolidation. Parameters, such confining pressure, shear strain, 
amplitude, and void ratio, are far the most important all cases. 

When field tests are conducted given location, void ratio and degree 
saturation may assumed more less the same. should recognized 
that when different types tests are conducted location certain factors 
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TABLE 3.—Values Dynamic Shear Moduli from Field Tests Different Sites 


Dynamic 
shear 
Dynamic modulus, 
shear 
modulus, kilonew- 
tons per 
kilonew- square 
tons per meter 
for 
Size meter confining Associated 
block pressure strain 
(3) (4) (5) (6) (7) 
Forced Block, 1.5 1.5 
1.26 
Shear Block, 1.5 
high 
load test 30.5 1.3 
6.5 
Forced Block, 1.5 
test high 
Forced Block, 1.5 5.8 
tion test high 
Cyclic Plate, 30.5 1.63 
plate 
load test 30.5 
Forced Block, 1.5 8.6 
Forced Block, 1.5 4.2 
Cyclic Plate, 30.5 1.8 


Site 
number 
(1) 
7 
load test 30.5 


Forced 
vertical 
vibration 
test 

Forced 
horizon- 
tal vibra- 
tion test 

Shear 
modulus 
test* 


Cyclic 
plate 
load test 

Forced 
vertical 
vibration 


tion test 
Free verti- 

cal vibra- 

tion test 


Free hori- 
zontal 
vibration 
test 

Cyclic 
plate 
load test 


Forced 
vertical 
vibration 
test 


Forced 
horizon- 
tal vibra- 
tion test 
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TABLE 3.—Continued 


Block, 1.5 
0.75 
0.70 
high 

Block, 1.5 
0.75 
0.70 
high 

Block, 1.5 
0.75 
0.70 
high 

Plate, 30.5 
30.5 

Block, 3.0 
1.5 
1.0 
high 

Block, 3.0 
1.5 
1.0 
high 

Block, 3.0 
1.5 
1.0 
high 

Block, 3.0 
1.5 
1.0 
high 

Plate, 


Block, 1.5 
0.75 
0.7 

Block, 


Block, 1.5 
0.75 
0.70 


1,310.8 
1,255.0 
1,143.0 
1,082.8 
1,498.0 
1,688.0 
1,562.0 
1,480.0 
1,674.0 
1,498.0 


2,443.0 
2,120.0 


1.1 
2.2 
2.8 
3.0 
3.6 
7.7 

9.8 


4.6 
3.06 
3.6 

1.06 


1.03 
1.1 
1.24 
1.45 
1.58 

1.3 
1.5 

4.6 
7.3 
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162.6 448.4 35° 
141.2 390.0 70° 
132.4 365.0 105° 
122.9 122.9 140° 
135.0 371.5 70° 
126.1 347.9 105° 
125.0 344.3 140° 
312.0 867.3 
198.3 166.0 
579.9 RBI 
574.9 RBI 
521.6 RB2 
test 494.0 RB2 
Forced 683.3 
horizon- 770.0 RBI 
tal vibra- 697.7 RB2 
691.2 RB2 
762.1 RBI 
683.3 RB2 
1,114.3 RBI 
967.1 RB2 
232.3 155.0 
204.4 136.3 
297.3 198.4 VP3 
223.0 148.8 
293.0 805.5 35° 
251.0 692.1 70° 
191.8 529.2 105° 
174.1 480.4 140° 
206.2 538.1 35° 
180.8 472.0 70° 
145.0 372.0 105° 
129.0 322.3 140° 
274.3 757.0 35° 
225.0 621.0 105° 
210.0 578.2 140° 
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Block, 


Free verti- Block, 1.5 
high 
high 
Free hori- Block, 1.5 2,001.0 
zontal 0.75 
vibration 0.7 
test high 
Block, 
high 


“Rayleigh wave test. 
angle setting concentric masses. 
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FIG. 2.—Shear Strain Levels Associated with Different In-Situ and Laboratory Tests 
(6) 


come into the picture and need considered. They are: (1) The depth 
the soil sampled; (2) strains induced the soil; and (3) significant variations 
local soil conditions, any. The conditions confining pressure occasioned 
applied loads and the depth the soil sampled are obviously different 
for different tests. For example, the effective depth below the surface for which 
determined the block vibration test and the wave propagation test are 
quite different, and accordingly the mean effective confinement for which 
has been determined these tests quite different. Similarly, the shear strains 
induced the soil different tests are quite different. Fig. shows the order 
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the shear strains induced soil when different laboratory and in-situ 
dynamic tests are conducted (6). 

previously examined, the confining pressure and shear strain level are 
different different tests, and because these parameters exert the largest 
influence the value this believed the main reason for the 
difference the values the dynamic soil constants obtained site 
different methods, and should taken into account choosing the value 
the dynamic shear modulus. This necessitates study how varies with 
confining pressure and shear strains. The shear modulus, has been found 
vary with confinement 


which dynamic shear modulus average effective confining pressure 
dynamic shear modulus average effective confining pressure 
and index power. The average effective confirming pressure, 


which effective vertical stress the point under consideration. Eq. 
true for sands and normally consolidated clays, and average value 
0.5 may used (23). 

The shear strains induced soil are not precisely known. the case 
wave propagation tests, the shear strain amplitudes are low and are assumed 
the order 10~°. The shear strains induced soil essentially depend 
upon the amplitude vibration settlement, which turn depends upon 
superimposed forces, the foundation contact area, and soil characteristics. The 
measured values amplitude settlement take care the factors affecting 
them. The shear strain amplitudes, y,, may taken equal the ratio 
the amplitude width (19,20). For the case vertical vibration 
tests, the ratio amplitude/width will give the normal strain, and the 
shear strain, y,, may obtained 


(18) 


(19) 
(20) 
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the loading condition corresponds the rest condition case’, then 

one considers extreme case incipient failure with the major 
principal stress 


and the value will depend upon the value the angle internal 
friction the soil, and the Poisson’s ratio, Typical values 
are listed Table 


TABLE 4.—Typical Values 


Angle internal 


For values and the range interest, reasonable assume 
that 

The shear strain from cyclic plate load tests may similarily obtained 
the ratio settlement the width the plate. The free vibration tests will 
induce shear strains, which will low compared those induced forced 
vibration tests resonance, and will higher than those the case wave 
propagation tests. The shear strains for these tests are taken 

the case block vibration tests, the mean effective confining pressure 
which corresponds with the value may taken the effective confinement 
that occurs depth equal width/2 below the base the block because 
the combined effect superimposed static loads the soil and the overburden. 
the case steady-state vibration tests, the effective confining pressure 
the case cyclic load tests, the confining pressure depth equal half 
the width the plate may taken also for block vibration tests. the 
case wave propagation tests, the average depth for which the measured 
wave velocity give representative value varies from L/8-L/10, which 
the distance between the point generation the waves and the observation 
point. The effective confining pressure can determined accordingly. 

using the aforementioned procedure, the values the effective confining 
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0.9317 1.014 1.35 
1.054 1.183 1.398 1.5 
1.097 1.216 1.41 
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pressure and shear strain the soil for which the test gives representative 
values were determined field tests for each case. using 0.5 
Eq. 15, the values obtained the tests were adjusted confining 
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FIG. 3.—Dynamic Shear Modulus, Versus Shear Strain Amplitude, (Site 
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FIG. 4.—Dynamic Shear Modulus, Versus Shear Strain Amplitude, (Site 


pressure 100 Any value mean effective confining pressure may 
used reference value for data reduction purposes and study the 
variation with shear strain amplitudes. The value finally adopted 
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analysis design must invariably corrected for the values confinement 
occuring given problem using Eq. with appropriate value 
The corrected values (reduced confining pressure 100 
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FIG. Shear Modulus, Versus Shear Strain Amplitude, (Site 
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FIG. 6.—Dynamic Shear Modulus, Versus Shear Strain Amplitude, (Site 


are listed Table Col. These values the dynamic shear modulus, 
versus the shear strain, y,, have been plotted Figs. 3-8, which curves 
are drawn for each site. was observed that regular variation with 
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FIG. 8.—Dynamic Shear Modulus, Versus Shear Strain Amplitude, (Site 
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strains obtained. Scatter the field data unavoidable due local variation 
soil conditions. This may construed added advantage these methods 
since the deviation from the average values also becomes known the designer. 
cases very large sites, may desirable subdivide the area into more 
than one zone necessary, provided economy and other considerations justify. 
will also noted from Figs. 3-9 that values for certain ranges 
strain were not obtained during the field tests. This may normally the case 
most the field investigations this type. However, more data may result 
only marginal readjustment the average curves and may not great 
practical significance. 

the case sites and where low strain values (y, 
were available from either the free vibration tests steady-state vibration 
tests, the ordinates the corresponding curve were divided the value 


RANGE FOR SANDS 


(SEED AnD HOR ISS, 1970) 


SHEAR STRAIN AMPLITUDE 


FIG. Versus Shear Strain Amplitude 


Fig. For site versus plot was extrapolated get the value 
shown the dotted line Fig. and the values 
with were obtained and plotted Fig. the case sites and 
Because the values low strain were not available sites and 
free vibration tests steady-state vibration tests wave propagation 
tests were not conducted, the variation G/G,,,, was not calculated. 
The curves obtained Seed and Idriss (23), which represent the range 
variation versus for sands based upon laboratory tests, are also 
shown Fig. can observed that the plots versus for 
the sites under consideration show good agreement with the laboratory test 
data. The differences are essentially because the test data Seed and Idriss 
(23) are for clean sands, whereas the soil the sites investigated the writers 
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generally consisted silty sands. felt, therefore, that the writers’ suggested 
approach, which based entirely upon field test data, may used study 
the variation with shear strains, and the values the shear modulus 
can chosen consistent with the shear strains induced the soil for the prob- 
lem under consideration. The values may then adjusted for the confine- 
ment effects corresponding the problem being investigated. 


The most important point needing consideration this stage the types 
tests that must conducted given site that one can study the variation 
with shear strains based upon in-situ test data and determination whether 
necessary conduct all the tests, namely, the wave propagation tests, 
steady-state vibration tests, block vibration tests, and cyclic load tests. 
following the suggested approach, has been observed that versus shows 
trend variation similar that obtained from laboratory test data. This, 
however, does not imply that only such tests the wave propagation tests 
should conducted and that the variation with should obtained 
comparing with the curves Seed and Idriss (23) for sands clays 
the case may be. The obvious reason that the variation with 
affected the soil type and this fact can established better investigation 
than estimation. The wave propagation test yields valuable information about 
the low strain moduli (y, and best suited for the purpose. 

For the purpose obtaining complete variation over the shear strain 
range interest, the data from the wave propagation test need supplemented 
other tests. Block vibration tests will, general, yield values soil moduli 
the range and cyclic loading tests can conducted obtain 
the values soil parameters strains larger than Therefore, conducting 
wave propagation tests, block vibration tests, and cyclic plate load tests, and 
following the suggested approach, variation over wide shear strain range 
can determined. Free vibration tests and steady-state vibration tests may 
also used wave propagation tests are not conducted. When block vibration 
tests are conducted site, the steady-state vibration tests and free vibration 
tests can made with nominal additional effort. 


approach has been suggested for determining the strain-dependent values 
dynamic shear modulus from field tests. The method suggested has been 
found rather suitable for cohesionless soils below the water table where 
difficult obtain undisturbed samples. The values picked directly 
from field tests can correct the confining pressure and strain amplitudes 
are correctly assessed. Otherwise, they will erroneous. The suggested approach 
attempt account for the effect significant parameters affecting the 
values when field test data are interpreted and may help arriving 
the realistic values for dynamic soil structure interaction. 
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Block Vibration Tests 

Exciting rotating mass mechanical oscillator operated 
5-hp compound-bound d-c motor. The speed the motor was varied 
with the use speed control unit matching capacity that also converted 
220 cps a-c supply into d-c supply. The oscillator motor assembly was 
mounted the top the block give sinusoidally varying force the desired 
direction. 

Transducers.—The response the block was monitored using Miller 
accelerometers with frequency response from The accelerometers 
were mounted the block with the sensitivity axis oriented appropriately 
relation the direction the vibrations being monitored. 

Recording Equipment.—The signal from the accelerometer was amplified 
through Brush universal amplifier and recorded Brush ink writing 
oscillograph with flat frequency response. 


Free Vibration Tests 


Exciting Equipment.—The block was excited into free vibrations striking 
the appropriate face the block with hammer. 

The monitoring and recording the vibrations were done manner similar 
the preceding for the block vibration tests. 


Steady-State Vibration Tests 

Excitation.—The block was excited the motor oscillator assembly described 
earlier that vibrated vertically particular constant frequency. 

Monitoring Wave Arrivals.—The wave arrivals were sensed with the help 
two geophones fixed known distance apart. The output the geophones 
was amplified through high gain d-c amplifiers, and recorded dual-channel 
ink-writing oscillograph through which paper moved constant speed 
125 mm/s. From this, the time interval could measured accuracy 


Cyclic Loading Tests 

Loading Arrangement.—The load the plate was applied through the reaction 
remotely controlled 50-ton capacity hydraulic jack, which acted against 
loaded platform. 

Measurement Settlement.—The settlement was recorded each the four 
corners the test plate with the help mechanical dial gages having accuracy 
count 0.01 mm. 
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The following symbols are used this paper: 


base area foundation; 

width foundation; 

coefficient elastic uniform shear; 

coefficient elastic nonuniform compression; 

distance between geophones; 

Young’s modulus; 

natural frequency, hertz; 

limiting natural frequency sliding, hertz; 

limiting natural frequency rocking, hertz; 

natural frequency vertical vibrations, hertz; 

dynamic shear modulus; 

dynamic shear modulus mean effective confining pressure 

dynamic shear modulus mean effective confining pressure 

dynamic shear modulus shear strain 

specific gravity particles; 

moment inertia base contact area about axis rotation; 
constant; 

wave length Rayleigh waves; 

mass moment inertia about axis through center gravity 
and perpendicular plane vibration; 

mass moment inertia about axis through center gravity 
base and perpendicular plane vibration; 

mass foundation and accessories; 

constant; 

velocity compression waves; 
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velocity Rayleigh waves; 
velocity shear waves; 
ratio 
bulk density soil; 
shear strain; 
normal strain; 
Poisson’s ratio; 
mass density soil; 
three principal stresses given point along axis and 
respectively; 
mean effective confining pressure; 
angle internal friction; 
phase angle; 
natural water content soil; 
natural circular frequency, radians per second; 
natural circular frequency sliding, radians per second; and 
natural circular frequency rocking, radians per second. 


963 


& 


JULY 1981 


CORRELATIONS FOR PILES SAND 


This paper summarizes the results study made improve the present 
state-of-the-art for design piles driven into cohesionless soils. Data obtained 
from full-scale load tests were used evaluate bearing capacity factors. These 
factors were then correlated with parameters which include depth pile 
penetration, pile diameter, and friction angle relative density the sand. 
Special attention was devoted the problem residual stresses for piles 
sand. The failure load from tension test was used represent the side load 
compression adjusted for residual stresses. New design correlations were 
developed which can used predict the bearing capacity piles sand. 


The static formula method (27) commonly used determine the ultimate 


axial bearing capacity pile, shown Fig. For this paper the 
equation and terminology used are 


2. = Q, Q, q.A, +f,A, 


which ultimate point load; ultimate side load; ultimate 
unit point resistance; ultimate unit side resistance; area pile point; 
and area pile shaft. 

Point Load.—The theoretical determination the point load has received 
extensive attention through the years. According (34) the theoretical 
approach solve this problem was initiated Caquot and Buisman the 
mid-1930s. Caquot and Buisman extended the classical work punching failure 
done Prandtl and Reissner approx earlier. Following the same basic 
approach, several different solutions were subsequently presented with different 
assumptions concerning the failure pattern, new empirical corrections were 
introduced. Among the many coniributors this area study were Terzaghi, 
Beer, and Jaky the 1940s; Meyerhof, Caquot, and the 1950s. 

Civ. Engrg., Texas A&M Univ., College Station, Tex. 77843. 

Assoc. Prof. Civ. Engrg., Universidade Federal Santo, Vitoria, Brazil. 

Note.—Discussion open until December 1981. extend the closing date one month, 
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somewhat different approach was taken Skempton, Yassin, and Gibson 
the early 1950s when they dealt with the soil failure induced the pile 
point special case the expansion cavity inside solid. This same 
idea was used again 1977 (37). latest approach involves 
the use complex soil parameters. 

all the theoretical solutions, the ultimate unit point resistance, 
given 


which effective unit weight the soil the pile point; smallest 
foundation dimension; effective overburden pressure the pile point 


Q = ULTIMATE BEARING CAPACITY 


Q, SIDE LOAD 
f, * UNIT SIDE RESISTANCE 
A, = SHAFT AREA 


FOR SAND: 
BK tane 


Py Ny 


UNIT POINT RESISTANCE 
* POINT AREA 


FIG. 1.—Ultimate Bearing Capacity—Static Method 


bearing capacity factors usually depending upon the soil friction angle 
and the assumed pattern mechanism failure. 

engineering practice Eq. usually simplified. the first place, since 
the study restricted sands (cohesionless soils), the second term the 
equation can eliminated. Secondly, when comparing the two remaining terms, 
seen that the first term relatively small and can neglected. The final 
simplification concerns the form the remaining term. Since most piles have 
circular square cross sections and the shape factor the same (14,33) 
reasonable use new bearing factor, that incorporates this constant 
shape factor. Therefore, the commonly used form (5) Eq. becomes 


most the theories the basic parameters, addition the pile geometry, 
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are the friction angle, which used determine the bearing capacity factor, 
N,, and the effective confining pressure the soil. theory considers the 
side resistance the soil along the pile shaft, possible interdependence 
between side and point resistances. Some the theories consider the soil 
rigid-plastic body (1,19,30) while others consider the soil compressible 
(3,37). 

All the bearing capacity theories require the evaluation for use 
Eq. summary the range values according the different 
theories presented Table evident that there are major deviations 
from one theory another, leading the conclusion that the true failure 
mechanism not, generally, well understood. Vesic 1977 (37) stated: 
the ultimate load quite difficult and general solution 
not yet and added: view the many uncertainties involved 


TABLE 1.—Bearing Capacity Factors for Deep Foundations (35,37) 


Approximate Values for 
Various Friction Angles, degrees 
Theories 
(1) 
Beer (1945) 
Meyerhof (1953): 
Driven piles 
Caquot-K érisel (1956) 
Brinch Hansen (1961) 
Skempton-Yassin-Gibson 
(1953) 
Brinch Hansen (1951) 
Berezantsev (1961) 
Vesic (1963) 
200 
Terzaghi (1943): 
General shear 
Localized shear 


analysis pile foundations, has become customary, and many cases 
mandatory, perform certain number full-scale pile load tests the 
site more important 

Side Load.—The theoretical determination side resistance for piles sands 
has received little attention the literature and only few methods (20,22) 
have been published. The interaction between the soil and the pile very 
complex and poorly understood. The determination the ultimate unit side 
resistance, f,, based the laws mechanics considering friction between 
solid surfaces. Fig. shows the system stresses considered. The magnitude 
commonly determined using 


which lateral earth pressure coefficient; average effective overburden 


73.3 
35.1 
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pressure along the segment pile being analyzed; and tan coefficient 
friction between the pile and the soil. The average effective overburden 
pressure, calculated from the knowledge the location the ground-water 
level (usually static) and the magnitude the unit weight the sand. Factors 
and tan need established order determine unit side resistance. 

Potyondy (23) determined the laboratory the coefficient friction between 
various construction materials and cohesionless soils different densities. This 
determination was made using direct shear tests. addition the determination 
the coefficient friction, tan between the soil and the pile material, 
the angle friction the sand, was also determined. Consequently, 
was possible develop relationship between and 

(37) proposed different approach for the determination tan The 
sand located the interface between the soil and the pile considered 
state ultimate failure for determination side resistance. Consequently, 
the angle friction between the pile and the soil, independent the 


soil element 
lundist urbed) 


FIG. 2.—Side Resistance Along Shaft 


initial soil density and pile material. Therefore can considered equal 
the residual friction angle the sand, The difference between the values 
determined Potyondy (23) and those suggested Vesic (37) not seem 
significant. 

Meyerhof 1959 (20) and Nordlund 1963 (22) dealt theoretically with 
the problem the determination the lateral earth pressure coefficient, 
The assumption was made both studies that the pile displaces the sand 
horizontal direction, without any vertical deformation. This displacement 
induces compaction the surrounding soil which maximum the pile-soil 
interface. Since the wall pushes against the sand and the horizontal movement 
large (equal the pile radius), theoretically possible for the magnitude 
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the lateral earth pressure coefficient, high the passive earth 
pressure coefficient. 

list typical values lateral earth pressure coefficients for pile foundations 
compiled Kezdi (14) presented Table should noted that, with 
the exception Ireland who presents range values, only one value 


TABLE 2.—Typical Values Earth Pressure Coefficients (14) 


Basis relationship Soil type Values 
(2) (3) 


Brinch Hansen theory sand 
Lundgren (1960) pile test sand 
Henry theory sand 
Ireland (1957) pulling tests sand 
Meyerhof (1951) analysis field data loose sand 


dense sand 
Mansur-Kaufman analysis field data silt 


(1958) (compression) 
0.6 


(tension) 
Lambe-Whitman 


(1969) 
Kézdi (1958) granular 


suggested each author; and consideration given parameters defining 
the soil-pile system such pile penetration. 


has been shown that the variation predicted values and 
wide that the choice one theory over another difficult exercise 
engineering judgment. The need for better understanding failure mechanisms 
and development acceptable theory has resulted number experimental 
investigations. 

Mayer and L’Herminier 1953 (17) reported work done others and 
recognized three distinct patterns for sand failure according depth penetra- 
tion. These patterns are general shear failure for very shallow foundations, 
localized shear failure for deep foundations, and Meyerhof type failure 
for intermediate penetrations. 1957 Tomlinson (31) noted that these failure 
patterns, with soil being displaced upwards and away from the pile, indicated 
tendency form gap between the pile shaft and the soil. Consequently, 
the active case lateral earth pressure could develop. 

1959 Meyerhof (20) indicated that, cohesionless soils, the effect pile 
driving general compaction resulting from permanent rearrangement and some 
crushing the soil particles. addition, the driving the pile the soil 
mass alters its state stress. Subsequently, 1966, Broms (4) reported 
the zones which undergo general compaction and alteration from the original 
stress. 


Measurements made Széchy 1961 (28), using laboratory model tests, 


[ 
0.8 
K, 
1.75 
0.5 
1.0 
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indicated that the idea general sand compaction not valid. Széchy found 
that highly compressed stress concentration areas are surrounded with 
areas stress-relief below the pile toe well around the pile shaft.’’ 
1964 Robinsky and Morrison (24), with more sophisticated measurements made 
loose and medium dense sands the laboratory, confirmed the earlier findings. 
According Robinsky and Morrison the over-all process sand compaction 
results seemingly erratic accumulation high and low density areas. 
thin sleeve loose sand created around the pile wall, which encircled 
cylinder dense sand. This cylinder, arching, prevents the development 
full lateral earth pressure the pile. Also, 1964, Robinsky and others 
(25) reported additional experimentation which indicated that with pile 
penetration the ‘‘arches’’ around areas loosened sand would built 
and broken down, inducing continuous changes mobilized side resistance. 

Kérisel 1961 (13) published noteworthy paper concerning laboratory 
experiments sand. was observed that different size piles the same sand 
attained maximum value unit load resistance which remained constant with 
increasing penetration. The depth where the constant value was attained varied 
with the pile diameter. The larger the pile diameter the deeper the required 
penetration. 

(36) 1970 confirmed, with field tests, the tendency for unit resistances 
increase with depth some limiting value. Vesic noted that even though 
the rate increase sharply decreases some depth, there was 
additional increase with further penetration. This depth was defined 
being between pile diameters for loose sands and for denser sands. 
1977, Meyerhof and Valsankgar (21) obtained additional laboratory evidence 
limiting value for unit point resistance. The study also showed that the 
depth for submerged sands 1.6 times greater than that for dry 
sands. This increased depth probably caused buoyancy effects. 

Biarez and Gresillon (2) 1972 reported laboratory experiments performed 
Grenoble, France. Pile models various diameters were penetrated into 
metallic rollers bidimensional problem), into glass spheres same diameter, 
and into sand. The sand was tested various densities and with different 
surcharges applied means air pressure (with membranes covering the 
sand). Once again limiting unit resistances were obtained. Also, was concluded 
that the compressibility the soil very important property the determination 
the bearing capacity deep foundations. Another important factor relation 
compressibility the crushing the sand grains noted DeBeer (6). 

The results experimental studies have indicated that certain parameters 
control the pile-soil behavior the pile tip. These parameters are: effective 
overburden pressure (related depth); soil displacement (related pile diameter); 
and sand relative density friction angle soil property). addition, the 
compressibility the soil related the field state stress influences the 
pile-soil behavior. this study, the compressibility the soil taken 
account empirically through combination the aforementioned parameters. 


Test INTERPRETATION 


Since experimental results are qualitative, there need for verification 
significant parameters interpretation full-scale field load tests. Data 
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from field load tests can used establish the significant parameters and 
develop appropriate correlations. For correlation purposes, necessary 
that the load tests sands conducted instrumented piles allowing 
measurement and separation the total load into its components: side load 
(distributed along the pile shaft), and point load (concentrated the pile point). 

The total force the top the pile usually applied means jack. 
This force can accurately measured means load cell. The force 
and the corresponding measured settlement the top the pile permit the 
development load-settlement curve shown Fig. 3(a). 

The distribution the total force between the pile point and the shaft can 
obtained means strain gages load cell placed close the pile 
tip. additional strain gages are placed along the pile shaft the side load 
distribution can calculated also, shown Fig. The difference 


Total Load, Q;. (tons) 


c 

a 

a 


strain gages 


FIG. Instrumented Pile Load Test Results: (a) Load Test; (b) Load 
Distribution in. 25.4 mm; 0.305 ton 8.90 kN) 


load determinations between two strain gages any pair measuring devices 
gives the load being transmitted friction from the pile the surrounding 
soil that pile section. From the strain gage, other measuring device, placed 
close the tip possible obtain both total side load and total point 
load. Then using Eqs. and and knowing the pile geometry, possible 
determine both the unit side and point resistances. 

When analyzing field load tests sand conducted different sites important 
problem that arises the definition failure load. Vesic (37) lists different 
criteria for definition failure load. For this study the Vesic (37) suggestion 
taking the failure load being the one which induces total settlement 
[see Fig. 3(a)] equal 10% the point diameter was used. cases where 
load transfer curves [see Fig. 3(b)] were not given for the established failure 
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Q 
Load Pile, (tons) 
200 100 200 
04 4 10 
8 \ > 20 
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bile 
a) b) 


Pile test 
(1) 
Arkansas 


Arkansas 


Arkansas 


Arkansas #10° 


Arkansas 


Arkansas 


Jonesville 


Jonesville 


Jonesville 


Low-Sill 


Low-Sill 


Low-Sill 


H-14 


Pile 
description 
(2) 
closed end 
steel pipe 


1.20 
closed end 
steel pipe 


1.50 
closed end 
steel pipe 


1.70 
closed end 
steel pipe 


precast sq. 
concrete 

(16 in. in.) 

pile 

14BP73 

equiv. dia. 

1.34 

Precast sq. 
concrete 

(18 in. in.) 


1.70 
closed end 
Steel pipe 
1.75 


1.42 


1.52 
closed end 
Steel pipe 
diameter 
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Pile 
penetration, 
feet 


average 
52.8 


average 

average 


average 


40.2 


average 


TABLE 3.—Test Pile 


Soil 
description 


(4) 


fine sand, 
fine sand with silt and clay, 
and 


fine sand, 
fine sand with silt and clay, 
and 


fine sand, 
fine sand with silt and clay, 
and 


fine sand, 
fine sand with silt and clay, 
and 


fine sand with silt clay 
laminations, and 


fine sand, 
fine sand with silt and clay, 
and 


fine and medium silty sand 
with lenses silt and 
clay, SP-SM 


silt 
fine and medium sand 


silt 

fine and medium sand with 
pocket gravel 

silt 

fine and medium sand 

silty sand ft; below 


fine and medium sand 
SW-SP 


and Soil Data 


Sand 
friction 
angle 


degrees 


Sand 
relative 
density 

percentage 
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Effective 
overburden 
pressure 
tons per 
square foot 


Total 


2. 


Failure Loads, 
tons 


Point 


(9) 


References 
(11) 
12, and 


and 


and 


and 


973 
(16) 
1.265 172 85° 87° 
0.774 
0.766 
0.774 
0.765 
0.593 
200 104° 96° 
1.184 
1.532 175° 68° 
0.757 
1.25* 
100 1.615 166° 
100 1.980 
1.25* 
100 1.632 196° 
100 2.013 
1.25* 
100 1.610 151° 180°" 
100 1.969 
347 214 133 
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TABLE 3.— 


H-15 


Gregersen precast cy- 


sandy fill until 5.6 ft; from 


lindrical 5.6 9.8 ft: fine sand, 
Gregersen concrete SP; below 9.8 ft: medium 
0.92 coarse sand, SP. 


“Unadjusted compression test. 
Adjusted compression test. 
test. 
sand layer. 
“Driven vibratory hammer. 

Pile diameter diameter the circle same pile area. 


load, was necessary extrapolate and develop load transfer curve corre- 
sponding the failure load. 

Instrumented tension tests performed part the Arkansas River Navigation 
Project (12) have shown, consistently, that resistive concentrated load exists 
the pile point. Also, when measurements load distribution along the pile 
were made after the compression loads test pile were released, was 
shown that the stresses along the pile were not zero. These observations confirm 
the existence residual stresses hammer driven piles, which can mask the 
real distribution the total load between side and point. Hunter and Davisson 
(12) concluded that measurements these load components (side and point) 
are seriously residual loads are not considered. However, 
the total loads both compression and tension are not changed. 

Even though the presence residual stresses hammer driven piles widely 
recognized (9, 10, 11, 12, 26, 29, and 38), method correcting for residual 
stresses has not been generally adopted. rational method for this purpose 
was presented early 1964 (7). Hunter and Davisson (12) advocated the 
use the same method 1969. Also, (38) 1978 suggested the use 
this method correcting for residual stresses. order use the method, 
necessary conduct tension test after the compression test conducted. 
detailed explanation the procedure used applying the method presented 
Hunter and Davisson (12). 

important consequence the adjustment loads for residual stresses, 
suggested Holloway, Clough, and that the adjusted compressive 
side load should equal the tensile load obtained pullout test. Robinsky 
and others (24,25) pointed out that the residual stresses caused the soil 
the pile driving induce instability the arching system which surrounds the 
pile wall. Since tension (pullout) test reverses the direction applied stresses 
the soil, very likely that the most unstable part the delicate arching 
system destroyed and more permanent side resistance attained. 
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Continued 


Test 


The purpose this study develop improved correlations between bearing 
capacity factors and pile-soil system parameters. order develop correlations, 
field load test data must obtained which include the required pile and soil 
data. Evaluation the bearing capacity factors, and used 
and requires knowledge the unit resistances, coefficient friction, and 
effective soil pressure. addition, for correlation purposes, necessary 
know the friction angle relative density the sand, both. 

For this study the evaluation the sand parameters was based Standard 
Penetration Test data conjunction with soil descriptions. Similar procedures 
were used obtain the unit weight the sand which needed for evaluation 
the effective overburden pressure. Also, information regarding the location 
the ground water had known. 

The average unit resistances, and can calculated explained 
the previous section from instrumented pile test data. Another way obtaining 
the side and point loads and corresponding resistances assume the tension 
test load equal the adjusted compressive side load proposed Holloway, 
al. (11). The point load obtained subtracting the tension test load from 
the total compression test load. For purposes data analysis this study, 
loads and unit resistances obtained this manner are referred ‘‘compres- 
data. 

Field load test data from tests were found the literature, but only 
them included both compression and tension test data. preliminary analysis 
all test data indicated that the best correlations (less data scatter) were achieved 
using data. Therefore, although both unadjusted compres- 
sion and compression/tension test data were analyzed, only the compres- 
test results are presented this paper. 

Table presents fairly detailed summary the test pile and soil data. 


1.822 22° 27° 
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this table pile diameter, the case piles with noncircular cross 
sections, diameter circle with the same cross-sectional area. should 
noted that the ground-water level all test sites was very close the 
ground surface. Table presents summary the unit resistances and bearing 
capacity factors which were determined for these tests. The ratio penetration 


TABLE 4.—Unit Resistances and Bearing Capacity Factors 


UNADJUSTED COMPRESSION TEST COMPRESSION/TENSION TEST 


1.05 
1.22" 
0.510 


127.0 


Pile Test 


JONESVILLE 


VESIC 


GREGERSEN 


"ha justed compression test data. 
Note: tef + 95.76 Ove. 


depth pile diameter (relative depth) also tabulated since was used 
the correlations presented this paper. 

When dealing with soil investigation considerable engineering judgment has 
used extrapolating soil properties from bore holes the test location. 


0.76 
1.45 
26.0 
6.49 
4 
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This becomes particularly difficult for soils varying composition, and sometimes 
the tested pile acts probe uncovering the influence minor occurrences 
the soil strata. That was the case for the Jonesville piles, where longer piles 
had less load capacity. There were indications that ‘‘clay greater 
depths created planes weakness which could explain the reduction capacity. 


TABLE 5.—Comparison Field K-Values 


0.7 
Unadjusted Adjusted Tension unadjusted 
(2) (3) (4) (5) 
Arkansas 
Arkansas 
Arkansas 
Arkansas #10 


Arkansas 


Low-Sill 
Vesi¢ 


° 
B Jonesville 
© Gregersen 


Bearing Capacity Factor, Ng 


FIG. 4.—Bearing Factor Versus 


For this reason only the shortest pile, Jonesville #1, was used for correlation 
purposes. 

The Low-Sill piles were tested site overlaid silt stratum about 
(15 thick. All measurements given refer the underlying sand stratum, 
with silt influencing only the effective overburden pressure and depth 
penetration. Meyerhof (21) suggests that, layered soils, penetration about 
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pile diameters one stratum enough mask the influence the thickness 
and strength the upper layer. 

The Gregersen piles represent tests made the same pile different 
penetrations. The driving residual stresses were analyzed and the measured 
loads adjusted (9). However, the method adjustment was different from 
that presented Hunter and Davisson (12). This last method would yield larger 
adjusted point load. 

Values from the Arkansas tests are summarized Table for the 
different loading conditions. Comparison these values shows clearly that 


Unit Point Resistance, q,, (tsf) 


Example Notation : 
(60) ~ Relative Density, 
(@) 34 + Friction Angle, ¢' 


— 
= 
> 


Arkansas 


Jonesville 
Low - Sill 


FIG. 5.—Ultimate Unit Point Resistance Versus 
95.8 


the values adjusted for residual stresses are very close the actual tension 
test values. fourth set values given Table These values were 
determined multiplying the unadjusted values 0.7. These values are also 
agreement with the tension test values. This comparison explains the commonly 
accepted practice (18) assuming that tension load 70% the unadjusted 
compression side load. 

Comparison values obtained from field test data and theories given 
Fig. Curves are plotted for the better known theories, covering wide 
range possible values. Also, the values obtained from field tests 


a 20 40 60 80 100 120 140 150 
20 —+ 3 36° 37° 
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are plotted for the data. Values were also determined 
for the unadjusted compression tests but there was considerably more data 
scatter. The results Fig. show that none the theories can used 
predict value. Even though the Terzaghi general shear theory 
seems favored, the measured values deviate widely from the theoretical 
values. This type comparison clearly shows the need for improved correlations 
which include all significant parameters involved the pile-soil interaction 
problem. 


According the static method determining pile bearing capacity the unit 
resistances, and are functions mutually independent factors: N,, 


Bearing Capacity Factor, Ng 


40 60 8 100 
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FIG. 6.—Bearing Factor Versus 


tan and effective pressure. Consequently, the bearing capacity factors can 
evaluated separately from the overburden pressure depth. However, 
experimental investigations have indicated that the bearing capacity factors 
and are dependent upon overburden pressure depth. result, there 
not strong justification for the individual determination Even 
so, order enhance the understanding pile-soil interaction, these factors 
were also included the correlations this study. 
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| | 
| ¥ | 
| F 
40) 
| | 


980 


JULY 1981 GT7 


The field load test data were used the development new correlations 
which relate the bearing capacity factors pile geometry parameters and sand 
properties. Several different combinations pile parameters and sand properties 
were investigated during the development these correlations. However, the 
correlations presented are considered the simplest and best for practical usage. 
The correlations were developed for the compression/tension test method 
interpreting pile load tests. 


Unit Side Resistance, f, (tsf) 


Example Notation : 
(50) + Relative Density, Dp 


Arkansas 


Jonesville 
Low 

Vesic 

Gregersen 


Relative Depth, 0/8 


10 


36 


FIG. 7.—Ultimate Unit Side Resistance Versus 
tsf 95.8 


Fig. presents the plot unit point resistance, g,, versus relative depth. 
Relative depth determined using the depth the pile point. The influence 
friction angle (or relative density) clearly defined that higher fraction 
angles (or relative densities) indicate higher unit point resistances. this case 
the influence the friction angle, compared the relative density, dominant. 
With minimal scatter, curves are developed for 30° (very loose sands), 
and 35° (dense sands). The other curves were developed interpolation 
and extrapolation. pile diameters penetration, has not reached 
constant value. reduced rate increase indicated below pile diameters 
penetration for loose sands, and below approximately pile diameters for 
denser sands. Below these critical depths the relationship between and relative 
depth nearly linear. 
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Fig. shows plot the bearing capacity factor (plotted log scale) 
versus relative depth. general, the values increase with increasing friction 
angle (and increasing relative density). The data permitted the development 
curves for 30° (very loose sand), 34°, and 36° (medium 
dense sand). The other curves were developed using extrapolation and 
interpolation. These curves are extended zero penetration (shallow foundations) 
accordance with the shapes similar curves determined experimentally (21,37), 
and the values for zero penetration were obtained using Terzaghi’s general 
shear theory. can seen that the values increase from zero penetration 


Combined Factor, K tané 
6 


General Equation: | 
=Log( ktané)- 


= -0: 


Example Notation : 
(50) + Relative Density, Op 
oO 33 + Friction Angle, ¢ 
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FIG. 8.—Combined Factor tan Versus 


maximum value roughly pile diameters. deeper penetrations the 
values seem decrease linearly with depth. 

Fig. shows the plot the unit side resistance, versus relative depth. 
Relative depth determined using the midpoint depth pile segment being 
analyzed. The effect increasing side resistance with increasing friction angle 
(or increasing relative density) indicated. Three well defined curves are 
developed: one curve for 30° (very loose sand); another for 33° 
(medium dense sand); and third for 37° (dense very dense sands). 
Below approximately pile diameters penetration for the loose sands, 
practically linear relationship and relative depth indicated. similar 
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relationship suggested for denser sands, deeper penetrations. 

Fig. shows the plot the combined factor tan (plotted log scale) 
versus the relative depth. Fig. presents the lateral earth pressure coefficient 
(plotted log scale) versus relative depth. The coefficient was determined 
using 0.8 times the peak friction angle. Three basic curves based friction 
angie (or relative density) are identified, and additional curves are obtained 
interpolation. For these semilog plots the relationships are linear and curves 
with constant slopes are indicated for each factor correlated. The equations 
for each basic curve are presented for each factor. The values correspond 


Lateral Earth Pressure Coefficient, K (5 = .8¢') 


40 59 


General Equation: 


Example Notation : 
(50) + Relative Density, Op 
© 33 + Friction Angle, ¢' 
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FIG. 9.—Coefficient 0.8 Versus B—Compression/Tension 


partially developed passive earth pressure coefficient shallower depths, 
and decrease order magnitude active coefficients the relative 
depth increases. 


result this study has been shown that the parameters which provide 
the best design correlations for piles sand are the relative depth (depth 
diameter ratio) and the sand friction angle. The relative depth probably related 
the effective overburden pressure the soil surrounding the pile, and the 
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volume displaced and deformed soil. All soils investigated this study have 
the same quartzose composition and most have fine gradation. Consequently, 
the sand friction angle, addition reflecting the shear strength, may 
related relative density (compressibility). 

For the unit point resistance, g,, constant values were not indicated 
penetrations diameters. And for the unit side resistance, constant values 
were not indicated penetrations diameters. The unit resistances showed 
pronounced increase with increase relative depth for shallow penetration. 
However, the rate increase average unit resistance with relative depth 
becomes smaller, and relatively constant, greater penetration (roughly between 
pile diameters and pile diameters). 
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FIG. Total Axial Load Error (from Versus B—Compres- 


The use the correlations for and recommended for estimating 
pile capacity. Fig. shows the results error analysis based predicted 
ultimate loads for pile load tests sand. Eight these piles were not 
used developing the design correlations. This analysis shows that all the 
loads were predicted within error band +30%, and the majority the 
predicted loads were within error band +20%. Load predictions were 
also made using the correlations for tan and but the errors were 
slightly larger. 
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The following symbols are used this paper: 


area pile point; 

area pile shaft; 

pile diameter (or, for noncircular piles, the diameter 
circle with the same area the pile cross section); 

cohesion soil; 

depth pile penetration; 

relative density sand; 

relative depth; 

ultimate unit side resistance; 
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lateral earth-pressure coefficient; 

back calculated value from unadjusted compression 
test, multiplied 0.7; 

combined bearing capacity factor for circular and square piles 
(includes shape factor); 

bearing capacity factors terms vertical overburden 
pressure; 

effective overburden pressure; 

average vertical effective overburden pressure along pile 
segment; 

effective vertical overburden pressure the pile point level; 
ultimate load carried the pile point; 

ultimate load carried pile shaft; 

total axial load pile; 

ultimate total axial load pile; 

ultimate unit point resistance; 

shape factors; 

effective unit weight soil; 

angle friction between pile and soil; 

peak friction angle soil based effective stresses; and 
residual friction angle soil, based effective stresses. 
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Note.—Discussion open until December 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. This paper part the Journal the Geotechnical Engineering Division, 
Proceedings the American Society Civil Engineers, Vol. 107, No. GT7, 
July, 1981. 
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Notes 


provide place within ASCE for publication technical ideas that have not advanced, 
yet, the point where they warrant publication Proceedings paper Journal, the 
publication Technical Notes was authorized the Board October 16-18, 
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RESPONSE CIRCULAR FOOTINGS 
VERTICAL VIBRATIONS 


Mysore and Asuri Sridharan” 


INTRODUCTION 


Sung (5) used half-space model find the displacement footing resting 
elastic half-space for various interface pressure distributions, viz. rigid 
base, uniform and parabolic. From the available theoretical and experimental 
investigations (1,3,4) can concluded that the variation contact pressure 
distribution should considered while solving problems dynamics 
foundations. this paper, analog solutions (2) are presented for the response 
circular footing resting elastic half-space with uniform and parabolic 
contact pressure distributions and subjected frequency dependent and fre- 
quency independent excitations. addition, analog solution rigid circular 
footing subjected frequency dependent excitation also presented. The results 
have been compared with the rigorous solution Sung (5) and the agreement 


Uniform Pressure spring constant, due uniformly 
distributed load given (5) 


which shear modulus; Poisson’s ratio; and radius the 
footing. The displacement functions, and for uniform pressure distributions 
given Sung (5) for four Poisson’s ratios, when multiplied the factor 
become almost independent The average value and Kf,/Ga 
are represented and F,, respectively, and are called displacement 
From these values and F,, the factor can expressed 
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which nondimensional frequency factor. suggested Lysmer and 
Richart (2), the value damping coefficient, should taken the average 
the values the range 0.3 0.8. The average value 
this range found equal 0.674. Using this value the analog for 
damping coefficient can expressed 


which mass density soil. Thus, the analog model for footing resting 


40 


FIG. 1.—Comparison Analog Solution with Elastic Half-Space Theory for Response 
Footing-Soil System with Uniform Pressure Distribution-Frequency Dependent 
Excitation 


elastic half-space with uniform pressure distribution the interface can 
written 


m 


function. Using the analog parameters for stiffness (Eq. and damping (Eq. 
3), the following parameters can written: 
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which and magnification factors for frequency dependent and 
independent excitations, respectively. 

Fig. compares the results the analog solution (Eq. with the rigorous 
solution Sung (5). The agreement found good. Similar good agreement 
was found for frequency independent excitation case also (Eq. 7). The magnifica- 
tion factor resonance and for frequency dependent and frequency 
independent excitations respectively are expressed (from and 7). 


UNIFORM PAR BOLIC 


FIG. 2.—Variation Magnification Factor Resonance with Modified Mass Ratio 
for Rigid, Uniform and Parabolic Pressure Distributions 


992 JULY 1981 


Substituting for (Eq. Eq. 


m 


The variation M,,) with shown Fig. Fig. also presents 
the variation M,,, and with obtained from the rigorous solution 
Sung (5) for three types pressure distributions. seen that there 


m~ m 
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FIG. 3.—Variation Frequency Factor Resonance with Modified Mass Ratio 
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good agreement between the results the analog solution and the rigorous 
solution. Fig. shows the variation frequency factor maximum magnification 
factor, with for both analog and rigorous solutions bringing out good 
agreement between them. 

Parabolic Pressure same procedure before has been 
followed for parabolic pressure distribution also. The spring constant, 
given (5) 


The expression for F,) can written 


The average value taken the range 0.3 0.8 found 


— EXACT SOLUTION 
ANALOG SOLUTION 


Analog Solution with Elastic Half-Space Theory for the 
Response Footing-Soil System with Parabolic Pressure Distribution-Frequency 
Dependent Excitation 
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The analog model for footing resting elastic half-space with parabolic 
pressure distribution the interface can thus written 


Using the analog parameters (Eq. 10) and (Eq. 12) the following parameters 
can written 


M, = 2\2 2431/2 


Fig. compares the results the analog model (Eq. 15), with the exact 
solutions Sung (5). The agreement found good. Similar good agreement 


inwhich 


and 


was found for frequency independent excitation (Eq. 16). The magnification 
factor resonance, and M,, can obtained from and 


cm m 


The variation M,,,(= with shown Fig. comparing the results 
the analog model with the results the rigorous solution (5) presented 
Fig. found good agreement. can also seen from Fig. 
that there good comparison the variation with obtained 
from the analog and rigorous solutions. 


Base Pressure 


Analog solution for frequency dependent excitation for rigid base pressure 
distribution has not been reported the literature. Using the values and 
given Lysmer and Richart (2) the equation for the magnification factor 
for rigid circular footing subjected frequency dependent excitation can 
written 


r 


The magnification factor resonance remains same for both frequency 
dependent and frequency independent excitation and can expressed 
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The variation with shown Fig. Fig. presents the variation 
with 


Analog models which have distinct advantages have been developed (Eqs. 
and 13) predict the response circular footing resting half-space, 
assuming uniform and parabolic pressure distributions below the footing, subject- 
frequency dependent and frequency independent exciting forces. The 
response rigid circular footing subjected frequency dependent excitation 
also found. The comparison the results from the proposed analog model 
with the results elastic half-space theory has shown good agreement. 
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RANDOM GENERATION 
STABILITY ANALYSIS 


and Charles ASCE 


INTRODUCTION 


Computerized solution the factor safety for slope stability via limiting 
equilibrium analysis and the method slices extremely common. However, 
computerized failure surface generation and analysis has been limited simple 
regular geometric shapes, viz., the circle. more complex surfaces need 
examined, they are manually generated prior computerized analysis. 
Computer program STABL (1,5), used routinely the Indiana State Highway 
Commission for yr, has unique searching routines which permit any failure 
surface shape and position generated and analyzed the computer. Dozens 
engineering firms and agencies over the world have found the totally 
computerized searching routines STABL valuable, and the object 
this paper acquaint others with them. 

STABL uses the simplified Bishop assumption resolve the familiar statical 
indeterminancy, and gives essentially the same factor safety for selected 
specific failure surface any other readily available routine (1,4). Users have 
applied STABL concurrently with their older approaches the stability problem 
order confident the factors safety generated, and detect any 
tendency for STABL slightly more less conservative. The uniqueness 
STABL lies its computerized searching routines, which allow hundreds 


positions and shapes failure surfaces both generated and analyzed 
for few dollars. 


TECHNIQUES 


The number slopes with subsurface profiles for which one can intuitively 
select the shapes and positions the more critical shear surfaces are few. 
some cases the positions may more obvious than the shape, e.g., along 
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layer seam low strength. Even then, may neither obvious where 
the critical shear surface enters and exits the layer, nor what path might 
take beyond these points. necessary try numerous shear surface shapes 
before the critical one identified. other cases, the shapes may more 
obvious than the positions, e.g., the circle for slopes roughly homogeneous 
and isotropic materials. The familiar computerized circular grid search evidence 
the usual need iterate for many potential circle positions. still other 
cases the soil profile complex that both the shape and the position are 
uncertain. these cases the options are to: (1) Manually input specific surfaces 
and analyze them manually computer; (2) use the appropriate STABL 
searching routines where surfaces are both generated and analyzed the 
computer. 

The STABL searching routines greatest interest are BLOCK and RANDOM, 
since these deal with the generation irregular shear surfaces. The CIRCLE 
searching routine within STABL (2,3,6) differs little from other readily available 
routines and will not described here. 


Use 


The key searching routines for irregular surfaces random generation 
these surfaces. This requires both random number subroutine and experience 
the application the random values. Since most computer systems have 
library function for generating random numbers, the first requirement easily 
met. Control Data Corporation has such function, named RANF(X), its 


6000 series system library, and this function used STABL. 

substantial amount research and experience was required bring the 
random surface generation concept practical use (1,5). Randomly generated 
shearing surfaces can course nonsensical both shape and position, 
which means that much the computer output wasted. avoid this, all 
surfaces are generated series straight lines which: (1) Initiate from 
logical range positions the boundary; (2) form shearing surface 
kinematically acceptable shape; (3) lie above lower boundary formed 
strong layer; and (4) terminate within logical range positions the boundary. 

Plotting routines are essential the exercise the necessary logic, and 
are used the following way STABL. Fifty 200 shearing surfaces are 
batch-generated and analyzed based upon set shape and position assumptions. 
plot this many surfaces has the appearance large ink blob, although 
does show the part the subsurface which has been searched. more 
instructive plot one that shows only the 5%-10% the surfaces with the 
lower factors safety. These surfaces form what termed zone”’ 
(5). 

The compactness the critical zone, its location within the zone containing 
all the surfaces generated, and the magnitude the range values for the 
factors safety the more critical surfaces, indicate the likelihood that 
shear surface exists with factor safety significantly lower than any calculated. 
the critical zone narrow, noncritical surfaces have been generated 
both sides the critical zone, and the values the more critical surfaces 
are nearly the same, the defined factor safety the most critical surface 
generated could accepted, with reasonable confidence. the other hand, 
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the critical zone wide, the more critical surfaces lie along one edge 
the zone containing all the surfaces generated, the magnitude the range 
factor safety values for the critical surfaces large, all these, 
the possibility shear surface with significantly smaller value for the factor 
safety may require generation additional surfaces. more trial shear 
surfaces are required, the limitations can revised restrict the additional 
surfaces particular zone interest. 

The necessity several man-machine communications for analyses complex 
slopes common all computer solutions. The searching and plotting routines 
STABL accommodate this need rather easily and efficiently. 


begin, line segment projected from initiation point direction 
chosen between two limits. The direction limits are either specified automati- 
cally calculated the program and are shown Fig. the latter case, 
the clockwise direction limit, a,, —45° with respect the horizontal, and 


Counterclockwise Limit 
Point 


— 


Limit 


FIG. Initial Line Segment 


the counterclockwise direction limit, a,, less than the inclination the 
ground surface the right the initiation point. 

The selection the inclination, which the initial line segment 
surface projects, determined randomly. random number, having value 
between and generated and squared. The squared random number and 
the difference the inclination the two direction limits, are multiplied. 
The result added the inclination the clockwise direction limit, a,, 
obtain the angle inclination the line segment: 


Squaring the value the random number generated, before multiplying 
with the absolute value the difference inclination the two direction 
limits, introduces bias the random direction selection process. 
Although selection any inclination the initial line segment within the two 
limits possible, more likely that the initial line segment located nearer 


B 
-45° 
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the clockwise limit. The primary reason for introducing the bias obtain 
good distribution completed surfaces. When the initial line segment was 
given equal likelihood being oriented any direction within the direction 
limits, poor distributions completed surfaces were obtained. 

The direction limits for successive line segments are established shown 
Fig. The counterclockwise limit, normally deflected 45° counter- 
clockwise from the projection the preceding line segment. for particular 
line segment the orientation the counterclockwise direction limit greater 
than 90° (measured with respect the horizontal), the inclination the 


Counterclockwise Limit 
45° 


Projection of Previous 
Line Segment 
2 


Clockwise Limit 
FIG. 2.—Direction Limits for Successive Line Segments Irregular Surface 


—-— Beyond Right Termination Limit 
——--— Below Depth Limitation 


FIG. 3.—Trial Failure Surface Acceptance Criteria 


counterclockwise direction limit, adjusted 90°. 
The clockwise deflection limit for the direction line segment, 
randomly selected for each shear surface generated, 


for particular line segment the inclination the clockwise direction limit, 
less than —45°, set this value. 


8, 
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The inclination line segment then established by, 


which the angular extent permissible directions, multiplied 
random number raised random power between and and added 
the inclination the clockwise direction limit, 

The foregoing procedure, although somewhat arbitrary, does produce irregular 
surfaces reasonable shape and position (5). Reverse curvature shear 
surface possible, but the frequency occurrence small. Unless very short 
line segments are used, ‘‘kinkiness’’ the resulting generated shear surfaces 
not problem. 

shown Fig. surface should intersect the ground surface short 
the left termination limit, rejected, and replacement surface generated. 
surface should intersect the ground surface beyond the right termination 
limit, uncompleted surface extends beyond this limit, the possibility 
successfully completing the surface without rejecting the entire surface 
checked. First, the violating line segment rotated the inclination the 
counterclockwise direction limit, 6,, which was used select it. this 
inclination the surface intersects the ground surface beyond the right termination 
limit, extends uncompleted beyond this limit, the surface rejected, and 
replacement generated. 

line segment extends below the maximum depth permitted, rotated 
its counterclockwise direction limit. still extends below the depth limit, 
the entire surface rejected, and replacement generated. not, 
adjusted the old inclination the line segment, and new inclination 
selected before. The procedure repeated necessary. 


zone has been fairly well defined previous usage the circular irregular 
surface generators, often advantageous use third generator. 

series boxes, minimum two, located along the path the zone 
investigated. The boxes are parallelograms with vertical sides. They can 
specified bracket the zone investigated shown Fig. point 
randomly selected from within each box. Any point within particular box 
has equal likelihood being selected. The points are then connected 
sequence with straight line segments (Fig. 4). 

Next, the active and passive portions the trial shear surface are generated 
shown Fig. The active and passive portions are composed line segments 
equal specified length. The inclination each line segment for the active 
portion the shear surface selected biased random fashion between 
direction limits inclined 45° and 90° 


For the passive portion the shear surface, the inclination each line segment 


selected biased random fashion between direction limits inclined —45° 
and 0°: 
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The same judgment, considered for critical shear surface searching with 
the irregular-shaped shear surface generator, can also applied this generator, 


intensive Search of Critical Zone Previously 
Defined by CIRCLE or RANDOM 


Search in Irreguior Weok Layer 


FIG. 4.—Sliding Block Generator Using more than Two Boxes 


Active 45° Direction 
Range 


Passive Surface 


Base 
FIG. 5.—Generation Active and Passive Portions Sliding Surface 


with regard the probable existence significantly more critical shear surface 
than those previously generated. 


Summary 


Random methods provide convenient means for generation trial shear 
surfaces for computerized slope stability analysis. Hundreds trial surfaces 
can examined reasonable cost. Three trial shear surface generators 
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have been developed STABL, and two have been described this paper: 
(1) Surfaces general irregular shape and (2) sliding-block-shaped 
surfaces with random orientation and location for the bases the central block 
and active and passive type wedges (BLOCK). 

All three generators have been routinely employed for approx the 
Indiana State Highway Commission, well dozens other users both 
the United States and abroad. Accordingly, they have been rather thoroughly 
debugged. all cases users have been advised run STABL 
with their older established methods, order gain confidence the newer 
searching modes observe any consistent tendency for STABL produce 
higher lower factor safety values, both. The details the generators 
are presented for the purpose encouraging engineers either: (1) Incorporate 
them their present systems; (2) expand and improve upon them more 
advanced systems. 
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TECHNICAL NOTES 


DILATION ANGLE AND LIQUEFACTION POTENTIAL 


Most our understanding the liquefaction phenomenon has come from 
laboratory cyclic loading tests. These tests indicate that the relative density 
the sand important parameter controlling its resistance liquefaction. 
the field difficult measure relative density and common practice 
infer the relative density and liquefaction resistance sand from penetration 
tests. The most commonly used relationship between penetration resistance and 
relative density the one proposed Gibbs and Holtz (5). More recently 
Seed, al. (9) have proposed that the penetration resistance corrected 
Penetration resistance tests cannot performed small laboratory samples 
and thus cannot correlated directly with laboratory cyclic load test behavior. 
would desirable have single parameter describing the initial state 
the sand which could obtained both the laboratory and the field and 
which would give measure liquefaction resistance. proposed herein 
that the dilation angle expansion rate the sand during shear such 
measurement. can readily obtained the laboratory from drained triaxial 
simple shear tests and the field can obtained from self-boring 
pressuremeter tests described Hughes, al. (7). 

liquefaction resistance chart similar that proposed Seed (8) but based 
dilation angle rather than corrected blow count presented. From this chart 
estimate the field liquefaction resistance saturated sands can obtained 
from dilation angles measured from self boring pressuremeter tests. 


Fig. shows stress-strain behavior Ottawa sand under drained simple 
shear conditions using constant vertical confining Over considerable 
range strain, both initially loose and dense samples undergo volume expansion 
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(dilation) and very large shear strains tend approach ultimate void 
ratio. However, the rate volume expansion with shear strain larger for 
the dense than for the loose sand. The dilation rate sand characterized 


the dilation angle (6), which the inverse sine the slope volume 
expansion curves Fig. 


simple shear tests were carried out using higher value confining pressure, 
results similar Fig. will obtained; but, given initial void ratio (or 
relative density D,), increase will result decrease and also 


Volumetric 


Shear (%) 


FIG. Change Behavior Ottawa Sand Drained Simple Shear 


the ultimate void ratio. The dilation angle thus function both relative 
density and confining pressure. 

clear from the results shown Fig. that the dilation angle 
sand closely related its relative density. The value was computed 
shear strain 10%, and shown plotted against the corresponding 
relative density Fig. linear relationship was obtained for the sand tested. 
This relationship, however, applies for vertical confining stress 200 kPa. 
pointed out earlier, the dilation angle also depends the level confining 
stress. The variation dilation angle with relative density and vertical confining 
pressure for Leighton Buzzard sand was studied analyzing data from 
comprehensive series simple shear tests Cole (3). was found that for 
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this sand also linear relationship exists between and each constant 
confining pressure. Furthermore, each relative density, decreases linearly 
and constant rate with confining pressure. The dilation angle confining 
pressure ton/sq (100 kPa) about 2-1/2° higher than ton/sq 
(200 kPa). assumed that Ottawa sand behaves manner similar 
Leighton Buzzard sand with respect the dependence confining pressure, 
linear relationship between dilation angle corrected confining pressure 


ton/sq (100 kPa) and relative density will obtained. This shown 
dotted Fig. 


Anatysis 


The liquefaction resistance saturated sand obtained from laboratory 
tests has been expressed functior relative density many researchers 
including Seed (8) and Finn and the first writer (4). performing drained 
Static simple shear tests the same sand over range relative densities 
and observing the volume change characteristics, its dilation angle function 
relative density can obtained. This allows the liquefaction resistance 
expressed terms dilation angle and thus obtaining situ measurements 


Ottawa Sand C-109 


a 
< 


Relative (%) 
FIG. between Dilation Angle and Relative Density 


dilation angles from self-boring pressuremeter tests estimate the 
situ liquefaction resistance saturated sand obtained. 


The liquefaction resistance Ottawa sand obtained from constant volume 
cyclic simple shear tests shown Fig. Here the cyclic stress ratio, 
cause liquefaction 10% double amplitude shear strain cycles shown 
function relative density. All tests were carried out using constant 
confining pressure 200 kPa. may seen that the liquefaction resistance 
increases with increasing relative density and very markedly for relative 
densities excess about 70%. 

Since from Fig. the dilation angle corrected normal pressure 
ton/sq (100 kPa) uniquely related the relative density, the liquefaction 
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40 60 80 
Relative Density, D, (%) 


Corrected Dilation (*) 


5 10 8 2025 40 
Corrected Blow Count,N, 


FIG. 3.—Resistance Liquefaction Sand Function Relative Density, Dilation 
Angle Penetration Resistance 


100 


Blow 


FIG. 4.—Relationship between blow count corrected ton/sq and Relative 
Density (after Gibbs and Holtz, 1957) 
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resistance can also related the dilation angle. This shown Fig. 
added scale for dilation angle v,, besides the relative density scale. This 
figureis thought give low conservative estimate the liquefaction resistance 
saturated medium sands. Thus obtaining the dilation angle the field 
described Hughes, al. (7) conservative estimate liquefaction resistance 
obtained. 

The chart shown Fig. was used estimate the liquefaction resistance 
dam British Columbia. This dam was constructed the hydraulic fill 
method and the core comprised nonplastic silt. Dilation angles measured 
from self boring pressuremeter tests and corrected effective confining stress 
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resistance corrected 
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@ Liquefaction — stress ratio bosed on estimated acceleration 
@ Liquefaction — stress ratio based on good acceleration data 
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© No liquefoction— stress ratio based on good acceleration data 


FIG. 5.—Correlation between Stress Ratio causing Liquefaction Field and Penetra- 
tion Resistance Sand (after Seed, al., 1975) 


ton/sq were the range From Fig. this material could 
expected have dynamic resistance ratio about 0.08. Laboratory cyclic 
load tests undisturbed samples obtained from locations close the site 
for pressuremeter tests, indicate that the dynamic resistance ratio 0.10, and 
thus the chart gives reasonable but conservative estimate liquefaction 
resistance for this case. 


The most commonly used relationship between relative density, blow count 
and confining pressure that proposed Gibbs and Holtz (5). For confining 
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pressure ton/sq ft, the relationship shown Fig. Using this 
relationship combined with the liquefaction resistance versus relative density 
curve Fig. the liquefaction resistance terms blow count normalized 
ton/sq obtained. This shown the solid line Fig. Field 
experience presented Seed (8) together with his liquefaction resistance curve 
(dashed line) also shown this figure and are seen very close 
agreement with the solid line for values less than about 20. 

Liquefaction resistance can correlated with relative density, corrected 
dilation angle corrected blow count and chart showing liquefaction resistance 
terms these three parameters shown Fig. indicates that liquefaction 
not likely occur regardless the stress level provided the relative density 
excess 75%, the corrected dilation angle greater than 16°, the 
corrected blow count exceeds 20. Such conclusion seems apparent from analyses 


field records liquefaction Seed (Fig. 5), Christian and Swiger (2) and 
Castro (1). 


ACKNOWLEDGMENTS 


The financial assistance Natural Science and Engineering Research Council 
Canada gratefully acknowledged. 


Castro, G., and Cyclic Mobility Saturated Journal the 
Geotechnical Engineering Division, ASCE, Vol. 101, No. GT6, Proc. Paper 11388, 
June, 1975, pp. 551-569. 

Christian, T., and Swiger, F., ‘‘Statistics Liquefaction and SPT 
Journal the Geotechnical Engineering Division, ASCE, Vol. 101, No. Proc. 
Paper 11701, Nov., 1975, pp. 1135-1150. 

Cole, L., Behavior Soils the Simple Shear Ph.D thesis 
presented Cambridge University, England, 1967, partial fulfillment the 
requirements for the degree Doctor Philosophy. 

Finn, L., and Vaid, P., Volume Cyclic Simple Shear 
Proceedings 2nd International Conference Microzonation for Safer Construction, 
Vol. 1978, pp. 839-851. 

Gibbs, J., and Holtz, G., Determining the Density Sands 
Spoon Penetration Testing,’’ Proceedings, 4th International Conference Soil 
Mechanics and Foundation Engineering, Vol. 1957, pp. 35-39. 

Hansen, B., Ruptures Regarded Narrow Rupture Zones; Basic Equation 
Based Kinematic Proceedings, Brussels Conference Earth 
Pressure Problems, 1958, Vol. pp. 39-48. 

Hughes, O., Wroth, P., and Windle, D., Tests 
Geotechnique, Vol. 27, No. 1977, pp. 455-477. 

Seed, B., Soil Liquefaction Effects Level Ground During 
presented the September 1976, ASCE Annual Convention 
and Exposition, held Philadelphia, Pa. (Preprint 2752). 

Seed, B., Mori, K., and Chan, K., Seismic History the 
Liquefaction Characteristics Sands,’’ Report No. EERC 75-25, Earthquake Engineer- 

ing Research Center, University California, Berkeley, Calif., 1975. 


JULY 1981 


DISCUSSION 


Note.—This paper part the Journal the Geotechnical Division, Proceedings 
the American Society Civil Engineers, Vol. 107, No. GT7, July, 1981. 


1009 


16356 GT7 


Discussions 


Discussions may submitted any Proceedings paper technical note published any 
Journal any paper presented any Specialty Conference other meeting, the Proceedings 
which have been published ASCE. Discussion paper/technical note open 
anyone who has significant comments questions regarding the content the 
note. Discussions are accepted for period months following the date publication 
paper/technical note and they should sent the Manager Technical and Professional 
Publications, ASCE, 345 East 47th Street, New York, N.Y. 10017. The discussion period may 
extended written request from discusser. 

The original and three copies the Discussion should submitted (220-mm) 
(280-mm) white bond paper, typed double-spaced with wide margins. The length 
Discussion restricted two Journal pages (about four typewritten double-spaced pages 
manuscript including figures and tables); the editors will delete matter extraneous the 
subject under discussion. Discussion over two pages long will returned for shortening. 
All Discussions will reviewed the editors and the Division’s Council’s Publications 
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manuscripts. 
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Discussion Claude Fetzer,’ ASCE 


The research described the subject paper should keen interest 
all engineers engaged the design, construction, and inspection dams. The 
research directed toward narrow fractures rock. Most published seepage 
and piping problems related foundation rock, including Wolf Creek and Teton 
dams, have been associated with relatively wide fractures and wide solution 
joints. Nevertheless, there need for research develop criteria for placement 
compacted soil, filter, and transition materials against fractured rock. 

The placement compacted soil against fractures rock somewhat 
analogous placement gravel packs around slotted well screens. From years 
experience criteria have been developed for selecting slot sizes compared 
the grain size the gravel pack. Materials ranging from coarse sand 
fine medium gravels are commonly used gravel packs. According Smith 
(11), the criteria select size well screen (slot size) that will retain 
90% more the gravel pack; other words, the width the screen slot 
set the D,, size the gravel pack. Many designers select the slot widths 
for natural formations the ranges the size the formation material. 
both cases recognized that the slot opening must smaller than the 
median size the sand gravel particles. 

Table the authors’ paper, most dimensions the largest grains collected 
from the effluent exceed the width the shim thickness. These data indicate 
that any criteria established for grain size versus width rock fracture should 
probably very conservative basis such untreated rock fracture 
should not wider than the D,, size the adjacent compacted soil, filter, 
transition material. This criteria would require sealing very narrow fractures 
even for sand filters. 

there only one chance treat the foundation rock, considered 
appropriate continue the current practice cleaning out all fractures 
rock depth times the width and refilling with sand cement grout. 
Where the fractures openings are too numerous treat individually, the 
appropriate sections the embankment can isolated from the bedrock opening 
continuous concrete facing was used extensively Patoka Dam (12). 


Smith, J., Selection Gravel Pack Key Successful The 
Johnson Drillers Journal, July-Aug., 1975, pp. 1-3. 
12. Fetzer, A., Dam USCOLD Newsletter, Nov., 1979, pp. 6-9. 


“May, 1980, Richard Goodman and Panchanatham Sundaram (Proc. Paper 
15433). 
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Discussion Reginald Barron,* ASCE 


The authors are commended for initiating studies this area 
engineering. Recent incidents indicate serious need for such information. 

The writer puzzled, however, the technique used the study. would 
appear that the use fixed openings would have been more suitable obtain 
data the piping—grain size characteristics various soils for various joint 
openings. The Waterways Experiment Station (Army Corps Engineers), has 
done some work this, but the goal was obtain design criteria for well 
screens. may that existing filter criteria could used for design purposes, 
but large degree conservatism would required because the wide 
variation magnitude joint openings. 

possible that the greatest danger piping soil the rock face 
not into the rock but seepage the joints direction parallel the 
face. opening exists downstream for the soil piped into (large voids 
the fill, such rock fill, surface rock face) then the flow could develop 
large pipes the earth fill without the fill moving into the rock joints. 

the first sentence the conclusions stated that there seem 
design standards for soils placed against rocks with open joints. This may 
true for grain size criteria, but certainly not for treatment the joints. 
Such treatments have been use for many years and consist of: (1) Cleaning 
the joints the extent practicable; and (2) fill the joints with erosion 
resistance material. For small joints filling slush grout often used; larger 
joints are filled with lean concrete. The fill placed over adjacent the 
jointed rock should fat clay, available. Upstream and downstream 
the impervious zone the rock face should covered with suitable filters. 
example such treatment that used for the Ririe Dam the Army Corps 
Engineers. 

The authors also state that earth was carried from the embankment into karstic 
openings Wolf Creek Dam. The writer not aware this. his knowledge 
pipe did develop the rock under the masonry dam and sink hole developed 
downstream the embankment wrap around and the right stilling basin 
wall. believed that this sink was caused the frequent strong surges 
the stilling basin. There was evidence some the alluvial soil being 
carried into the karstic joints and that, given sufficient time, piping could have 
occurred. was for this reason that the concrete wall was installed the 
embankment and the karstic joints the upper rock foundation. course, 
the problem piping into karstic openings quite different from that piping 
into open rock joints, but the treatment the same; clean out and fill with 
concrete. 


Engr., Horseshoe Road, Guilford, Conn. 06437. 
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Members, ASCE 


The writers appreciate the interest shown this topic Fetzer and Barron. 
Both the discussers pointed out the similarity between placing compacted soil 
against discontinuities rocks, and placing gravel packs around slotted well- 
screens. However, doubtful whether the criteria used for gravel packs 
would sufficient prevent erosion very fine particles from the compacted 
zone into the fractures, since gravel packs customarily exclude fines. Fetzer 
has indicated apparent anomaly between the shim thickness and the maximum 
size particles eroded the experiments. Most the larger particles eroded 
were plate-shaped gypsum fragments whose third dimension, which could not 
measured directly, was presumably smaller than the shim thickness. The 
writers agree with Fetzer that conservatism wise when placing compacted 
soil against rock-face with fracture traces, since sieve analysis does not 
usually give the smallest dimension particle. Regarding the use continuous 
concrete facing, the first writer recommended similar treatment several 
damsites Colombia, where silt was placed against decomposed granite founda- 
tions. 

Regarding the question raised Barron the experimental technique, the 
use shims insured fixed openings for smooth fractures, but rough fractures 
shims could not used. The writers thank Barron for pointing out the different 
types failure mechanisms that could occur the compacted fill without 
the fill material being eroded into the rock fractures. Both types mechanism 
should studied. 

The writers concur with Barron that, wherever possible, fracture openings 
should filled with either grout lean concrete. However, not known 
the writers such treatment was adopted for Teton Dam foundation and 
abutments. Barron was correct asserting that the foundation problem the 
Wolf Creek Dam was not simply result erosion (13). 

The writers wish add that further fundamental studies are necessary 
gain understanding erosion compacted soil into open rock fractures. 


13. Couch, B., and Ressi Cervia, L., cutoff wall installed through 


dam construction Civil Engineering ASCE, Vol. 49, No. Jan., 1979, pp. 
62-66. 


Errata.—The following correction should made the original paper: 


Page 486, line Should read Creek Power Canal Embankment, Washing- 
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The authors’ claim that general constitutive model exists which can 
adequately predict the time-dependent deformation cohesive soil under 
arbitrary three-dimensional systems would disputed this writer 
and various co-workers for whom this has been major research topic for 
more than yr. Various aspects this research have been published 
That both the authors’ and this writer’s programs have had similar origins 
evidenced comparing the authors’ Eq. for linear elasticity with the writer’s 
equations, (e.g., for axisymmetric stress) (33) 


and 


The authors Eq. for pseudo-elasticity, may compared with this writer’s 
procedure (39) establishing time-dependent relationships between the three- 
dimensional stress tensor 


and the principal strain tensor which subdivided follows 


The subdivision Eq. similar, but not identical that adopted the 
authors, which based volumetric and deviatoric strains. The strains 
are defined those which can only develop under drained conditions. The 
Strains can develop under both drained and undrained conditions. More 
importantly, the latter are independent the size specimen tested, whereas 
the former occur rate which determined consolidation conditions. 
This system believed preferable that the authors, when measurements 
made one size specimen have applied specimens different 
size elements field stratum (see, e.g., 44). 

The authors’ confirmatory work, reported the paper all seems have 
been carried out specimens the same size those which were used 
for the evaluation the parameters required. This liable obscure 

1980, Edward Jr. and James Mitchell (Proc. Paper 15488). 
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number basic errors associated with predictive procedures. The success 
failure the latter the ultimate test the validity constitutive model. 
The authors appear restrict their treatment the rather specialized axisym- 
metric condition, stating 628, that ‘‘attempt account for the influence 
would unnecessary complication this stage the development 
the 

The axisymmetric condition, although convenient for making measurements, 
poor representation common types deformation such plane strain. 
this writer’s experience that satisfactory correction factors based 
comparison between the extended forms the axisymmetric and the three- 
dimensional tensors can result greater accuracy. 

The authors’ assumption (found 615), that was independent 
both deviatoric and volumetric stress’’ surely invalid, bearing mind that 
essentially relationship between stress and strain. There wealth 
evidence indicate otherwise, Refs. 24, 25, 27, 41, 43, and inclusive 


consolidation 
specimens 
radial drainage 

constant throughout test 


FIG. 21.—One-Dimensional Consolidation (Radial Drainage) Showing Variation 
Lateral Stress (Hanrahan, 40) 


(also Ref. from the authors’ paper). Admittedly, there are special occasions 
when does appear (42) that independent load. these circumstances 
usually found that after certain loading period, the lowered moisture-content 
(and, thereby increased stiffness) the soil will compensate for the increased 
stress. This balance can, however, upset slight changes conditions 
test, such size specimen, length drainage path, stress history, duration 
loading time, etc. 

The authors’ statement, cited 613, that ‘‘one-dimensional compression 
represents deformation constant deviatoric stress level, assumed 
that K,, the coefficient lateral earth pressure rest, 
serious error. Thus 


K,= (by definition) 


peat initial 4857 
\ WwW s 140% | 
0-8 6, s 10 psi. (69 kPa) 
organic silt initial 46% 
0-2 
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Transposing 


Differentiating 


Thus, variable, which varies with time, depth, and rate which 
proportional the rate decrease pore pressure. This variation calls 


into question the validity the Terzaghi equation one-dimensional consolida- 
tion 


Clearly, the Biot and Terzaghi equations can only reconciled which 
also been shown experimentally, Fig. (32,40) using special (radially-drained) 
triaxial consolidation test. 

major criticism the paper that only the axial strain appears have 
been considered the formulation the model. attention has been given 
the intermediate and minor principal strains. Under certain important conditions 
restrained deformation, such plane strain and one-dimensional consolidation, 
the interactive effect the components volumetric and deviatoric strains, 
exercises fundamental influence rate porepressure reduction and settle- 
ment. For example, one-dimensional consolidation, the horizontal component 
the deviatoric strain constrained equal and opposite the corresponding 
component the volumetric strain, effect the required condition 
zero lateral strain. Similarly with plane strain, there must corresponding 
equalization the intermediate strain components. both cases, has been 
shown (44) that continuous alteration the appropriate stress (i.e., 
one-dimensional, plane strain) necessary achieve the desired condition 
zero net strain. Furthermore, knowledge the complete stress system 
essential, all times, the strain the other principal directions 
properly evaluated. not clear how such evaluation possible with 
the authors’ model. 


24. Adams, T., Comparison Field and Laboratory Measurements Peat,’’ Proceed- 
ings, 9th Muskeg Conference Peat, National Research Council, 1964. 

25. Berry, L., and Poskitt, J., Consolidation Geotechnique, 22:1, 
London, England, 1972, pp. 27-52. 


(13) 
she Biot (26) equation, assuming highly idealized soil properties, 
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Closure Edward Kavazanjian, ASCE 
and James ASCE 


The writers thank Hanrahan for his penetrating discussion which serves 
focus attention several aspects their model. His work this area will 
valuable for further refinement the model, which now progress. 
The governing equations the writers’ and the discusser’s models are indeed 
somewhat similar. There does, however, appear fundamental difference 
between our postulate immediate and delayed contributions both the 
volumetric and deviatoric components deformation and Hanrahan’s assumption 
one deformation component which can develop only under drained conditions, 
and second component which can develop under both drained and undrained 
conditions. 

The writers wish reemphasize that the work presented this paper 
represented not final ‘‘general constitutive model which can adequately 
predict the time-dependent deformation cohesive soil under three-dimensional 
systems but the first generation the development such model. 
Two types assumptions were made: (1) Those simplify the numerical 
formulation this first generation model; and (2) those which would minimize 
the number parameters required accurately describe soil behavior. The 
formulation has been based philosophy that useful model will one 
that simple possible and uses well-known and easily measured parameters. 

The assumption that the coefficient secondary compression, constant 
and independent both volumetric and deviatoric stress assumption 
the latter type. Certainly one could introduce variety special functions 
closely describe the variation with both volumetric and deviatoric 
stress based upon the results battery laboratory tests. believed 
that the error introduced assuming constant would not usually 
engineering significance and overshadowed the benefit derived from the 
simplicity determining the model parameters for engineering analysis. 
this judgement proves error, the model flexible enough allow 
for insertion stress time-dependent relationship. 

Specimens only one size were used the testing program reported 
the paper. Tests different size specimens are indeed desirable when testing 
the validity the model. However, the use standard size specimens for 
determination soil parameters is, within the context the model, the most 
desirable procedure, because the amount back calculation and curve fitting 
required determine the model parameters minimized. Subsequent testing 
predictions for soil specimens different size order. 

The writers’ statement that ‘‘one-dimensional compression represents de- 
formation constant deviatoric stress level, assumed that K,, the coefficient 
lateral earth pressure rest, constant’’ should amended say 
secondary compression. thank the discusser for 
pointing out this error. Furthermore, are quite aware that even during 
secondary compression the assumption constant not generally valid, 
and have used our model (9) predict the variation with time during 


Prof. Civ. Engrg. Stanford Univ., Stanford, Calif. 94305. 
and Chmn., Dept. Civ. Engrg., Univ. California, Berkeley, Calif. 94720. 
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EDGAR KAOLIN 


PREDICTION 


FIG. 22.—Comparison Observed and Predicted Effective Stress Paths for Edgar 
Plastic Kaolin 
EDGAR PLASTIC KAOLIN 


PREDICTION 
TEST RESULTS 


DEVIATOR STRESS PSI 


STRAIN DIRECTION, PERCENT 


FIG. 23.—Comparison Predicted and Observed Stress-Strain Curves for Edgar 
Plastic Kaolin 
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secondary compression observed Lacerda (47) for San Francisco Bay mud. 

Hanrahan states that major criticism the paper that only axial strain 
appears have been considered the formulation the Predictions 
both volumetric and axial strains were presented the paper. predictions 
radial strain were presented because for the case triaxial compression 
the radial strains can easily derived from the axial and volume strains and 
thus present them would redundant. Triaxial stress states were used 
simplify the mathematical formulation this first generation model and facilitate 
validation testing. the time since the work presented the paper was completed 
the model has been extended consider arbitrary three-dimensional stress states. 
Figs. and show comparisons between the observed soil behavior true 
triaxial tests done elsewhere and behavior predicted under conditions 
using the immediate component the three-dimensional model. Work the 
model continuing with the immediate goal installing the current version 
the constitutive model within the framework finite element program 
order predict the load-deformation behavior full size engineering 
structures. 


46. Kavazanjian, Jr., E., Mitchell, K., and Bonaparte, R., Predic- 
tions Using General Phenomenological presented the May 28-30, 1980, 
NSF/NSERC Symposium Constructive Relationships for Soils, held McGill 
University, Montreal, Canada. 

Lacerda, J., Relaxation and Creep Effects the Deformation 
thesis presented the University California, Berkeley, Calif., 1976, partial 
fulfillment the requirements for the degree Doctor Philosophy. 


Mountain 
Prosect: 


The authors are commended for presenting the performance records 
for the subject project. Actual performance records are most helpful engineers 
engaged designing and constructing similar projects. 

The open joints encountered the drilling and grouting program indicate 
again the need for the design engineers monitor this work. The close-order 
drilling for the grouting work should always expected reveal defects 
the foundation that were not found the more widely-spaced, design-stage 


“June, 1980, William Swiger, Philip Wild, and Thomas Lamb (Proc. Paper 
15506). 
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borings. The decision relocate portion the west dike excellent 
example the design continuing through construction, and the extra costs 
involved avoid future problems associated with the open joints are undoubtedly 
fully justified. Even the unfilled voids the open joints could have been 
filled with nonshrink grout, the high gradients across the grout 
curtain could have caused the sand remaining the joints removed 
ever enlarging pipes. 

The potential for removal soil fines from open joints rock beneath 
upper reservoir embankments similar that for abutments valley bottom 
embankments. both cases, there the potential for full reservoir head 
develop the reservoir side and for ground water appreciable 
depth the downstream side. This situation ready made for the under-seepage 
transport the soil fines unfilled voids lower level carry them 
out springs which develop downstream. Detection developing pipe 
virtually impossible. The grain-by-grain removal starting from the downstream 
side unlikely cause detectable changes the piezometer readings until 
the pipe breaks through the reservoir. The first warning that pipe has 
developed likely the breakout large spring downstream. Unless 
the reservoir can drawn down very quickly, failure the embankment 
apt occur. The rapidity the failure would depend the foundation conditions 
and the composition the embankment. concrete-faced, coarse rockfill 
embankment would less affected pipe the foundation than would 
all earth embankment central-core rockfill embankment. 

The paper makes reference seepage losses upper reservoirs several 
other pumped storage projects. The state-of-the-art for upper reservoirs would 
considerably enhanced the engineers and owners all upper reservoirs 
would present papers detailing the geologic conditions, the design, the construc- 
tion, the operational problems, the remedial measures, and the lessons learned. 

For comparing settlement records with other projects, would helpful 
the authors their closure would present density and gradation data for 
the core, transitions, and rockfill zones. 


Closure William Philip Fellows ASCE, 
and Thomas Lamb,’ ASCE 


wish acknowledge with sincere thanks the thoughtful discussion 
Fetzer. are completely agreement with his statements regarding the need 
for close and continued supervision foundation preparation and grouting for 
dams and with his comments regarding subsurface erosion. Subsurface erosion 
and internal migration are easily overlooked design and have been the source 
significant problems several recent projects. 

Fetzer requested additional information density and gradation for the 
materials the upper dam. The core was obtained from local small drumlin. 


Pres., Stone Webster Engrg. Corp., 245 Summer Street, Boston, Mass. 02107. 
Pres., Stone Webster Engrg. Corp., 245 Summer Street, Boston, Mass. 02107. 
Soils Engr. Stone Webster Engrg. Corp., 245 Summer Street, Boston, Mass. 
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typical New England till grading from cobbles and gravels sands, 
silts, and slight amounts clay. Specifications for the material required 
minimum 60% pass screen and not less than 18% the fraction 
passing the No. screen pass No. 200 screen. Maximum particle size 
permitted the core was in. The material was very well graded within these 
limits. Tests taken the borrow area indicated that fines passing the No. 200 
screen varied from 18%-21% the fraction passing the No. Compaction 
was 50-ton rubber tired rollers. The material was placed in. lifts. 
Moisture content varied from optimum approx wet optimum. The 
required density was not less than 95% standard maximum density, ASTM 
D-698 (American Society for Testing and Materials). The average field density 
obtained was 98.5% for unit density 126.4 

Filter material was obtained crushing selected rock excavated from the 
intake channel. The material was graded from in. fines with not more 
than passing No. 200 screen. Filters were placed lifts and compacted 
minimum 95% standard. Average field density was approx 98%, 
with unit weight 130 

There transition zone fine rock grading from in. in. 
maximum size immediately against the filter. This was compacted 2-ft lifts 
10-ton vibratory roller. The remainder the shells was placed 4-ft 
lifts and compacted three passes 10-ton vibratory roller. Gradation 
was quarry run rock 36-in. maximum size with oversize material placed 


the face the Field density tests were not made the rock filled 
zones. 


Errata.—The following corrections should made the original paper: 


Page 677, Paragraph line Should read the main cavern, tailrace 
tunnel, instead the main cavern, tailrace, tunnel 

Page 678, line Should read fill was placed 4-ft (1.2-m) thick 
instead fill was placed 3-ft (0.9-m) thick 
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